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I.  Abstract 


The  goal  of  this  research  program  has  been  to  explore  and  understand  the 
interaction  of  electromagnetic  radiation  with  the  low  lying  excitation  spectra 
in  condensed  matter  and  on  clean  metal  surfaces.  Both  coherent  and  incoherent 
optical,  infrared,  and  far  infrared  techniques  have  been  required  to  probe  the 
dynamics  of  these  many  body  systems  both  in  and  out  of  thermal  equilibrium. 
The  formation  of  persistent  spectral  holes  is  one  new  aspect  of  solid  state 
spectroscopy.  Another  is  the  discovery  of  large  far  infrared  nonlinearities  in 
high  mobility  two  dimensional  electron  gas  semiconducting  systems. 
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II.  Objectives 

Novel  high  resolution  electromagnetic  wave  techniques  have  been  used  in 
the  optical,  infrared  and  far  infrared  spectral  regions  to  explore  the  electronic 
states  at  metal/dielectric  interfaces.  Because  infrared  surface  plasmons  on 
metal  surfaces  propagate  for  many  wavelengths,  a  measurement  of  the 
transmission  of  these  surface  excitations  has  proven  to  be  a  sensitive  probe  of 
the  surface  itself.  Both  broadband  and  single  frequency  generation  techniques 
have  been  developed.  Reconstructed  surfaces  as  well  as  surfaces  covered  with 
a  chemisorbed  atomic  monolayer  or  a  thin  dielectric  or  molecular  film  have 
been  investigated  with  these  new  methods. 

Persistent  IR  spectral  holes  provide  a  new  high  resolution  technique  with 
which  to  explore  the  properties  of  defects  in  solids.  The  inhomogeneously 
broadened  absorption  lines  associated  with  the  electronic  states  in 
semiconductors  are  in  general  not  coincident  with  the  frequencies  of  ir-line 
tunable  gas  laser;  consequently,  the  high  power  techniques  of  incoherent 
saturation  and  transient  hole  burning  cannot  be  applied  to  identify  the  dominant 
relaxation  process  in  most  electronic-defect-lattice  combinations.  With  low 
power  tunable  diodes  two  techniques  are  potentially  available:  electronic 
fluorescence,  which  gives  the  excited  state  lifetime;  and  persistent  hole 
burning  which  not  only  provides  a  high  resolution  probe  of  the  defect  but  also 
identifies  the  excited  state  dephasing  time  of  the  electronic  state. 

Measurements  have  been  made  on  the  far  infrared  nonlinear  properties  of 
high  mobility  2-DEG  samples.  One  goal  has  been  to  drive  the  collective  2-DEG 
electronic  mode  far  from  equilibrium  with  high  intensity  radiation  tuned  to  the 
natural  quantum  transition  frequencies.  Large  changes  are  observed.  Very 
little  other  work  has  been  carried  out  on  the  nonlinear  high  frequency 
properties  of  2-DEG  systems  and  it  represents  an  exciting  area  for  future 
studies. 


III.  Accomplishments 


A.  Introduction 

This  program  started  with  the  application  of  three  different  IR  techniques 
to  probe  the  electronic  and  vibrational  properties  of  deep  level  defects  in 
semiconductors,  namely,  fourier  transform,  diode  laser,  and  saturation 
spectroscopy.  The  overall  theme  of  this  program  has  been  to  investigate  the 
stability  and  dynamics  of  specific  lattice-defect  systems  by  probing  the 
energy  level  spectra  using  these  modern  high  resolution  spectroscopic 
techniques.  Although  the  spectra  for  most  of  the  deep-level  defects  in  Si  were 
investigated  initially  in  the  1960's(1),  the  motivation  for  our  work  came  from  a 
number  of  studies  in  the  1970's  and  early  1980's  on  related  solid  state  systems 
which  showed  that  persistent  spectroscopic  changes  could  often  be  induced  in 
defect  systems  (2'6).  The  persistent  effect  is  usually  characterized  by  a  very 
small  quantum  efficiency,  i.e.,  it  takes  .JO6  transitions  before  the  defect 

becomes  trapped  in  another  state  which  is  not  resonant  with  the  laser 
frequency;  on  the  other  hand,  because  of  the  high  brightness  of  the  cw  probe 
and  the  integrating  nature  of  the  effect,  appreciable  conversion  can  be  observed 
on  minute  time  scales.  Dynamical  mechanisms  which  produce  this  optical 
effect  continue  to  be  identified  for  an  ever  increasing  number  of  inorganic  and 
organic  crystals  and  glasses  (7*9).  Persistent  IR  spectral  hole  burning  also  has 
been  observed  for  molecular  defects  in  noble  gas  and  alkali  halide  crystals 
where  no  electronic  excitation  is  involved  and  only  vibrational  degrees  of 
freedom  are  photoexcited  with  infrared  radiation  (10). 

With  complex  defects  a  variety  of  inequivalent  defect  orientations  may 
provide  the  necessary  low  temperature  traps  preventing  the  return  of  the 
defect-lattice  system  to  the  equilibrium  ground  state  but  there  is  now 
concrete  evidence  that  e’ven  a  point  defect  in  an  ionic  crystal  can  have  more 
than  one  low  energy  elastic  configuration  available  to  it (11,12).  Perhaps  some 
of  the  diverse  results  found  in  the  persistent  hole  experiments  with  other 
kinds  of  solids  are  a  consequence  of  these  little  understood  low  lying  elastic 
states  which  also  would  inhibit  a  direct  return  to  the  ground  state 
configuration. 

For  covalently  bonded  crystals  only  one  observation  of  persistent  spectral 
hole  burning  in  diamond  had  been  reported  (12).  The  zero  phonon  lines  of  four 
different  defect  centers  all  showed  this  effect  indicating  that  the  phenomenon 
was  a  rather  general  one  for  this  host.  Silicon  is  one  of  the  best  controlled 
hosts  for  defect-lattice  work  hence  it  provides  an  interesting  medium  in  which 
to  study  the  stability  and  dynamics  of  defect-lattice  system.  Since  persistent 
IR  spectral  holes  (PIRSH's)  provide  a  new  way  of  probing  these  states  it  is 
important  to  determine  the  conditions  on  PIRSH  production  in  the  SI  host.  In  a 
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general  sense,  one  expects  that  three  requirements  must  be  met.  First,  there 
must  be  several  ground  state  configurations  of  the  total  system,  and  the  IR 
absorption  energies  from  these  ground  states  must  differ  by  more  than  the 
laser  linewidth.  Second,  there  must  exist  an  IR  pumping  pathway  that  connects 
these  nearly  degenerate  "ground"  state  configurations.  Finally,  the  relaxation 
among  the  "ground"  states  must  be  slower  than  the  excited  state  lifetime.  If 
all  of  these  conditions  are  met  then  persistent  spectral  holes  may  result. 

We  have  succeeded  in  observing  persistent  spectral  holes  for  a  few 
Si-defect  systems  but  unlike  diamond  most  combinations  do  not  show  long  time 
hole  burning  effects.  To  understand  the  differences  between  these  two  hosts  it 
has  been  necessary  to  fabricate  a  large  number  of  defect  systems.  Fourier 
transform  IR  spectroscopy  has  not  only  provided  the  first  high  resolution  data 
on  many  of  these  centers  but  also  operated  as  a  broad  band  probe  beam  in  a  new 
class  of  pump-probe  experiments  on  defect  centers.  By  combining  IR  laser  and 
broad  band  techniques,  we  made  a  number  of  discoveries  including  the  first 
observation  of  motional  tunneling  at  deep  level  electronic  defects  in  elemental 
semiconductors. 

B.  IR  Saturation  Measurements  on  Peep  Levels 

Although  high  power  saturation  studies  of  shallow  levels  in 
semiconductors  show  that  impurity  lifetimes  in  the  range  from  1  ns  to  1  jis  can 
be  obtained  at  low  temperatures  (U*17\  before  our  work  no  similar 
measurements  on  deep  levels  had  been  reported.  The  main  reason  for  this  lack 
of  data  has  been  that  high  power  tunable  lasers  are  not  available  in  the  IR 
region  so  that  even  a  simple  saturation  measurement  requires  that  the 
absorption  spectrum  of  the  defect-lattice  system  accidentally  be  coincidence 
with  the  fixed  frequencies  available  from  a  gas  laser  ( such  as  C02).  We  have 

found  two  absorption  Iin6s  which  satisfy  this  requirement:  one  is  the  1s(A1)-> 

1  s(Tt )  transition  of  the  Se  donor  in  AISb(18)  and  the  other  is  the  IsfA^-^sO^) 

transition  of  H:S  donor  in  Si.  We  review  now  the  important  findings(19)  for  the 
Se  donor  in  AlSb. 

At  low  temperature  the  5.6  cm'1  wide  1s(A1)->1s(T1)  transition  is 
straddled  by  two  C02  laser  lines.  The  line  breadth  can  be  interpreted  in  one  of 

two  ways:  either  the  line  is  homogeneously  broadened  and  the  low  temperature 
linewidth  is  a  measure  of  the  intrinsic  multiphonon  decay  lifetime  or  the  line 
is  inhomogeneously  broadened  due  to  crystal  strains  so  the  low  temperature 
width  obtained  from  linear  spectroscopy  provides  no  intrinsic  information 
about  the  lifetime.  For  the  former  case  the  laser  intensity  where  saturation 
should  commence  can  be  estimated  since  both  the  absorption  cross  section  and 
the  linewidth  are  known.  The  calculated  number,  3.7  MW/cm2,  represents  the 
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maximum  possible  value  for  this  system.  The  saturation  parameter  must  be 
smaller  in  the  case  of  inhomogeneous  broadening  since  the  excited  state 
lifetime  is  necessarily  longer.  An  experimental  measurement  permits  one  to 
distinguish  between  these  two  possibilities.  We  find  that  the  IsO^)  lifetime 

is  very  short  between  1  and  5  ps  which  is  about  three  orders  of  magnitude 
smaller  than  found  for  shallow  levels  in  three-five  compounds.  This  result 
shows  for  the  first  time  that  the  increase  in  the  electron-phonon  coupling 
constant  for  these  localized  states  more  than  balances  against  the  highly 
forbidden  multiphonon  decay  required  to  relax  them. 

Similar  results  have  been  obtained  for  the  H:S  donor  in  Si. 

C.  Tunneling  at  Deep  Level  Acceptor  Complexes  in  Si 

Some  unusual  properties  of  hydrogen  or  lithium-related  shallow  centers  in 
Ge,  which  have  been  identified  by  FIR  photo-electric  spectroscopy,  were 
originally  ascribed  to  the  dynamic  tunneling  of  the  light  atom  among  its 
equivalent  sites  in  the  lattice(20'22).  Recently  it  has  been  shown  that  all  of 
these  acceptor  centers  with  the  possible  exclusion  of  the  Cu-H2  do  not  behave 

as  tunneling  centers(23).  For  the  donors  other  models  exist(24)  to  explain  their 
unusual  stress  dependent  IR  spectra  and  measurements  are  in  progress  to 
determine  which  model  is  correct(23).  The  current  status  is  that  none  of  the 
shallow  centers  unambiguously  display  dynamic  tunneling  behavior. 

With  IR/FIR  spectroscopy  we  have  carried  out  a  systematic  study  in  Si  on 
a  series  of  isomorphic  defects{25,26)  Be-H,  Be-D  and  Be-Li.  Dynamical 
tunneling  is  observed  for  the  acceptor  complexes  composed  of  substitutional  Be 
plus  H  or  D. 

The  acceptor  A(Li,  Be),  an  example  of  a  non-tunneling  center,  produces  an 
acceptor  hydrogen  series(27)  around  800  cm'1.  Temperature  dependent 
measurements  show  tha't  the  rg  ground  state  is  split  into  two  Kramers  doublets 

(  aE  „  1 1 .3  cm'1 )  due  to  the  symmetry  lowering  (Td  ->  C3v )  produced  by  the 

Li+  ion3.  The  hole  absorption  spectra  of  the  A(H,  Be)  and  A(D,  Be)  centers  are 
different  and  show  three  new  distinguishing  features.  1)  When  H  ->  D  the  eight 
transitions  all  display  a  large  isotope  shift  of  7.8  cm*1  to  higher  frequencies. 

2)  A  weak  high  frequency  replica  of  these  transitions  is  observed,  even  at  the 
lowest  temperature  of  1 .7  K,  with  a  shift  of  38.8  cm'1  for  A(H,  Be)  and  16.2 
cm'1  for  A(D,  Be).  3)  With  increasing  temperature  two  additional  replicas 
appear  in  sequence  on  the  high  frequency  side  of  the  main  lines.  The  spectra 
for  the  A(D  ,Be)  center  at  a  number  of  temperatures  are  shown  in  Fig.  2  of  Ref. 
() .  I  R  spectra  for  A(H  ,Be)  show  quite  similar  temperature  dependences  but 
the  spectral  changes  occur  at  somewhat  higher  temperatures. 

These  high  resolution  temperature  dependent  measurements  of  the 


infrared  and  far  infrared  spectra  associated  with  the  acceptor  complexes  Be-H 
and  Be-D  in  Si  provide  a  direct  identification  of  motional  tunneling  at 
electronic  defects  in  semiconductors.  The  large  isotope  shift  of  7.8  cm'1  in  the 
IR  demonstrates  that  the  dynamical  properties  of  hydrogen  play  an  important 
role  in  the  determination  of  the  electronic  energies.  The  observation  of  two 
sets  of  optical  transitions  from  the  ground  to  the  excited  states  with  different 
tunneling  configurations  provides  the  first  direct  determination  of  the 
tunneling  parameter.  At  least  three  levels  of  the  ground  state  manifold  have 
been  identified  for  A(H,Be)  and  A(D,Be)  from  the  temperature  dependent  studies 
of  IR  and  FIR  spectra  which  should  be  compared  to  the  theoretical  predictions 
of  five  manifolds  with  dynamic  tunneling  and  two  without.  Our  discovery(28) 
that  large  isotope  effects  can  be  observed  for  deep  level  tunneling  centers  in  Si 
(density  _  1016  cm'3  )  opens  up  the  experimental  possibilities  so  that  now  a 

variety  of  spectroscopic  tools  can  be  employed  to  probe  the  dynamical 
behavior  of  these  unusual  electronic  defects. 

D.  Nonlinear  Properties  of  Cyclotron  Resonance  in  a  Two 
Dimensional  Electron  Gas 

Although  the  d.c.  electric  field  dependence  of  the  nonlinear  mobility  of 
electrons  in  a  two  dimensional  electron  gas  (2-DEG)  structure  has  been  well 
characterized  because  of  applications  to  fast  field  effect  transistors(29,30), 
information  on  the  nonlinear  mm  wave  response  of  these  high  mobility  systems 
was  lacking  until  the  work  initiated  by  Schlesinger  and  Sievers(31).  The  inital 
measurements  were  made  while  one  of  us  (AJS)  was  on  sabbatic  leave  at  the 
high  intensity  near  mm  wave  laboratory  of  Dr.  F.  Keilmann  at  the  Max  Plank 
Institute  in  Stuttgart.  This  opportunity  was  used  both  to  learn  about  new  high 
power  techniques  in  this  frequency  region  and  to  initiate  an  investigation  into 
the  nonlinear  mm  wave  properties  of  semiconductor  interfaces  which  support  a 
2-DEG.  Our  preliminary  measurements  on  the  far  infrared  intensity  dependence 
of  the  magneto-transmissivity  of  a  2-DEG  at  a  GaAs/(AI0  3Ga0  7)As  interface 

has  been  described  in  Ref.  (31).  The  data  demonstrate  that  high  intensity 
cyclotron  resonance  provides  a  new  method  for  identifying  and  studying 
specific  hot  electron  effects  in  these  systems. 
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The  possibility  of  observing  quantum  sue  effects  in  very  small  metal  particles  by  absorption  spec* 
troscopy  is  reexamined.  According  to  the  Gor’kov-Eliashberg  theory,  nonquadratic  frequency 
dependence  of  the  absorption  coefficient  persists  at  sufficiently  low  frequencies  even  when  there  is  a 
size  distribution. 


Early  experimental  studies  of  the  absorption  of  far- 
infrared  radiation  by  very  small  metal  particles1,2  were 
motivated  in  part  by  the  possibility  of  observing  quantum 
size  effects  (QSE’s)  predicted  by  the  theory  of  Gor’lcov 
and  Eliashberg3  (GE).  According  to  GE,  the  signature  of 
QSE’s  is  the  presence  of  oscillations  superimposed  on  the 
quadratic  frequency  dependence  of  the  absorption  coeffi¬ 
cient.  No  evidence  for  QSE’s  in  the  far  infrared  has  been 
observed.  The  negative  results  are  consistent  with  the 
corrected4-6  GE  theory  if  experimentally  realizable 
particle-size  distributions  are  taken  into  account.  The  os¬ 
cillations  are  washed  out  by  “destructive  interference.” 

The  purpose  of  this  note  is  to  reexamine  the  possibility 
of  observing  QSE’s  by  far-infrared  or  microwave  spectros¬ 
copy.  The  GE  theory  is  emphasized  over  other  models7 
due  to  the  richness  of  phenomena  it  predicts.  In  particu¬ 
lar,  nonquadratic  frequency  dependence  in  the  low- 
frequency  absorption  coefficient,  which  remains  even  if 
there  is  a  size  distribution,  is  pointed  out  for  the  first 
time. 

Gor'kov  and  Eliashberg  calculated  the  frequency- 

dependent  electric  susceptibility  of  a  small  metal  particle 

by  linear  response  theory.  Their  result  is 

*(<u)=*o+QFU>) ,  (1) 

where  Xo=AkFx1/2(hriaB  is  the  static  susceptibility  and 

n  139  A 

0  HQOir2  kFaB  ' 

Here  x  is  the  particle  diameter,  kF  is  the  Fermi  wave  vec¬ 
tor,  A =m*/m  is  the  specific-heat  effective-mass  ratio,- 
and  aB  is  the  Bohr  radius.  F(<o)=F\  +  /F2  is  the  dimen¬ 
sionless  integral, 

F{w)  =  l{fuo)lD(EF)  £  (3, 

£2-(*o-hO+)2 

where  D(EF)  is  the  density  of  electron  states  at  the  Fermi 
energy  and  R[E)  is  a  two-level  correlation  function. 
Specifically,  R[El—E1)dEldE1  is  the  probability  that  an 
energy  level  of  the  particle  is  within  each  of  the  intervals 
(£i,£i  +d£, )  and  (E1,Ei+dEl),  regardless  of  the  posi¬ 


tions  of  any  other  levels.  GE  borrowed  the  function  R  (£) 
from  the  statistical  theory  of  levels,*  which  was  developed 
by  Wigner,9  Dyson  and  Mehta,10  and  others  for  applica¬ 
tion  to  highly  excited  atomic  nuclei.  In  this  model,  the 
levels  of  the  particle  are  treated  statistically  by  means  of 
an  ensemble  average  over  all  possible  Hamiltonians  con¬ 
sistent  with  general  symmetries  (e.g.,  time-reversal  invari¬ 
ance,  rotational  invariance,  etc.).  There  are  three  cases 
and  hence  three  distinct  ensembles,  called  orthogonal, 
symplectic,  and  unitary  after  the  groups  of  transforma¬ 
tions  under  which  the  ensembles  are  invariant  The  model 
takes  into  account  correlations  between  levels  due  to  level 
repulsion.  The  expressions  for  F(co)  in  GE  are  incorrect, 
but  the  integrals  were  recently  reevaluated.6  The  correct¬ 
ed  results  predict  oscillations  of  the  far-infrared  absorp¬ 
tion  coefficient,  largest  for  the  symplectic  ensemble,  su¬ 
perimposed  on  the  quadratic  background  predicted  for 
electric  dipole  absorption  by  Drude  metal  particles. 

In  the  far  infrared,  both  the  electric  dipole  and  magnet¬ 
ic  dipole  terms  of  the  Mie  series"  make  important  contri¬ 
butions  to  the  absorption  coefficient  of(a  dilute,  spatially 
random  collection  of  small  (x  <  100  A)  metal  particles 
embedded  in  a  transparent  host^Both  the  GE  model  and 
quantum-box  models7  predict  a  quartic  frequency  depen¬ 
dence  for  the  magnetic  dipole  absorption  for  frequencies 
sufficiently  greater  than  the  Kubo12  gap 

6=  IZirflVm'k/x3 , 

in  contrast  to  the  quadratic  frequency  dependence  ob¬ 
tained  for  particles  of  Drude  metal.  Dyson’s  three  ensem¬ 
bles  become  identical  for  fko»  5  if  the  oscillatory  parts 
are  neglected.  Let  f +/r2  and  msmi+'Mj  be  the 
complex  dielectric  function  and  magnetic  permeability  of 
a  particle  and  let  the  corresponding  quantities  for  the  host 
be  unity.  If  the  volume  fraction  of  metal,  /,  in  the  com¬ 
posite  satisfies  /« 1  and  6j«f|,  Mi— Mi«*»  and 
4jtYo»1,  the  low-frequency  absorption  coefficient  in  the 
GE  model  is 


a=  18ir/v 


(1  ,  M2 

(4  irA'o)2  +  9 


(4) 
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where  v=<u/2nr  is  the  frequency  in  wave  numbers.  The 
two  terms  represent  electric  and  magnetic  dipole  absorp¬ 
tion,  respectively.  Note  that  the  GE  model  requires1  that 
a «2vp/x,  where  vF  is  the  Fermi  velocity.  In  the  far  in¬ 
frared 


fe\+^xv)2t 


(5) 


for  small  spheres.11  Thus,  for  the  GE  model14  under  the 
specified  conditions. 


<r 

/ 


139 ir3 


«*. _2  ,  139V3  aV_4 
M/*  V  +  4500  aBK  V 


(6) 


where  k,=#/mc  is  the  Compton  wavelength  of  the  elec¬ 
tron.  The  GE  model  predicts  the  same  frequency  and  size 
dependence  for  the  electric  dipole  term,  but  higher  powers 
of  v  and  x  for  the  magnetic  dipole  term,  compared  with 
the  Drude  model.  Quantum-box  models7  lead  to  the  fre¬ 
quency  and  size  dependence  given  by  Eq.  (6),  but  to  slight¬ 
ly  different  values  for  the  coefficients. 

As  an  example,  consider  Au  particles15,16  with 
kF=  1.21x10*  cm-1  and  A = 1.1.  Equation  (6)  for  the 
GE  model  predicts 


4  =2.02x  10~3—  +  5.68X  HT22*5*4  . 
/  x 

For  particles  of  Drude  Au 

4  =  1.59X10-3—  +2.21  XlO-11*1*2  . 
/  x 


(7) 


(8) 


In  Eqs.  (7)  and  (8),  a  and  v  are  in  cm-1  and  x  is  in  A. 
The  electric  dipole  terms  are  comparable  in  magnitude. 
To  observe  the  quartic  frequency  dependence  of  the  mag¬ 
netic  dipole  term  predicted  by  the  GE  model,  while  satis¬ 
fying  the  assumptions  of  the  model,  the  Au  particles 
should  be  a  few  hundred  angstroms  in  diameter. 

A  quartic  frequency  dependence  of  the  absorption  coef¬ 
ficient  in  the  far  infrared  can  be  interpreted  as  a  quantum 
size  effect.  However,  the  magnetic  dipole  term  was  ob¬ 
tained  by  substitution  without  justification  of  the  dielec¬ 
tric  function  into  the  classical  expression  rather  than  by 
explicit  quantum-mechanical  calculation.17 

If  the  GE  model  in  fact  describes  the  behavior  of  small 
metal  particles,  there  is  a  quantum  size  effect  at  very  low 
frequencies  (v«S)  due  to  level  repulsion  that  does  not 
arise  in  the  box  model.  For  /« 1,  the  electric  dipole  ab¬ 
sorption  coefficient  for  uniformly  distributed  particles  of 


one  size  is 


/ 


IT1 

— Cv3/52  (orthogonal  ensemble) 

-i^-jr5Cv6/65  (symplectic  ensemble) 
2  __ 

— j-Cv4/83  (unitary  ensemble)  , 


(9) 


where 


C  = 


417  A 


1 


50  kFaB  (3+4irX0>2 


(10) 


and  6  is  the  Kubo  gap  in  wave  numbers.  Note  that  each 
ensemble  has  a  different  exponent  for  the  frequency 
dependence  and  that  the  predictions  for  all  three  ensem¬ 
bles  differ  from  the  Drude  model. 

This  effect  is  present  even  when  there  is  a  size  distribu¬ 
tion.  Let  the  subscript  i  denote  the  components  of  a 
discrete  size  distribution.  For  /« 1  and  4vXa  » 3  (to 
make  the  integrals  easy  to  evaluate). 


/v  417 


\Kf  t 


(11) 


with  2///  =/»  xi  tl,e  diameter  of  particles  of  type  i,  and 
Fa  the  result  of  integration  over  one  of  Dyson’s  ensem¬ 
bles.6  For  convenience,  replace  the  sum  by  an  integral 
over  the  log-normal  distribution,  i.e.. 


^  ,  ,exp[--f(ln2o>]  r. 

]-dyap 


JZll 


2(  ln'cr) 


+  3 y 


(12) 


where  y  =ln(x/x),  and  x  and  a  determine  the  mean  size 
and  width.  The  log-normal  distribution  often  provides  a 
good  fit  to  measured  histograms.16  After  integration,  for 
v«8  (here  6  is  the  Kubo  gap,  in  wave  numbers,  for  a 
particle  of  diameter  X),  the  reduced  absorption  coeffi¬ 
cient,  A(z)==axkf//v2,  which  would  be  constant  for 
low-frequency  electric  dipole  absorption  by  well-separated 
particles  of  Drude  metal,  is 


A  (z)s 


vBz  exp[  8(  In2o  )z]  ( orthogonal ) 
■jj-^expJ-^OnV)]  (symplectic) 

/?z2exp[-y-(ln2o)]  (unitary). 


(13) 


where  B=(  139ir1/24)(a*Ac)  and  z=2iri'/S=2irfko/6. 

In  the  Mie  theory,  which  applies  to  the  dilute  limit,  the 
electric  dipole  absorption  coefficient  is 


a=18i7V6o/32 


fjfV 


J  (2€0+Ciy)2  +  «iy 


(14) 


where  e0  is  the  dielectric  constant  of  the  supporting  host 
and 

fy  =* ly  +  —  1  +4ir[Yoy  +CFy(a>)]  . 

Figures  1  through  3  show  log-log  plots  of  the  low- 
frequency  absorption  coefficient  obtained  from  Eq.  (14) 
for  small  Au  particles  in  the  GE  model  when  the  oscilla¬ 
tions  due  to  discrete  levels  are  washed  out  by  a  log-normal 
size  distribution.  Although  the  calculations  were  done 
specifically  for  Au  particles  of  20-A  mean  diameter,  the 
plots  in  dimensionless  units  are  essentially  “universal”  and 
should  apply  to  particles  of  any  metal  with  reasonable 
mean  diameter.  The  reduced  absorption  at  low  frequen¬ 
cies  should  appear  as  a  gap  in  a  linear  plot.  Figure  4 
shows  the  far-infrared  absorption  coefficient  predicted  for 
the  throe  ensembles  as  well  as  for  the  Drude  model  for 
7=20  A  Au  particles.  According  to  this  plot  there  is  a 
change  in  absorption  due  to  discrete  levels  for  a  log¬ 
normal  distribution  with  o=I.2  in  the  far  infrared.  The 
three  ensembles  show  distinct  non- Drude  power  law 


12 


32 


POSSIBILITY  OF  OBSERVING  QUANTUM  SIZE  EFFECTS  . . . 


1953 


FIG.  1.  Reduced  Mie  electric  dipole  absorption  coefficient  vs 
reduced  frequency  predicted  by  the  Gor’kov-Eliashberg  model 
for  polydisperse  Au  particles  with  /cf=  1.21x10*  cm-1  and 
A  =  1.1  for  the  orthogonal  ensemble.  The  supporting  dielectric 
is  vacuum  («o=  1.0).  The  log-normal  size  distribution  is  used, 
with* =20  A. 


FIG.  2.  Reduced  Mie  electric  dipole  absorption  coefficient  vs 
reduced  frequency  predicted  by  the  Gor’kov-Eliashberg  model 
for  polydisperse  Au  particles  with  kf=  1.21x10*  cm-1  and 
A  =  1.1  for  the  symplectic  ensemble.  The  supporting  dielectric 
is  vacuum  («o=l.0).  The  log-normal  size  distribution  is  used, 
with  X=20  A. 


FIG.  3.  Reduced  Mie  electric  dipole  absorption  coefficient  vs 
reduced  frequency  predicted  by  the  Gor’kov-Eliashberg  model 
for  polydisperse  Au  particles  with  kf=  1.21X10*  cm-1  and 
A  =1.1  for  the  unitary  ensemble.  The  supporting  dielectric  is 
vacuum  (e0=  1  -0).  The  log-normal  size  distribution  is  used,  with 
X=20  A. 


Frequency  (cm*1 ) 


FIG.  4.  Frequency  dependence  of  the  electric  dipole  absorp¬ 
tion  coefficient  for  a  log-normal  distribution  of  Au  particles 
predicted  by  the  GE  model  with  k,=  l. 21x10*  cm-1,  A  =  l.l, 
/* 0.01,  X  =  20  A,  and  <7=  1.2.  Also  shown  is  the  prediction  for 
particles  of  Drude  Au  with  boundary  scattering  (r=*/2ty). 
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behavior  for  v  <  5  and  quadratic  frequency  dependence  for 
v>6.  The  small  discrepancies  between  these  curves  and 
the  predictions  of  Eqs.  (7)  and  (13)  are  due  to  the  simpli¬ 
fying  assumptions  (4jrf0»3,  2f0+«t j  »*2j)  made  to  ob¬ 
tain  the  asymptotes  analytically.  The  “gap”  due  to 
discrete  levels  (Fig.  4)  is  not  abrupt  because  the  low- 
frequency  absorption  follows  a  power  law  rather  than  an 
exponential,  and  the  size  distribution  introduces  broaden¬ 
ing.  Surprisingly,  a  becomes  constant  (not  shown  in  the 
figure),  rather  than  quadratic,  in  frequency  according 
to  Eq.  (14)  for  sufficiently  high  frequencies 
(v/8»-fifkfX2 3).  However,  the  assumptions  of  the  GE 
model  break  down  in  this  region,  so  this  prediction  cannot 
be  taken  seriously.  #  _ 

For  J=20  A  and  x  =  50  A  Au  particles,  5=218  and 
13.9  cm-1,  respectively,  so  it  is  advantageous  to  use  very 
small  particles  to  search  for  nonquadratic  frequi ,  cy 
dependence.  The  effect  would  be  difficult  to  pin  down  for 
any  x,  since  the  high  power  laws  imply  very  small  absorp¬ 
tion  coefficients  for  small  v.  Also,  to  extract  the  exponent 
from  data,  the  baseline  for  a  must  be  known  accurately. 
Figure  4  shows  that  the  predictions  for  the  three  ensem¬ 
bles  are  not  much  different  in  the  far  infrared  when  linear 
axes  are  used.  The  most  convincing  verification  of  the 
GE  theory  would  be  to  transform  one  ensemble  into 
another  and  observe  the  change  in  frequency  dependence. 
For  example,  the  symplectic  ensemble  (a~v6 7 * * 0)  applies  to 
a  metal  with  strong  spin-orbit  coupling,  but  in  a  strong 
magnetic  field  the  unitary  ensemb *le  (a—v4)  is  appropri¬ 
ate.  The  transformation  could  be  observed  by  dividing  the 
spectra  obtained  for  a  specific  sample  by  Fourier- 
transform  spectroscopy  (FTS)  with  and  without  an  ap¬ 


plied  magnetic  field.  The  sensitivity  of  the  FTS  technique 
is  maximized  for  this  measurement  because  no  properties 
of  the  sample,  including  its  position,  are  changed,  except 
for  the  applied  field.  If  additional  sensitivity  is  required, 
a  strong  coherent  source  with  discrete  lines,  such  as  a  far- 
infrared  laser,  might  be  used.  Samples  with  dilute  con¬ 
centrations  of  well-separated  particles  are  required  to 
avoid  interactions  between  particles,  including  the 
enhanced  far-infrared  absorption1’2, 19,20  associated  with 
clusters  of  particles,  which  might  overwhelm  the  small  ef¬ 
fect  under  consideration  here.  A  small-particle  composite 
with  well-characterized  morphology  that  can  be  studied 
by  far-infrared  techniques  was  recently  developed,20  but  a 
search  for  quantum  size  effects  using  this  material  would 
be  difficult  due  to  absorption  by  the  host. 

We  conclude  the  following. 

(1)  Oscillations  in  the  frequency  dependence  of  the  far- 
infrared  absorption  coefficient  a  predicted  by  the  GE 
model  are  unobservable  with  presently  available  samples 
due  to  the  width  of  the  size  distribution.6 

(2)  Quartic  frequency  dependence  of  a  in  the  far  in¬ 
frared,  if  observed,  can  be  interpreted  as  a  quantum  size 
effect. 

(3)  The  GE  model  can  be  tested  by  looking  for  nonqua¬ 
dratic  frequency  dependence  (i.e.,  an  energy  gap)  in  a  at 
low  frequencies  for  very  small  particles. 
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Two-dimensional  electron  gas  in  ln0.53  Ga0  47  As/lnP  heterojunctions  grown 
by  atmospheric  pressure  metalorganlc  chemical-vapor  deposition 
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High-quality  modulation  doped  In^si  Ga^  As/lnP  heterostructures  have  been  grown  by 
atmospheric  pressure  metalorganic  chemical-vapor  deposition  (MOCVD)  using  solid 
trimethylindium  source.  The  two-dimensional  nature  of  electrons  bound  in  the  In^j  Gao.47  As/ 

InP  heterojunctions  is  proved  by  a  Shubnikov-de  Haas  effect  experiment.  Electron  Hall  mobilities 
as  high  as  12000, 83000, 98000,  and  92000  cmW  s  at  300, 77, 40,  and  4.2  K  are  obtained, 
respectively.  The  electron  effective  mass  is  measured  to  be  m  J*  =  0.043  m0  by  cyclotron 
resonance  experiments  on  the  samples  with  two-dimensional  electron  sheet  concentrations  of 
(3.0-3.7)X  10' '/cm2.  From  far-infrared  impurity  absorption  data  the  ionization  energy  of  the 
residual  donors  in  the  MOCVD- grown  In^jj  Gao.47  As  is  determined  to  be  2.95  meV. 


I.  INTRODUCTION 

Two-dimensional  electron  systems  bound  in  semicon¬ 
ductor  surfaces,  heterojunction  interfaces,  or  quantum  wells 
have  been  attracting  great  attention  for  their  unusual  phys¬ 
ical  properties  such  as  quantum  hall  effect,1  and  because  of 
their  promise  for  applications.2  Advanced  silicon  MOS  tech¬ 
nology  has  made  the  electron  system  in  this  structure  avail¬ 
able  for  scientific  investigations,  which  confirmed  for  the 
first  time  the  two-dimensional  nature  of  the  electrons  in  the 
inversion  layer  at  the  surface.3  However,  it  has  been  found 
that  electron  mobility  in  this' structure  is  limited  by  surface 
roughness  scattering,4  which  is  difficult  to  eliminate. 
GaAlAs/GaAs  heterojunction  interfaces  grown  by  molecu¬ 
lar  beam  epitaxy  or  metalorganic  chemical-vapor  deposition 
(MOCVD)  are  abrupt  and  smooth  on  the  monolayer  scale, 
and  mobility  of  the  two-dimensional  electron  gas  (2-DEG)  is 
limited  mainly  by  ionized  impurity  scattering  at  low  tem¬ 
perature.  Hence,  modulation-doped  GaAlAs/GaAs  hetero¬ 
structure  exhibited  4  K  electron  mobilities  as  high  as  2  X 106 
cmVV  s  under  illumination,3  a  result  which  stimulated  in¬ 
tensive  research  on  high-speed  devices.  The  lattice-matched 
I%5i  Gao.47  As/lnP  heterostructure  is  another  interesting 
materials  system.  For  In^  Gao.47  As,  room-temperature 
electron  mobility  as  high  as  12000  cm2/V  s  is  available,  and 
peak  electron-drift  velocity  is  higher  and  negative-differen¬ 
tial  mobility  is  larger  than  GaAs.6  This  makes  InGaAs  espe¬ 
cially  attractive  for  high-speed  devices.  Furthermore,  its 
band  gap  is  wide  enough  for  device  application  at  room  tem¬ 
perature  and  its  optical  properties  are  suitable  for  fabricating 
optoelectronic  devices.  So  far,  the  only  reports  on  2-DEG  in 
InGaAs/InP  heterostructures  grown  by  MOCVD  have 
come  from  the  French  group  at  Thomson-CSF. 7,3  They  used 
a  low-pressure  MOCVD  system  with  liquid  triethylindium 
(TEI)  as  the  indium  source  to  get  sharp  heterointerfaces  and 
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to  avoid  adduct  formation.  Their  highest  77  K  mobility  is 
60000  cm2/V  s9  and  reported  highest  4.2  K  mobility  is  5 1000 
cm2/V  s.7  Besides,  some  groups  are  working  for  high-mobil- 
ity  InGaAs/InP  structures  by  VPE.  For  example,  M.  Ta- 
kikawa  etal.  obtained  2-DEG  mobilities  of 9400, 71200,  and 
106000  cm2/V  s  at  300, 77,  and  4, 2  K,  respectively,  by  chlo¬ 
ride  transport  VPE. 10 

We  report  here  the  first  observation  of  2-DEG  forma¬ 
tion  in  Ino.53  GA0.47  As/lnP  heterojunctions  grown  by  atmo¬ 
spheric  pressure  MOCVD  using  more  stable  solid  trimethy¬ 
lindium  (TMI)  as  the  indium  source.  Various  experimental 
techniques  were  applied  to  characterize  the  2-DEG  at  the 
heterointerface.  Standard  Hall  measurements  were  carried 
out  to  obtain  the  mobility  and  the  sheet-carrier  density  n,  as 
a  function  of  temperature.  Shubnikov-de  Haas  experiments 
were  also  performed  to  determine  independently  the  .sheet 
density  and  to  confirm  the  two-dimensional  nature  of  the 
electron  gas  at  the  heterointerface.  Cyclotron  resonance  ex¬ 
periments  at  various  magnetic  fields  using  far-infrared  Four- 
ier-transform  spectroscopy  determined  the  effective  mass 
for  the  2-DEG.  Unlike  laser  spectroscopy,  Fourier-trans- 
form  spectroscopy  enables  one  to  view  the  entire  resonance 
in  frequency  space  and  determine  the  linewidths  at  different 
fields  unambiguously.  In  addition,  far-infrared  absorption 
spectra  of  the  donor  impurities  in  the  MOCVD-grown 
Ino  ,}  Ga^  As  resulted  in  the  determination  of  the  ioniza¬ 
tion  energy  of  the  residual  donors. 


II.  EXPERIMENTS 

Our  measurements  were  carried  out  on  two  typical  sam¬ 
ples.  Sample  1  is  a  double  heterostructure  (DH)  grown  on  an 
iron-doped  semi-insulating  InP  substrate.  A  0.6^jm-InP 
buffer  layer  is  grown  on  the  substrate,  followed  by  a  i.Q-ftm- 
thick  Ino.9]  Ga<,  «7  As  epilayer  and  a  0.3-jrm-InP  cap  layer. 
The  electron  concentrations  of  the  InP  epilayers  are 
~  1 X  1016  cm-3  and  the  estimated  electron  concentration 
for  the  InGaAs  is  —  2.7  X 1014  cm  “  *.  Sample  2  is  the  same  as 
sample  1  except  without  InP  cap  layer.  The  sample  struc- 
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FIG.  I.  Schematic  diagram  of  the  sample  structures  and  parameters. 


turn  are  shown  schematically  in  Fig.  1.  Growth  of  the  sam¬ 
ples  was  carried  out  by  atmospheric  pressure  MOCVD  with 
a  specially  designed  reactor  configuration  to  enhance 
growth  efficiency  and  interface  sharpness. 1 1  TMI  was  used 
as  the  In  source  with  sublimation  temperature  as  low  as 
12.1  *C.  Attention  was  paid  to  changing  the  gas  flows  when 
growing  interface  regions;  the  details  of  the  growth  proce¬ 
dure  are  described  elsewhere. 12  It  is  worth  pointing  out  that 
our  samples  are  a  little  negatively  mismatched,  with  a  lattice 
mismatch  A a/a  =  —  0.04%  at  room  temperature  for  sam¬ 
ple  1.  This  was  done  to  assure  exact  lattice  matching  at  the 
630X1  growth  temperature  where  excited  atomic  motion 
easily  forms  mismatch  dislocations. 

Hall  measurements  were  carried  out  by  standard  Van 
der  Pauw  technique,  using  an  Oxford  instruments’  CF204 
continuous  flow  cryostat  in  the  temperature  range  from  2.7- 
300  K.  The  sample  was  made  clover-leaf  shaped,  and  tin  dots 
were  alloyed  at  420  *C  for  2  min.  to  make  good  Ohmic  con¬ 
tacts  to  the  2-DEG.  Hall  measurements  were  done  with  a 
magnetic  field  of  2  kG,  keeping  the  sample  in  the  dark. 

Shubnikov-de  Haas  effect  experiments  were  carried  out 
at  1.3— 1.8  K  on  the  samples  immersed  in  pumped  liquid  4He. 
The  sample  was  bar-shaped  with  six  tin  Ohmic  contacts  on 
it  A  constant  current  perpendicular  to  the  magnetic  field, 
was  passed  through  the  sample,  and  magnetoresistance  and 
Hall  resistance  were  measured  through  the  respective  pairs 
of  contacts.  A  superconducting  magnet  provided  a  continu¬ 
ously  variable  magnet  field  from  0  to  7.4  T.  The  angle 
between  the  sample  surface  normal  and  the  magnetic  field 
was  varied  from  0  to  90  *C  by  tilting  the  sample. 

Cyclotron  resonance  experiments  were  made  using  far- 
infrared  Fourier-transform  spectroscopy.  A  lamellar  inter¬ 
ferometer  was  used  in  the  frequency  range  from  5-50  cm-1 
and  a  Michelson  interferometer  was  used  from  50-250 
cm-1.  The  far-infrared  source  was  a  mercury  arc  and  the 
detector  a  germanium  bolometer  cooled  to  0.3  K  by  pumped 
JHe.  Samples  were  held  at  1.3  K  by  immersing  them  in  a 
pumped  liquid  4He  bath.  The  magnetic  field  supplied  by  a 
superconducting  magnet  was  varied  from  0  to  6  T.  The  sam¬ 
ple  transmissivity  was  measured  as  a  function  of  frequency 
for  given  magnetic  fields. 


III.  RESULTS  AND  DISCUSSION 
A.  Low-fMd  transport  proparttm 

For  the  double  heterostructure  sample  there  are  two 
heterojunction  interfaces  which  may  form  two  separated 
quantum-well  channels.  To  determine  the  contribution  of 
each  channel  to  the  total  mobility  and  sheet-carrier  concen¬ 
tration,  the  Hall  measurement  was  repeated  after  selectively 
etching  away  the  cap  InP  layer  with  HQ.  Within  the  — 10% 
experimental  error,  the  Hall  mobilities  and  sheet-carrier 
concentrations  were  the  same  before  and  after  the  etching. 
The  Hall  measurement  results  of  sample  1  are  shown  in  Fig. 
2.  In  the  high-temperature  region  from  100  to  300  K  the 
mobility  decreases  approximately  as  T“ 1 75  because  of  two¬ 
mode  LO  phonon  scattering.  The  mobility  peak  occurs  at 
about  40  K.  The  Hall  data  at  some  important  temperatures 
are  tabulated  in  Table  I.  We  believe  these  are  the  highest 
mobilities  reported  to  date  for  2-DEG  in  Io.3J  Gao  47  As/InP 
grown  by  MOCVD. 

Generally,  the  measured  Hall  data  are  effective  values 
when  a  sample  has  several  conduction  channels.  At  tempera¬ 
tures  below  about  7  K  all  carriers  in  the  bulk  material  are 
frozen  out,  leaving  only  the  2-DEG  still  conducting.  Because 
the  cap  layer  removal  did  not  affect  the  sheet-carrier  concen¬ 
tration  and  mobility,  the  cap  layer  interface  apparently  did 
not  form  a  2-DEG;  this  is  possibly  because  the  surface  bar¬ 
rier  depletion  layer  extended  to  the  interface,  and  prevented 
the  conduction  band  edge  from  being  below  the  Fermi  ener¬ 
gy.  Hence,  only  the  buffer-layer  interface  formed  a  quan¬ 
tum-well  conduction  channel,  and  this  channel  contains  a 
sheet  2-DEG  density  of  n,  —  3.0  X 10' '/cm2.  Our  repeated 
experiments  showed  that  InP  layers  less  than  2.0-2. 5  fim 
thick  grown  with  similar  conditions  on  semi-insulating  InP 
substrates  were  completely  depleted  by  the  surface  barrier 
and  the  substrate-epilayer  interface  barrier.  If  the  surface 
barrier  of  InP  is  —0.25  eV  and  the  interface  barrier  is  —0.6 
eV,13  then  respective  depletion  layer  widths  are  —0.8  and 
—  1 ,2/rm  for  our  doping  level.  This  corroborates  our  conclu¬ 
sion  that,  for  our  samples,  cap  and  buffer  InP  layers  are 
completely  depleted. 

For  a  triangular  well,  the  lowest  quantum  state  energy  is 
(cf.  Stem14) 
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FIG.  2.  Temperature  dependence  of  the  electron  mobility  and  iheet -elec¬ 
tron  density  of  sample  I,  IrwGeo.,  At/lnP  double  hetetoatructure. 
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TABLE  I.  Mobilities  sod  sheet-carrier  concentrations  of  sample  1  at  some 
important  temperatures. 


T|K) 

300 

77 

40 

4.2 

/i(cm  VV  s) 

12000 

S3000 

98000 

92000 

/1,1cm-1) 

4.4x10“ 

3.8X10" 

3.8X10" 

3.0X10" 

E0  =  5(33  in)in  =  48.0  meV.  (1) 

Here  we  used  our  cyclotron  resonance  effective  mass  of  the 
2-DEG  as  =  0.043  and  e,  =  14.  The  Fermi-level 
height  of  16.7  meV  above  E0  is  deduced  by  computing  the 
occupancy  of  this  level  n,  using  the  density  of  states 
D(E)  =  m^/tr/P,  and  neglecting  the  occupation  of  higher 
subbands.  Hence,  the  Fermi  level  is  only  64.7  meV  above  the 
bottom  of  the  well.  This  means  the  interface  quantum  well  is 
shallow,  and  most  of  the  conduction-band  discontinuity  is 
above  the  Fermi  level.  Hence,  the  conduction  band-edge  pic¬ 
ture  of  sample  1  can  be  drawn  as  shown  in  Fig.  3. 

At  high  temperatures,  because  the  dominant  scattering 
is  two-mode  LO  phonon  scattering,  the  mobility  of  both  the 
2-DEG  and  the  InGaAs  bulk  carriers  is  the  same  value  of 
1200  cmVV  sat  300  K  for  the  best  samples.9-12,15  We  assume 
that  at  200  K,  the  2-DEG  and  the  InGaAs  bulk  carriers  still 
have  the  same  mobility.  In  a  two-conduction  channel  case 
the  effective  Hall  concentration  and  mobility  are  expressed 
by'6 

««*  =  (« i  +  «^i)  J/(n +  n^),  (2) 

M*  =  ("i/r?  +  n*i21)/njil  +  nji2),  (3) 

where  the  subscripts  1  and  2  denote  the  bulk  and  the  2-DEG, 
respectively.  If/i,  =  p2  at  200  K,  then  =  n,  +  n2.  From 
the  measured  Hall  data  we  therefore  get  the  Ino,j  Gao  47  As 
bulk-carrier  concentration  (#i£°°)  to  be  (2. 7-3.0) X  10'Vcm5. 
If  the  surface  depletion  barrier  of  InGaAs  is  not  zero,  this 
estimation  is  lower  than  the  real  value.  At  the  temperature 
(—40  K)  of  peak  mobility,  the  donors  in  InGaAs  are  com¬ 
pletely  ionized;  hence,  the  bulk-electron  concentration  is 
nearly  the  same  as  that  at  200  K,  sinqe  the  sheet  density  of 
the  2-DEG  is  a  weak  function  of  temperature,  we  can  take 
n,  =0.8x  10u/cm2,  n2  =  3.0x  10'Vcm2,  and  from  the 
measured  Hall  data  using  Eqs.  (2)  and  (3)  we  get  pb  =  65080 
cmVV  s,  and  p2deo  —  103600  cmVV  s.  Because  the  elec¬ 
trons  in  the  2-DEG  are  transferred  from  the  InP  buffer  layer 
to  the  interface  quantum  well,  the  separation  of  the  2-DEG 


FIO.  3.  Conduction  band-edge  bending  picture  of  sample  I  drawn  schema¬ 
tically  according  to  the  analysis  in  the  text. 


from  the  parent  ions  reduces  the  ionized  impurity  scattering. 
This  enhances  the  low-temperature  mobility  of  the  2-DEG 
to  a  value  much  higher  than  that  of  the  carrier  in  bulk 
I1io.jjGao.47  AS,  as  shown  previously  for  AlGaAa/GaAs 
heterostructures.17 

B.  Shubnlkov-d*  Haas  affect  measurements 

The  existance  of  a  2-DEG  may  be  inferred  from  Fig.  2. 
Direct  proof  is  given  by  Shubnikov-de  Haas  effect  measure¬ 
ments  shown  in  Fig.  4.  These  measurements  were  done  on 
sample  2.  The  sheet-carrier  concentration  of  the  sample  is 
3. 1 X 10* 1  /cm2  at  4.2  K.  by  Hall  measurement.  Figure  4  dem¬ 
onstrates  pronounced  Shubnikov-de  Haas  oscillations  of 
magnetoresistance  and  shows  the  strong  dependence  on 
the  angle  8  between  the  magnetic  field  direction  and  the  lay¬ 
er  normal.  Also  shown  are  corresponding  multistage  Hall 
resistances  pn  for  8  =  55*.  Figure  5  shows  the  reciprocal 
magnetic  field  corresponding  to  the  p„  oscillation  maxima 
as  a  function  of  the  Landau-level  index  N  for  8  =  0  and  55*. 
The  oscillations  are  periodic  in  B  ~ 1  and  follow  the  depen¬ 
dence 

_1_ _  N_  le cos 8 
B~  n,  h 

This  proves  the  two  dimensionality  of  the  electron  gas.  From 
Eq.  (4)  and  Fig.  5  the  sheet  density  of  the  2-DEG  is  calculat- 


FIG.  4.  Results  of  Shubmkov-de  Hus  effect  measurement  on  sample  2,  also 
shown  are  Hall  resistance  p„ . 
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«d  to  be  n,  »  3JX  10“/cm2,  in  food  agreement  with  the 
Hall  data.  In  Fig.  4  there  is  additional  structure  at  low  field 
which  we  (hall  diacuM  elsewhere. 

(X  Cyclotron  resonance  experiments 

1.  Cyclotron  resonance  of  the  2-OEG 

Examples  of  cyclotron  resonance  spectra  for  a  few  fixed 
fields  are  shown  in  Fig.  6.  In  the  figure,  the  main  absorption 
peaks  correspond  to  the  cyclotron  resonance  of  the  2-DEG, 
while  the  smaller  peaks  to  the  right  correspond  to  the  cyclo¬ 
tron  resonance  absorption  of  the  donor-bound  electrons. 
The  frequencies  of  these  peaks  are  plotted  as  a  function  of  the 
magnetic  field  in  Fig.  7.  It  is  clear  that  the  absorption  peak 
frequency  shifts  proportionally  to  the  magnetic  field,  and  no 
deviation  from  linearity  was  observed  in  the  range  of  the 
magnetic  field  used,  indicating  no  coupling  with  LO 
phonon.  Indeed,  in  the  ternary  IncjGA^,  As  the  GaAs- 
like  LO  mode  and  InAs-like  LO-mode  frequencies  are  27S 
and  233  cm  ~ 1 ,  respectively.  '*• 19  The  InP  substrate  TO-mode 
frequency  is  307  cm-1,  so,  it  should  be  possible  to  observe 
the  anticrossing  of  the  2-DEG  cyclotron  resonance  at  233 
cm-1,  only  if  we  raised  the  magnetic  field  to  about  11  T. 
From  the  slope  of  Fig.  7,  using  the  relation  ac  =  eB  /m  *c,  we 
can  extract  a  cyclotron  resonance  effective  mass  of  mg* 
=  0.043  m0  for  the  2-DEG.  The  measurements  on  other 
samples  showed  that  this  value  corresponds  to  2-DEG  densi¬ 
ties  N,  =  (3.0-3.7)X  10ll/cm2,  which  are  measured  by 
Shubnikov-de  Haas  measurements.  The  band  edge  effective 
mass  of  electrons  in  In,,,,  Ga<,  47  As  is  m*  —  0.041  The 

increase  of  the  2-DEG  effective  mass  can  be  attributed  to 
nonparabolicity  of  the  conduction  band.21 

For  the  absorption  spectra  in  Fig.  6,  the  linewidth  &&> 
can  be  expressed  by 

A*u  =  2/rf,  (5) 

if  we  assume  the  line  shape  is  Loren tzian.  Over  the  range  of 
magnetic  fields  used,  the  linewidths  of  resonances  (FWHM) 
are  Aru  —  10-11  cm~\  and  are  nearly  independent  of  the 
field.  Thus,  Eq.  (3)  leads  to  a  relaxation  time 
r#  =  (1.0  ±  0.1)X  10-12  s  for  the  2-DEG.  If  we  denote  the 
low-field  relaxation  time  as  rc,  which  is  deduced  from  the 
relation/*  =  ere/m *,  then  the  ratio  r,/rc  =0.5  for  the  range 
of  magnetic  field  in  Fig.  6. 

2.  Donor-bound  electron  cydotron  resonance 

The  upper  line  in  Fig.  7  represents  donor-bound  elec¬ 
tron  cyclotron  resonance  in  the  InGaAs  bulk  region,  which 
is  caused  by  the  transition  ls-*{2p,m  =  +  1)  of  the  hydro- 
genic  donor-bound  electrons.22,23  In  the  region  where  the 
magnetic  field  satisfies  (1/2)  1iue/R*>l,  where  R  •  is  the 
effective  Rydberg  of  the  donors  in  InGaAs,  this  line  is  ap¬ 
proximately  parallel  to  the  2-DEG  cyclotron  resonance  line. 
If  we  extrapolate  this  impurity  absorption  line  to  zero  field 
and  approximately  take  this  value  as  (3/4)  R  *,  then  the  ioni¬ 
zation  energy  of  the  residual  donor  in  the  MOCVD  grown 
InGaAs  is  R  *  =  2.95  meV.  According  to  the  quasihydro- 
genic  model  of  the  donor  in  semiconductors 

13.6<mVm0)(l/€?).  (6) 
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If  we  take  m*  =  0.43  er  —  14,  then  R  *  =  2.98  meV,  in 
good  agreement  with  the  experimental  result.  However,  the 
residual  donor  impurity  itself  has  not  been  identified  yet. 

IV.  CONCLUSION 

A  high-qualtiy  2-DEG  bound  in  an  Ino  jjGA^o  As/ 
InP  heterojunction  was  obtained  for  the  first  time  by  atmo¬ 
spheric  pressure  MOCVD  using  solid  TMI  as  the  indium 
source.  Low-field  Hall  mobilities  as  high  as  /ijoo  =  12000, 
Hn  —  83000,  /* 42  s  92000  cmVV-s  have  been  obtained 
with  the  maximum  mobility  98000  cm2/V  s  at  around  40  K. 
The  two-dimensional  nature  of  the  electrons  bound  in  the 
buffer  layer  heterointerface  has  been  proved  by  Shubnikov- 
de  Haas  effect  experiments,  and  the  2-DEG  sheet  density 
was  measured  to  be  n,  =  3.2x  lO'Vcm2  for  sample  2,  in 
good  agreement  with  the  results  of  4.2-K  Hall  measurement. 
The  effective  mass  of  the  2-DEG  is  measured  to  be  mg* 
=  0.043  m0  for  2-DEG  densities  (3.0-3.7)X  10n/cm2  by  a 
cyclotron  resonance  experiment  using  far-infrared  Fourier- 
transform  spectroscopy.  The  relaxation  time  of  the  2-DEG 
was  estimated  to  be  (1.0  ±  0.1)X  10~12  s  from  the  cyclotron 
resonance  linewidth.  The  ionization  energy  of  the  residual 
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FIG.  7.  Cyclotron  resonance  absorption  peak  frequency  vt  magnetic  field 
for  lamplc  I. 

donors  in  MOCVD- grown  As  was  measured 

from  the  donor-bound  electron  cyclotron  resonance.  The 
measured  ionization  energy  2.95  meV  is  in  good  agreement 
with  the  theoretical  result  of  the  quasihydrogenic  model. 
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Under  certain  circumstances  the  electric  dipole  interaction  can  reduce  or  eliminate  the  effects  of  spectral 
inhomogeneous  broadening  on  an  infrared  absorption  line  of  a  thin  film.  The  oscillating  electric  dipole  mo¬ 
ment  must  be  oriented  perpendicular  to  the  film  surface,  and  the  dipole  interaction,  as  measured  by  the 
frequency  shift  that  it  induces,  must  be  larger  than  the  width  of  the  inhomogeneous  frequency  distribution. 

An  example  of  this  interaction-induced  line  narrowing  is  presented  for  a  vibrational  mode  of  a  KRe04  Him, 
and  the  relevance  of  this  mechanism  to  two-dimensional  electron-gas  subband  transitions  and  adsorbate  vi¬ 
brational  modes  is  considered. 


Infrared-active  modes  oriented  perpendicular  ( i )  to  the 
plane  of  a  thin  film  may  assume  a  collective  nature  due  to 
the  long-range  (r"2)  electric  dipole  forces  within  the  film. 
For  these  perpendicular  modes,  frequencies  tend  to  be  shift¬ 
ed  upward  and  mode  intensities  do  not  obey  ordinary  sum 
rules.  Analogous  effects  are  found  in  a  number  of  systems 
including  thin  films,1'2  inversion  layers  in  semiconduc- 
tors,1-4  and  adsorbate  systems.5"1  Recent  advances  in  the 
spectroscopy  of  adsorbate  vibrational  modes  at  monolayer 
and  submonolayer  coverages’- 10  as  well  as  the  heightened 
interest  in  two-dimensional  (2D)  electron  systems  enhance 
the  importance  of  understanding  these  effects. 

In  this  paper  we  study  the  v3  vibrational  mode  of  a 
KRe04  thin  film.  This  system  has  a  strong  dipole-dipole  in¬ 
teraction  and  exhibits  unusual  nonlinear  behavior  which  has 
been  discussed  elsewhere."  Here  we  concentrate  on  linear 
properties  and,  in  particular,  show  that  the  perpendicular 
mode  linewidth  is  narrowed  considerably  due  to  the  electric 
dipole-dipole  interaction.  The  possibility  that  this 
interaction-induced  line  narrowing,  for  which  there  is  a 
direct  analogy  in  magnetic  resonance,"  has  been  observed 
previously  in  the  two-dimensional  electron  gas  on  helium" 
is  also  explored  briefly. 

Let  us  begin  with  the  experimental  result.  Noncrystalline 
or  polycrystalline  KReO«  films  are  thermally  evaporated 
from  a  powder  onto  various  metal  or  transparent  substrates. 
In  Fig.  1  the  transmission  spectrum  of  a  =  1000  A  KRe04 
film  on  a  KCI  substrate  is  shown  for  TM  polarized  radiation 
incident  at  45°.  The  polarization  and  angular  dependencies 
of  these  absorption  line  strengths  show  clearly  that  the 
lower  frequency  (=908  cm"1)  feature  is  due  to  parallel 
modes  while  the  sharper  line  at  =951  cm"1  is  associated 
with  the  perpendicular  modes.  As  we  expect  the  film  to  be 
nominally  isotropic,  the  substantial  difference  in  the  width 
and  shape  of  the  perpendicular  and  parallel  mode  absorp¬ 
tions  is  surprising.  We  will  explain  this  difference  without 
invoking  any  anisotropy,  as  a  consequence  of  the  collective 
(i.e.,  dipole-coupled)  nature  of  the  perpendicular  mode  os¬ 
cillations  in  the  thin-film  system.  (An  excellent  introduc¬ 
tion  to  the  properties  of  the  vibrational  modes  of  thin  films 
and  other  smaller  systems  is  given  in  Ruppin  and  Engl- 
man.2) 

The  system  under  study  consists  of  a  thin  film  of  thick¬ 


ness  d,  sandwiched  between  a  substrate  (index  of  refraction, 
nj)  and  air  («i-  1).  As  we  are  interested  in  obtaining  only 
the  9-0  (optical)  response,  we  can  model  the  film  by  a 
complex  dielectric  function  «2(<*<).  (To  obtain  finite  q 
response  more  sophisticated  approaches,  such  as  the  one 
described  by  Persson  and  Ryberg,1  must  be  used.)  For  such 
a  three-layer  system  the  linearized  equation  for  the  change 
in  the  substrate  reflectance  due  to  the  presence  of  the  film 
is  given  in  Ref.  14.  From  the  appropriate  Fresnel  equation15 
one  can  derive  the  analogous  equation  for  the  change  (due 
to  the  film  presence)  in  the  transmittance  of  TM  polarized 
radiation 

—  2o>— CoICii  cos20 ImUjfw)  1 
T  c 

+  C4  sin20  Im[  — t{  '(o>)  ])  ,  (1) 


FIG.  1.  Absorption  lines  associated  with  the  mode  of  Re04 
in  a  KRe04  film  are  shown  for  S-45*  The  broad  resonance  near 
910  cm"'  is  due  to  oscillations  polarized  in  the  film  plane,  while  the 
sharper  and  weaker  line  near  950  cm"1  is  due  to  perpendicular  os¬ 
cillations. 
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where  to  and  9  are  the  frequency  and  angle  of  incidence  of 
the  incident  radiation,  and 


\  2  1 

l/Z 

Co- 

"1 

-  — ~  sinJ0 
i  "J 

+  «jCOS0 

(2) 

(3) 


Ci-nfni  . 


(4) 


Although  not  really  necessary  for  calculation  Eq.  (1)  is  heu- 
ristically  useful  since  the  term  containing  lm(«2(<u)l  is  asso¬ 
ciated  with  absorption  due  to  electric  dipole  oscillations  in 
the  plane  of  the  Him  (parallel  modes),  while  the 
Im(  —  '(tu)l  term  is  due  to  the  perpendicular  modes.  We 
will  now  examine  the  resonances  in  these  response  func¬ 
tions  for  several  simple  models. 

For  a  single  mode  Lorentz  oscillator16  (frequency  <o„ 
damping  y,  mode  strength  toj)  with  static  screening  e„ 

«2(<*»)  —  tm  +  to}/(to1  —  to1  —  iyto)  ,  (5) 


lm[«2(«<)]  has  a  resonance  at  the  frequency  to,  of  width  y 
corresponding  to  the  absorption  due  to  the  parallel  mode. 
On  the  other  hand,  the  longitudinal  response  function 


Im(-<2_l(w)] 


Im 


Op/ 


to1  +  (lip!  —  to1  +  i  yto 


(6) 


which  is  associated  with  the  perpendicular  oscillations,  has 
its  resonance  at  the  higher  frequency 


toi  —  to,  +  <u,7«.  . 


(7) 


With  respect  to  the  parallel  mode  the  perpendicular  mode 
resonance  is  (1)  shifted  upward  in  frequency  by  approxi¬ 
mately 


2<_w  T 


(8) 


(2)  weaker  in  intensity  by  **.,  and  (3)  equal  in  linewidth. 
(1)  and  (2)  are  consistent  with  the  experimental  results 
(Fig.  1);  however,  (3)  is  not  since  in  the  experiment  we 
find  that  the  parallel  mode  absorption  is  significantly 
broader  than  that  associated  with  the  perpendicular  mode. 

In  order  to  explain  this  linewidth  discrepancy  let  us  con¬ 
sider  the  introduction  of  inhomogeneous  broadening  to 
t}(to),  the  simplest  case  of  which  involves  a  two  oscillator 
dielectric  function: 


«»<«>-«-+  i  -2—^ — •  (9) 

In  this  case  limlt-fm)),  as  shown  in  Fig.  2,  has  a  doublet 
resonance  with  peaks  at  u>,,  and  tu,2  corresponding  to  ab¬ 
sorption  due  to  parallel  modes.  Surprisingly,  however, 
Im(  -  tf  '(to)  I  does  not  mimic  this  doublet  structure,  but 
has  only  a  single  sharp  resonance  at  a  frequency  between 
«ii  +  (<oJi  +wj:)/(m  and  tol}  +  (to,,  +topi)/tm.  This  tenden¬ 
cy  of  the  resonance  in  Im[  -  «2~  '(to) )  to  remain  narrow 
despite  the  inhomogeneous  broadening  in  «2(<u)  is  found  to 
prevail  whenever  the  perpendicular-parallel  mode  splitting 
( )  is  greater  than  the  separation  of  the  parallel  modes 

(w,,  —  to,}). 

It  is  now  clear  why  a  transmission  spectrum  calculated  us¬ 
ing  Eqs.  (1)  and  (9)  can  compare  favorably  to  the  data 
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FIG.  2.  Resonances  of  the  two  mode  Lorentz  oscillator  dielectric 
function  are  shown  as  a  function  of  frequency.  In  Im(-c-1)  the 
double  resonance  put  into  c  has  coalesced  into  a  single  line. 


shown  in  Fig.  1.  The  striking  point  is  not  that  we  can  fit  the 
data  with  Lorentz  oscillators,  but  that  the  (doublet)  struc¬ 
ture  introduced  to  fit  the  parallel  modes  does  not  produce  a 
corresponding  structure  in  the  perpendicular  absorption  line. 

Departing  from  the  data  of  Fig.  I,  we  can  explore  this 
electric  dipole-dipole-interaction  narrowing  by  considering 
an  inhomogeneously  broadened  line  with  a  Gaussian  profile, 
i.e., 


2  r- 

*2(10) -•«..  + -7=—  J  — - - - dto  , 

v  ir A  ^  "•  <u  ““  cu  —  /  yo) 


(10) 


Here  to,  is  the  center  frequency  of  the  Gaussian  packet  and 
A  and  y  are  the  inhomogeneous  and  homogeneous 
linewidths.  respectively.  (We  assume  A  >  y.) 

The  resonance  in  Im[e(r*i)],  which  is  associated  with 
parallel  mode  absorption,  always  exhibits  the  full  inhomo- 


FIG.  3.  Resonances  of  a  Gaussian  distribution  of  Lorentz  oscilla¬ 
tors  are  shown.  The  resonance  in  lm(<>  is  Gaussian  for  all  a.  In 
contrast,  the  resonance  in  Im(-e_l)  crosses  over  from  a  broad 
Gaussian  to  the  narrow  Lorentztan  in  the  neighborhood  of  a  —  l,  as 
illustrated. 
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gcneous  width  A,  as  shown  in  Fig.  3.  The  width  and  shape 
of  the  resonance  in  Im[  —  <f ’(<u ) ].  however,  depends  on 
the  dimensionless  ratio  of  coupling  strength  to  inhomo- 
geneily: 

a  — at,/ A  . 

For  the  strongly  coupled  system  (a  »  1)  the  resonance 
exhibits  the  homogeneous  width  (y)  and  shape,  while  in 
the  weak-coupling  regime  (a«  1)  the  line  is  Gaussian  of 
width  A.  An  abrupt  crossover  occurs  near  a  -  1.  The  line 
shape  and  width  of  perpendicular  modes  thus  depend  criti¬ 
cally  on  a. 

This  line  narrowing  may  actually  have  been  observed  for 
the  perpendicular  (subband)  excitations  of  electrons  bound 
to  the  surface  of  He.  In  this  system  inhomogeneous 
broadening  can  be  introduced  in  a  controlled  manner  via  a 
parallel  magnetic  Held.  It  has  been  observed  that  with  in¬ 
creasing  areal  electron  density  a  remarkable  narrowing  of 
the  resonant  absorption  line  occurs.13  Although  an  explana¬ 
tion  in  terms  of  changes  in  the  velocity  autocorrelation 
function  has  been  put  forward,  dipole-dipole-interaction- 


induced  narrowing  should  also  be  considered,  especially 
since  the  observed  line  narrowing  occurs  for  increasing  a 
near17  a  -  l. 

In  summary,  we  have  shown  that  substantial  differences 
observed  in  the  linewidth  and  shape  of  parallel  and  perpen¬ 
dicular  modes  of  KrReO.*  film  can  be  accounted  for  by  a 
dipole-dipole-interaction-induced  narrowing  of  the  perpen¬ 
dicular  mode.  In  this  collective  mode  individual  molecules 
oscillate  in  phase  at  a  single  frequency  despite  having  differ¬ 
ing  natural  frequencies  and  thus  the  inhomogeneity  inherent 
in  the  system  is  obscured.  This  interaction-induced  line 
narrowing  can  occur  in  a  wide  variety  of  systems  involving 
two-dimensional  arrav  of  interacting  dipole  oscillators,  in¬ 
cluding  electrons  on  .  ,  planar  arrays  of  metal  particles,  in¬ 
version  layers,  and  adsorbate  layers,  as  well  as  in  thin  films. 
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Using  the  infrared-absorption  technique  of  surface-electromagnetic-wave  spectroscopy  (SEWS)  we 
have  studied  the  effects  of  H2  and  D2  chemisorption  on  the  SEW  attenuation  coefficient  a  in  the 
10-lim  region  between  165  and  350  1C.  The  change  in  a  at  room  temperature  has  also  been  mea¬ 
sured  for  N2,  02,  and  CO  adsorption.  The  coverage  (0)  dependence  of  a  for  N2,  02,  and  CO  is  fair¬ 
ly  simple  [either  a  monotonic  increase  (N2)  or  a  peak  at  some  intermediate  coverage  (02,  CO)]  and 
points  to  changes  in  free-carrier  (FC)  surface  scattering  as  the  dominant  a -changing  mechanism. 

Much  richer  structure  in  a<0)  upon  H2  and  D2  adsorption  is  observed:  A  small  peak  followed  by  a 
i  sharp  dip  at  9=0.42  monolayer  (ML)  (6  =  2  monolayers  =  saturation),  a  broad  maximum  at 

0=  1.28  ML,  and  a  value  at  0=2.0  ML  higher  than  for  the  clean  surface.  For  0>O.42  ML  the 
variations  are  dominated  by  reconstruction-induced  changes  in  FC  scattering  from  the  surface.  The 
difference  in  a  between  9=  1.28  ML  and  saturation  indicates  that  the  6=  1.28  ML  disordered  sur¬ 
face  phase  and  the  6=2.0  ML  p{  lx  1)-D  phase  are  temperature  independent  between  165  and  300 
K.  Adsorbate-induced  quenching  of  the  W(100)  0.3-eV  intrinsic  surface  states  also  contributes  to 
the  changes  in  a;  the  associated  loss  of  oscillator  strength  is  most  apparent  in  a  strong  frequency 
dependence  of  a  relative  to  the  clean  surface  for  D2  adsorption.  In  addition,  it  is  possible  that  a  de¬ 
crease  in  relaxation  of  the  W(100)  surface  as  9  approaches  0.42  ML  for  H2  or  D2  has  some  effect  on 
the  SEW  attenuation  by  influencing  the  amount  of  FC  interband  scattering  at  the  metal-vacuum  in¬ 
terface. 


I.  INTRODUCTION 

In  1952  Holstein'  showed  with  a  simple  physical  argu¬ 
ment  that  the  infrared  (ir)  absorptivity  at  a  metal  surface 
is  larger  if  the  conduction  electrons  colliding  with  the  sur¬ 
face  are  diffusely  rather  than  specularly  reflected.  The 
next  year  Dingle  further  developed  the  theory  of  the 
anomalous  skin  effect2,3  from  the  formal  equations  of 
Reuter  and  Sondheimer4  and  applied  his  theory5  to  a  large 
group  of  optical  data  which  resulted  in  qualitative  agree¬ 
ment  as  long  as  the  electronic  scattering  from  the  surface 
was  assumed  diffuse.  Other  experimental  results  on  the 
absorption  of  electromagnetic  radiation6-"  also  indicate 
diffuse  surface  scattering.  However,  Bennett  et  a/.,*2-14 
using  supersmooth  substrates  for  ultrahigh  vacuum 
(UHV)  evaporated  Au,  Ag,  and  A1  films  measured  in  a 
dry  N2  atmosphere,  were  able  to  obtain  ir  reflection  re¬ 
sults  which  were  consistent  with  specular  scattering. 
Furthermore,  by  increasing  the  roughness  of  the  sub¬ 
strate14  they  were  able  to  increase  the  value  of  the  diffuse 
scattering  and  thereby  decrease  the  reflectivity.  Their  re¬ 
sults  suggest  that  a  sample  of  sufficient  quality  can  exhib¬ 
it  specular  scattering. 

Although  the  early  prejudice  was  that  only  totally  dif¬ 
fuse  scattering  was  possible,  in  spite  of  a  small  amount  of 
experimental  evidence  to  the  contrary,15  many  more  re¬ 
cent  experiments16-22  show  that  specular  scattering  is 
possible  for  atomically  smooth  surfaces  and  that  contam¬ 
ination  due  to  chemisorption,  which  is  necessarily  present 
in  all  but  careful  UHV  studies,  need  not  inhibit  all  of  the 
specular  electronic  reflection  from  the  surface.  The  prob- 
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lem  has  been  and  continues  to  be  to  identify  specularity  in 
a  quantitative  manner.  With  the  development  of  UHV 
surface  probes  has  come  the  realization  that  surface 
reconstruction  is  a  natural  phenomenon  which  can  be  trig¬ 
gered  in  a  variety  of  ways.23  In  this  paper  we  describe  our 
infrared  investigation  of  free-carrier  scattering  at  a  metal 
surface  as  a  function  of  adsorbate-  and  temperature- 
induced  surface  reconstruction. 

The  W(  100)-H  surface  has  been  studied  by  a  variety  of 
UHV  techniques  which  has  enabled  a  fairly  detailed 
structural  picture  to  develop.  In  spite  of  the  large  quanti¬ 
ty  of  data  about  this  system  no  study  has  directly  dealt 
with  the  interaction  of  the  tungsten  conduction  electrons 
with  the  surface  upon  H2  adsorption.  Although  optical 
absorption  studies  have  been  done  on  this  system  they 
have  either  concentrated  on  energies  above  0.6  eV  where 
interband  transitions  dominate  the  optical  response,24,25  or 
they  have  addressed  only  the  question  of  the  vibrational 
response  of  the  adsorbate-adsorbent  system.26,27  Here  we 
use  the  ir  absorption  technique  of  surface  electromagnetic 
wave  spectroscopy  (SEWS)  in  the  125  meV  (1000  cm-1) 
region  to  study  the  surface  scattering  phenomenon. 

We  find  that  for  H2  and  D2  adsorption  the  largest  con¬ 
tribution  to  the  change  in  absorptivity  in  this  energy  range 
can  be  attributed  to  variations  in  the  free-carrier  scatter¬ 
ing  from  the  surface  caused  by  adsorbate-induced  surface 
reconstruction.  In  addition  there  is  a  contribution  due  to 
the  quenching  of  intrinsic  surface  states  (ISS)  located  near 
0.3  eV  which  are  inherent  to  this  surface.  In  the  next  sec¬ 
tion  we  discuss  background  information  on  free-carrier 
contributions  to  the  SEW  attenuation  coefficient  with  an 
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emphasis  on  the  surface-assisted  ir  absorption  process.  In 
Sec.  Ill  the  experimental  details  are  presented.  The  re¬ 
sults  for  the  change  in  the  SEW  attenuation  coefficient 
are  described  in  Sec.  IV:  Room  temperature  (RT)  expo¬ 
sure  and  coverage  dependence  for  H2  and  D2  is  presented 
along  with  the  frequency  dependence  and  temperature 
dependence  of  D2-induced  variations.  In  addition,  room 
temperature  exposure  results  are  presented  for  adsorption 
of  Nj,  02,  and  CO.  In  Sec.  V  we  discuss  these  results  and 
show  that  changes  in  the  bulk  Drude  parameters  in  a 
small  region  (several  lattice  constants)  near  the  surface  are 
too  small  to  account  for  the  data.  We  then  discuss  the 
variations  in  attenuation  in  terms  of  the  0.3-eV  ISS  and 
the  surface-assisted  absorption  process.  The  concluding 
section  summarizes  our  results  on  electronic  surface 
scattering  from  the  W(100)  surface. 

II.  FREE-CARRIER  CONTRIBUTIONS  TO  THE  SEW 
ATTENUATION  COEFFICIENT 

A.  Drude  model 

The  ability  of  any  metal  such  as  tungsten  to  absorb  in¬ 
frared  radiation  at  room  temperature  depends  upon  the 
behavior  of  the  conduction  electrons  within  a  skin  depth  5 
( ss  100  A )  of  the  surface.  More  specifically,  in  the  Drude 
model,  the  conduction  electron  contribution  aM  to  the 
SEW  attenuation  coefficient  a  for  a  metal-vacuum  inter¬ 
face  is  related  to  the  relaxation  time  r  of  the  charge  car¬ 
riers  within  5  by 
12 


co 


—  UO«0}„), 
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where  co  is  the  frequency,  cop  the  plasma  frequency  (both 
in  rad/sec),  and  c  the  speed  of  light.28  In  the  infrared 
near  room  temperature  there  are  three  contributions  to 
1/t  which  must  be  considered:  the  electron-electron 
scattering  rate  1/t,.,,  the  electron-phonon  scattering  rate 
1/t*,  and  an  effective  scattering  rate  due  to  collisions  with 
the  surface  1/t,.  These  processes  are  often  treated  as  in¬ 
dependent  so  that 

—  +  —  +  — ; 

T  T,„  T*  T, 

hence,  we  first  consider  each  one  separately. 

1.  Electron-electron  scattering 

The  electron-electron  scattering  rate  has  been  theoreti¬ 
cally  considered  by  Gurzhi29,30  and  Lawrence.31  Both  au¬ 
thors  have  derived  expressions  for  the  scattering  rate  of 
the  form 
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Gurzhi’s  order-of-magnitude  expression  has  K~cop/ 
)2  which  for  tungsten  is  ss2Xl014  s-leV-2. 
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TABLE  I.  Room-temperature  ir  parameters  for  W. 


Parameter 

Value  at  '000  cm'1 

Pa c* 

5.59  pH  cm 

feu/ 

7.0  eV 

Vf/cC 

4.0X  10“ J 

Vp/(0 

63  A 

6“ 

140  A 

1/T,./ 

1.4X1013  s’1 

1/t/ 

5.3X  10n  s'1 

1/T,* 

( 1  —  p)1.6x  1011  s"1 

4T  h 

70 

'Reference  41. 

'Reference  34. 

bReference  32. 

fSee  text. 

Reference  52. 

•Reference  1. 

dReference  50. 

hReference  62. 

Experimentally,  the  co2  and  T2  signatures  allow  1/t,., 
to  be  determined  from  temperature-dependent  emissivity 
measurements.  For  the  transition  metals  Mo,  V,  Ta,33 
and  W  (Refs.  33  and  34)  1/t,.,  appears  to  be  substantially 
larger  than  either  of  the  theoretical  estimations  predict. 
In  particular,  for  tungsten  Wojcik  el  al ,34  have  deter¬ 
mined  a  value  of  AT  =  1.4x  10lft  s-leV-2.35  At  room 
temperature  and  1000  cm-1  this  yields  a  value  of 

l/r,„  =  1.4x  1013  --1  ; — ~ 


s  1  which  implies  a,.,  =0.15  cm' 


2.  Electron-phonon  scattering 

At  zero  frequency  near  room  temperature  electron- 
phonon  scattering  is  the  dominant  relaxation  mechanism 
for  the  free  carriers  and  provides  the  major  contribution 
to  the  dc  resistivity  pdc=4ir/[w^T<0)].36  For  0  an  ex¬ 
tra  decay  mechanism  exists  due  to  the  spontaneous  emis¬ 
sion  of  phonons  which  is  absent  in  the  zero-frequency  lim¬ 
it.37-39  This  spontaneous  emission  process  can  produce  a 
significant  addition  to  the  dc  relaxation  rate  1/t*(0)  for 
fuo>lkBT  and  feu>2ks0D  (0D  denotes  the  Debye 
temperature  of  metal)"?9 

In  the  limits  fko»kBT  and  fuo»kBQD  Holstein’s37-38 
and  Gurzhi’s39  results  indicate  at  high  temperatures 
(  r»0fl)  that  1/t*=  1/t*(0).  But  as  the  temperature  is 
lowered  the  spontaneous  emission  process  gains  increasing 
importance  in  the  ir  absorptivity  until  7=0,  at  which 
point  it  constitutes  the  only  lattice  vibrational  decay 
mechanism  of  the  electrons.  For  temperatures  which  are 
not  too  low  [T>Qd/ 4)  the  relationship  between  the  ir 
and  dc  relaxation  rates  is  well  represented  by  a  second- 
order  expansion  in  (&D/T)2: 


(1) 
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3240 

7 

For  tungsten  ©p  =310  K,30  so  a  12%  correction  to  the  dc 
rate  results  in  the  ir  at  room  temperature. 

Since  l/r*(0)  and  1/t,.,(0)  are  the  only  significant  RT 
contributions  to  p^,  the  electron-phonon  relaxation  rate 


Table  I  for  cop  (Ref.  32)  and  other  ir  parameters  for  W.]  can  be  calculated  from  p,*.,41  1/t,.,(0),  and  Eq.  (1) 


Lawrence's  expression,  which  works  fairly  well  for  the  no¬ 
ble  metals  Cu  and  Ag,  yields  values  of  K  <  10i5  s-1  eV-2. 


Thus  l/r*  =  5.3x  1013  s  1  which  implies  a*=0.55 


cm  1  =:  3.7a,.,  at  1000  cm 
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3.  Electron-surface  scattering 

In  the  surface-assisted  optical  absorption  process  it  is 
important  to  distinguish  between  scattering  events  which 
conserve  the  tangential  velocity  u,|  of  the  impinging  elec¬ 
trons,  which  is  defined  as  specular  scattering,  and  those 
events  which  do  not  since  conserving  events  contribute 
negligibly  to  the  overall  absorption.'  The  causes  of  non- 
specular  scattering  can  be  conveniently  divided  into  two 
categories. 

The  first  category,  which  was  considered  by  Holstein 
and  Dingle,  is  composed  of  causes  which  are  classified  as 
surface  disorder  and  leads  to  diffuse  scattering.  For  a  sin¬ 
gle  crystal  examples  of  such  inhomogeneities  are  missing 
or  extra  atoms  in  the  top  layer,  steps  in  the  surface  due  to 
slight  misalignment  along  a  crystallographic  axis,  or,  as  is 
relevant  to  our  discussion,  adsorption  of  molecules  onto 
the  surface. 

Nonspecular  scattering  also  arises  because  of  the  com¬ 
plex  nature  of  real  Fermi  surfaces  and  is  present  even  in 
single  crystals  with  perfect  surfaces.43  -46  We  refer  to  this 
process  as  surface  umklapp  scattering.  To  illustrate  this 
for  tungsten  a  slice  of  the  Fermi  surface,47  normal  to  the 
(001)  direction,  is  drawn  in  Fig.  1.  It  is  composed  of  elec¬ 
tron  “jacks”  and  hole  “octahedrals”  and  “ellipsoids.”  An 
electron  located  on  the  Fermi  surface  at  point  A  will 
strike  the  ( 100)  metal-vacuum  interface  at  an  incident  an¬ 
gle  of  45*.  If  the  sample  surface  has  the  same  periodicity 
as  the  lattice  then  the  electron  must  be  reflected  to  either 


FIG.  1.  A  slice  of  the  W  Fermi  surface  perpendicular  to  the 
(001)  direction.  An  electron  on  the  Fermi  surface  at  point  A 
(the  vector  indicates  the  direction  of  the  electron’s  velocity)  in¬ 
cident  on  the  (100)  metal-vacuum  interface  can  be  reflected  to 
either  point  B  or  C  with  no  change  in  tangential  crystal  momen¬ 
tum.  Scattering  into  B  is  specular  while  scattering  into  C  is  a 
surface  umklapp  process. 


Frequency  (cm  ’) 

FIG.  2.  a  /a+  versus  frequency.  The  dashed  lines  show  vari¬ 
ous  contributions  in  the  Drude  model  (a*,a,^)  with  Holstein’s 
approximation  for  diffuse  surface  scattering  (a,d).  The  solid 
lines  are  curves  calculated  from  Dingle’s  anomalous  skin  effect 
equations. 


point  B  or  C  since  the  tangential  crystal  momentum  fik  „ 
is  necessarily  conserved.  If  it  emerges  in  the  state  B  then 
the  scattering  is  specular  since  the  tangential  velocity  is 
also  conserved.  Interband  scattering  from  the  free-carrier 
(FC)  electron  band  to  the  FC  hole  band  into  the  state  C, 
however,  results  in  a  large  change  in  the  velocity  parallel 
to  the  crystal  surface,  which  contributes  to  ir  absorption. 
An  atomically  perfect  surface,  therefore,  need  not  appear 
perfectly  specular. 

The  overall  scattering  behavior  of  conduction  electrons 
from  the  surface  is  phenomenologically  described  by  the 
Fuch’s  specularity  parameter  p.48  If  all  the  electrons  are 
reflected  with  no  tangential  velocity  change  then  p=  1. 
p= 0  indicates  that  the  electrons  are  randomly  scattered 
into  all  possible  scattering  angles.  Other  types  of  behavior 
can  be  characterized  by  a  value  of  p  between  0  and  1 
which  indicates  the  effective  specularity  of  the  outgoing 
electrons.  For  example,  if  all  the  electrons  arc  reflected 
into  angles  near  the  specular  direction  then  p  is  close  to  1 
even  though  there  is  no  true  specular  scattering  occur¬ 
ring.49 

In  order  for  the  amount  of  specularity  to  have  a  major 
effect  on  the  absorptivity  of  an  optically  thick  sample  two 
conditions  must  apply.  The  first  is  that  the  mean-free 
path  l0  of  the  electrons  must  be  at  least  comparable  to  the 
classical  skin  depth  5  (Ref.  50)  of  the  metal.  If  / 0«6 
then  bulk  scattering  dominates  the  absorptivity.  The  oth¬ 
er  condition  is  that  the  energy  ftco  of  the  impinging  radia¬ 
tion  be  below  the  region  where  interband  transitions  dom¬ 
inate  the  optical  response  of  the  metal.  For  good  conduc- 
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tors  at  room  temperature  the  range  between  100  and  4000 
cm-1  generally  satisfies  these  two  conditions,91  below 
which  l0  «8  and  above  which  interband  transitions  have 
the  largest  effect. 

Within  this  range  of  frequencies  there  are  two  regions 
separated  by  8 =vF/eo  (Ref.  52)  ( vF  is  electronic  velocity 
at  Fermi  surface).  If  the  skin  depth  8  is  much  larger  than 
the  distance  traveled  by  an  electron  in  one  cycle  of  the 
field  then  the  dielectric  function  is  local  in  space  and  the 
surface  and  bulk  contributions  can  be  separated.1  In  this 
case  the  surface  contribution  to  the  free-carrier  response  is 
characterized  by  a  scattering  rate 


1  1 1  » 3  Vf 
--a 


For  tungsten  this  equals  (1— p)1.6xl013  s-1  which  im¬ 
plies  a,=(l— p)0.16  cm-1  at  1000  cm-1.  In  the  diffuse 
limit  a,irfsr0.3a^.  In  Fig.  2  a*.  ad+at+,  and 
a*+ae-*+as,d  normalized  by  at  are  plotted  versus  fre¬ 
quency. 


B.  Anomalous  skin  effect 

Since  for  tungsten  at  1000  cm-1,  8a=2 vF/co,  Holstein’s 
local  description  of  the  surface-assisted  absorption  process 
is  not  quite  valid.  Also,  because  SEW's  are  a  p- polarized 
(TM)  wave  the  effect  on  the  absorptivity  of  microscopic 
fields  located  at  the  surface93  must  be  considered  along 
with  the  nonlocality  of  the  dielectric  function.  Fortunate¬ 
ly,  as  Kleiwer  and  Fuchs  have  shown, 94,99  for  co.  and  r 
appropriate  for  W  near  v— 1000  cm-1  the  surface  im¬ 
pedance  z  for  non-normal  incidence  is  independent  of  the 
polarization  of  the  radiation  and  is  essentially  equivalent 
to  that  given  by  the  theory  of  the  anomalous  skin  effect  as 
developed  by  Reuter  and  Sondheimer  (RS).4 

In  the  frequency  region  of  interest  Dingle2,3  has  calcu¬ 
lated,  using  RS’s  equations,  an  expression  for  the  surface 
impedance  z  which  is  appropriate  for  W  at  room  tempera¬ 
ture.  Using  the  general  relationship  that  a 
=  —  2 ((o/c)rnx„  (Ref.  56)  where  rn  (xH )  is  the  real  (imagi¬ 
nary)  part  of  the  normalized  surface  impedance 
z,  =z/(4w/c)  we  have  plotted,  also  in  Fig.  2,  a /a+  for 
several  values  of  p  for  the  frequency  region  covered  in  our 
measurements.  In  the  b»vF/u>  limit  the  results  for  p= 0 
and  1  approach  the  Drude-Holstein  results  for 
(a*+a,.,  +ald)/a4  and  (ad+a,_,)/a+,  respectively. 


III.  EXPERIMENTAL 


tion  absorption  spectroscopy  when  using  a  strong  ir 
source  such  as  a  laser.97 

The  sample  is  mounted  on  a  Huntington  600-XYZTRC 
manipulator  in  the  UHV  system  and  is  thermally  connect¬ 
ed  to  a  liquid-N}  cooling  assembly  as  shown  in  Fig.  3. 
Two  feedthroughs  on  the  manipulator  flange  are  connect¬ 
ed  on  the  inside  via  a  stainless-steel  tube  (/)  through 
which  liquid-N2  flows.  At  the  lowest  point  of  the  tube  a 
copper  block  (h)  is  clamped  onto  it.  Attached  to  the 
copper  block  two  0.025-cm-thick  flexible,  annealed  copper 
strips  (e)  run  down  to  the  sample  mount  where  they  are 
sandwiched  in  place.  The  sandwich  is  necessary  to  electri¬ 
cally  isolate  the  sample  during  cleaning  or  thermal 
desorption  runs  [since  it  is  placed  at  high  voltage  via  tan¬ 
talum  leads  (/)]  from  the  liquid-N?  tube  which  is  at 
ground.  The  sandwich  consists  of  a  stainless-steel  block 
(j)  which  is  connected  to  and  held  by  the  manipulator 
shaft  (g),  the  copper  cooling  strips  (e),  a  1.9-cm- 
diamx0.2-cm-thick  sapphire  flat  id),  a  0.025-cm  Ta 
sheet  piece  (c)  indirectly  connected  to  the  sample,  and  a 
macor  piece  (/')  for  insulating  the  screws  ik)  which  are 
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A.  Sample  details 

The  tungsten  sample  (6.1  cmX0.6  cmx0.17  cm)  was 
oriented,  spark  cut,  and  then  polished  on  both  of  its  larg¬ 
est  sides.  Two  gratings,  used  for  coupling  free-space  radi¬ 
ation  to  the  SEW’s,  were  then  etched  into  one  of  the  pol¬ 
ished  surfaces  with  standard  photolithographic  techniques 
utilizing  type  J  photoresist  and  a  5:1  mixture  of  HF  and 
HNOj  are  the  etchant.  The  32-pm  spacing  gratings  are 
an  estimated  10  pm  deep  and  are  separated  by  5  cm.  The 
large  distance  between  the  gratings  enhances  the  already 
high  sensitivity  SEWS  affords  over  conventional  reflec- 


FIG.  3.  A  front  sketch  of  the  W  sample  mount.  The  sample 
(a)  is  thermally  connected  to  the  liquid-Nj  cooling  feedthrough 
tube  (/)  via  Ta  rods  (6),  Ta  sheets  Ic),  sapphire  flats  (d),  Cu 
sheet  ( e ),  and  a  Cu  block  ih).  Macor  pieces  ( i)  and  the  sapphire 
flats  electrically  isolate  the  sample  from  ground  during  e-beam 
heating  which  is  provided  by  the  Ta  filament  in)  and  the  high- 
voltage  Ta  leads  ( /).  Temperature  is  monitored  with  a  W-Re<3 
at.  %)— W-Re(25  at.  %)  thermocouple  ( m ). 
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used  to  squeeze  the  sandwich  together  (only  shown  dashed 
on  left-hand  side).  The  Ta  strips  are  spotwelded  to  the  Ta 
rods  (b)  which  are  electron-beam  welded  to  the  ends  of 
the  tungsten  sample  (a).  By  maintaining  a  constant  flow 
of  cold  N2  gas  (or  liquid)  through  the  stainless-steel  tube 
the  sample  cools  from  2500  K  to  room  temperature  in  8 
min  and  to  165  K  in  25  min. 

Sample  temperature  is  measured  with  a  W-Re<3 
at. %)— W-Re<25  at. %)  thermocouple  (tc)  (m)  which  is 
spotwelded  to  the  center  of  the  top  edge  of  the  sample. 
Above  298  K  standard  tables  were  used  to  determine  the 
temperature,  but  below  this  the  tc  was  calibrated  against  a 
Chromel-Alumel  tc.  This  calibration  differed  by  no  more 
than  ±2.0*C  from  the  curve  published  by  Sandstrom  and 
Withrow.5* 

Before  any  measurements  were  made  the  sample  was 
baked  in  IX  10"7  Torr  of  02  for  about  20  h  at  1520  K 
(with  periodic  flashes  to  2300  K)  via  electron  bombard¬ 
ment  heating  from  a  Ta  filament  (n)  directly  beneath  the 
sample.  After  this  treatment,  the  low-temperature 
(v/2x'/2)R45*  low-energy  electron  diffraction  (LEED) 
pattern  was  observed,  indicating  that  this  procedure  was 
sufficient  to  remove  any  C  from  the  surface.59  Before 
each  experimental  run  the  sample  was  again  flashed  to 
2300  K  to  remove  any  contaminants,  and  the  low- 
temperature  LEED  pattern  was  periodically  observed  as  a 
check  against  subsequent  C  buildup. 

B.  Optical  layout 

The  32-/zm  spacing  of  the  gratings  was  chosen  so  that 
the  SEW’s  and  free-space  radiation  couple  near  an  angle 
of  45*  at  1000  cm"1  as  shown  in  Fig.  4.  S-polarized  (TE) 
radiation  from  a  COz  laser  (either  15C02  or  l2C02)  is  sent 
through  a  chopper  ( C)  and  then  a  half-wave  plate  (WP) 
which  is  set  to  produce  p-polarized  (TM)  radiation.  (Be¬ 
cause  SEW’s  are  a  TM  wave  they  only  couple  to  TM 
free-space  radiation.)  Most  of  the  radiation  passes 
through  a  beam  splitter  (BS)  and  enters  the  vacuum 
chamber  through  a  BaF2  window  (W\).  A  40-cm  focal 
length  (FL)  lens  (L2)  is  used  to  focus  the  radiation  onto 
the  input  grating  where  the  SEW’s  are  generated.  The 
surface  waves  then  travel  the  5  cm  across  the  W(100)  sur¬ 
face  to  the  output  grating  where  they  are  converted  back 
to  free-space  radiation  which  travels  through  another 
BaF2  window  (  W2)  and  impinges  on  a  liquid-N2-cooled 
HgCdTe  photoconductive  detector  ( D2). 

Although  the  free-space  radiation  generated  at  the  out¬ 
put  grating  is  well  collimated  its  angle  of  departure 
changes  as  the  frequency  of  radiation  is  changed.  Z.3  is  a 
10-cm  FL  lens  which  images  the  output  grating  onto  D2 
and  hence  minimizes  the  need  to  realign  the  detector 
whenever  the  laser  frequency  is  changed.  The  mirror 
( M3)  mounted  on  a  translation  stage  also  minimizes  prob¬ 
lems  in  maximizing  the  detected  signal  whenever  a  new 
laser  line  is  selected  since  translating  M3  perpendicular  to 
the  axis  of  L2  changes  the  angle  of  incidence  without 
altering  the  position  of  the  incoming  radiation  on  the 
sample. 

The  radiation  which  does  not  pass  through  the  beam 
splitter  is  focused  onto  a  pyroelectric  detector  (Dl)  by  LI. 
The  output  signals  of  D\  and  D2  are  sent  to  separate 


FIG.  4.  Top  view  of  the  optical  setup.  A  chopped  C02  laser 
beam  is  sent  into  the  sample  chamber  through  a  BaF2  window 
( 1VI)  and  strikes  the  input  grating  of  the  sample.  A  HgCdTe 
detector  ( D2)  is  used  to  measure  the  transmitted  intensity  of 
free-space  radiation  generated  at  the  output  grating  by  SEW's 
which  have  traveled  across  5  cm  of  the  front  (100)  surface. 
LEED  can  be  performed  on  the  rear  surface  as  ir  measurements 
are  made  on  the  front.  The  ratio  of  the  phase  detected  signals 
from  D2  and  the  reference  pyroelectric  detector  D\  are  recorded 
on  a  chart  recorder. 


lock-ins  tuned  to  the  chopper  frequency  and  their  outputs 
are  ratioed  and  recorded. 

C.  Experimental  procedure 

A  typical  experimental  run  consists  of  the  following: 
The  sample  assembly  is  First  cooled  down  as  far  as  possi¬ 
ble  to  insure  as  rapid  cooling  as  possible  from  high  tem¬ 
peratures.  The  sample  is  then  flashed  to  2300  K  and  al¬ 
lowed  to  cool  to  the  desired  adsorption  temperature  at 
which  time  the  Ta  filament  is  turned  on  and  its  current 
adjusted  so  that  it  radiatively  stabilizes  the  sample  tem¬ 
perature.  Gas  is  then  leaked  into  the  chamber  (typically 
at  5  X  10-9  Torr  initially)  and  the  SEW  signal  is  moni¬ 
tored  as  a  function  of  exposure.  The  change  in  the  SEW 
attenuation  coefficient  a  between  two  coverages  0,  and 
Of  is  calculated  from  Aa  =  —  [1/(5  cm)]ln(  1 +  A///0) 
where  A /  is  the  change  in  the  detected  intensity  J0  be¬ 
tween  0,  and  0/. 

IV.  RESULTS 
A.  H2  and  D2  adsorption 

The  features  in  A  a  versus  exposure  are  highly  non¬ 
monotonic  and  are  qualitatively  the  same  at  all  frequen¬ 
cies  studied  from  886  to  1088  cm"1.  Figure  5  shows  two 
typical  room-temperature  runs,  one  for  H2  and  one  for 
D2.  Initial  absorption  results  in  a  small,  but  significant, 
rise  in  a  followed  by  a  sharp  minimum  at  an  uncorrected 
gauge  exposure  of  0.100±  0.004  (0.140  ±0.006)  langmuir 
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Exposure  (10~6  Torr  sec) 


FIG.  5.  Change  in  a  versus  exposure  for  H2  and  D2  at  a  laser 
frequency  of  1032  cm-1  and  a  temperature  of  298  K.  The  expo¬ 
sure  is  a  corrected  value  obtained  by  multiplying  the  ion  gauge 
pressure  by  2.3. 


(L)  for  H2  (D2)  (1  Ls  10_®  Torrsec).  This  is  followed  by 
a  rise  to  a  broad  maximum  at  0.45  ±0.04  (0.59  ±0.05)  L. 
(These  exposure  values  are  averages  of  all  room- 
temperature  runs.)  As  the  surface  becomes  saturated  a 
decreases  to  a  final  value  above  the  zero  coverage  value 
with  aHj  >  ctDj.  The  stronger  absorption  for  H2  is  mainly 
due  to  the  WrH  symmetric  stretch  (v,)  which  has  been 
previously  reported.27  Note  that  the  ratio  of  exposures  at 
minimum  a  between  D2  and  H2  runs  is  1.40±0.09,  in 
good  agreement  with  previous  results  that  the  sticking 
coefficient  is  the  same  for  both  isotopes.®0 

In  order  to  convert  exposure  to  coverage  we  use  thermal 
desorption  measurements  which  show  that  the  coverage  at 
the  minimum  is  0=O.42±O.O4  ML.  Using  this  as  a  cali¬ 
bration  point  and  the  fact  that  the  sticking  coefficient  is 
well  represented  by  s=C(  1—0/2)  (Refs.  60  and  61) 
where  C  is  a  constant,  the  coverage  for  a  given  exposure  is 
determined.  A  a  vs  D2  and  H2  coverage  is  shown  in  Fig.  6 
for  the  same  data  shown  in  Fig.  5.  The  broad  maximum 
now  occurs  at  a  coverage  of  0=  1.28  ±0.08  ML. 

Figure  7  shows  the  experimental  results  for  Aa/cr*  for 
D2  adsorption  between  various  features  (zero  coverage, 
0=0.42  ML  dip,  0=  1.28  ML  peak  and  saturation)  as  a 
function  of  frequency.  Data  for  D2  instead  of  H2  is 
presented  in  order  to  eliminate  complications  caused  by 
the  Vj  mode  absorption  which  is  centered  at  1070  cm-1 
for  the  p  ( l  X 1  )-H  state.  Note  that  the  data  labeled  B  and 
C  actually  show  a  decrease  in  |  Aa/a*  |  as  a  increases 
while  the  data  labeled  A,  D,  and  E  show  ai\  increase. 

The  same  basic  features  in  A  a  versus  coverage  for  D2 


Coverage  B 

FIG.  6.  Change  in  a  versus  coverage  for  H2  and  D2  at  a  laser 
frequency  of  1032  cm-1  and  a  temperature  of  298  K.  Coverage 
is  determined  from  partial  coverage  thermal  desorption  mea¬ 
surements. 


FIG.  7.  Change  in  a  versus  frequency  for  D2  adsorption  be¬ 
tween  various  initial  (02)  and  final  (0/)  coverages.  Curve  A, 
0,  =0.42  ML,  0/=  1.28  ML;  curve  B,  Q,~ 0  ML,  0/=1.28 
ML;  curve  C,  0,  =0  ML,  0/  =  2.O  ML;  curve  D,  0,  =  1.28  ML, 
0/=2.O  ML;  curve  E,  0,  =0  ML,  ©,  =  0.42  ML.  The  solid 
lines  are  fits  to  the  theory  of  the  anomalous  skin  effect  [Eq.  (4)], 
while  the  dot -dashed  lines  are  linear  guides  to  the  eye. 
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FIG.  8.  Change  in  a  between  various  coverages  as  a  function 
of  temperature.  Curve  A,  0(=O.42  ML,  0/=1.28  ML;  curve 
B,  0,  =0  ML,  0/=  1.28  ML;  curve  C,  0,  =0  ML,  0/  =  2.0  ML; 
curve  D,  0f  =  1.28  ML,  0/=2.O  ML;  curve  E,  0(=  0  ML, 
0/=O.42  ML. 


Exposure  (10-6  Torr  sec) 


FIG.  9.  Aa  versus  exposure  for  Nj,  02,  and  CO  at  room  tem¬ 
perature.  The  exposure  is  uncorrected  in  all  cases. 


and  H2  are  present  from  16S  to  330  K.  The  temperature 
dependence  of  these  changes  for  D2  is  shown  in  Fig.  8. 
We  have  drawn  curves  through  the  data  labeled  A  and  E 
as  guides  to  the  eye.  The  curves  through  B  and  C  are 
linear  fits  below  300  K  while  a  constant  fit  is  shown  for 
D.  Above  RT  there  is  no  data  for  6  =  2.0  ML  since  D2 
partially  desorbs  above  300  K. 


B.  Other  gases 

The  rather  complex  nature  of  the  coverage  dependence 
of  A a  for  H2  and  D2  prompted  a  look  at  the  adsorption  of 
other  gases  to  see  what  similarities  and  differences  might 
exist.  Figure  9  shows  RT  results  for  adsorption  of  N2,  02, 
and  CO  onto  the  W(100)  surface  as  a  function  of  expo¬ 
sure.  In  each  case  the  coverage  dependence  of  the  absorp¬ 
tion  is  much  simpler  than  for  H2  or  D2. 


V.  DISCUSSION 


A.  Changes  in  modified  Drude  parameters  at  the  surface 

Although  tungsten  is  a  transition  metal,  bulk  ir  date  are 
well  described  by  an  appropriate  plasma  frequency  a)p,  re¬ 
laxation  time  r,  and  high-frequency  interband  contribu¬ 
tion  e0  (Ref.  62)  in  a  modified  Drude  model.  One  possible 
consequence  of  chemisorption  is  the  alteration  of  these  pa¬ 
rameters  in  a  small  region  (1  or  2  atomic  layers)  near  the 
surface.  Since  the  three-layer  model  of  McIntyre  and 
Aspnes61  works  well  in  interpreting  the  results  of  the  W2- 
H  vj  mode  absorption,  we  use  it  here  to  show  that  changes 
in  the  SEW  attenuation  coefficient  a  due  to  such  effects 
are  too  small  to  account  for  the  data. 

For  a  sufficiently  thin  ( d  /25  « 1 )  layer  of  thickness  d, 
the  attenuation  due  to  the  layer  can  be  written  as 


ai=&ir*v2d  Im 
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where  em  (e;)  is  the  dielectric  function  of  the  metal  (layer) 
and  v  is  the  frequency  (cm-1).64  If  we  assume  a  modified 
Drude  response 
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for  the  metal  (k  =  m)  and  the  layer  (k  =/),  then  the  above 
effects  correspond  to  altering  a)p  t  rh  and  e0  !. 

The  change  in  a  due  to  modifying  the  plasma  frequen¬ 
cy  in  a  thin  layer  calculated  from  Eq.  (2)  and  Eq.  (3)  is 
given  by 


Aa  =  —  a. 


co 


p,m 


d_ 
25  ' 


Assuming  that  one  conduction  electron  is  tied  up  in  each 
W— H  bond  this  yields  Aa=1.5xI0_:9  cm-1  at  1032 
cm-1.  While  this  can  account  for  the  size  of  A  a  at  sa¬ 
turation  it  fails  at  lower  coverages.  Nor  does  it  seem  like¬ 
ly  that  the  many  variations  in  the  sign  of  Aa  as  0  in- 
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creases  can  be  explained  since  all  of  the  H  is  located  at  the 
bridge  site61  implying  that  the  adsorbate-substrate  bond* 
ing  is  basically  the  same  at  all  coverages.  . 

If  instead  fA/  is  changed  and  <a£/  is  held  at  a  fixed 
value  on  the  order  of  ai^m,  then  a  is  altered  even  less.  In 
this  case 

A  _®i_A  ±  ' 

Aa_3a<*0J1  „  Al  0,1 26 

»7x  10_4Ae^/-^-  cm-1  . 

20 

Since  d/(28)=  10-2  and  the  maximum  in  Af0i/  =  102 
(Ref.  66)  we  have  result  which  is  2  orders  of  magnitude 
too  small. 

The  largest  excursion  in  a  as  the  electron-phonon  or 
electron-electron  interaction  is  altered  in  a  thin  layer  near 
the  surface  occurs  as  a)Tt  increases  from  —  7  X 10-2  to 
—  10  and  results  in  Aa=3x  10-2  cm-1  for  d= 3  A  (=2 
atomic  layers).  While  this  is  nearly  large  enough  to  be  a 
possible  explanation,  more  reasonable  variations  in  r;,  on 
the  order  of  Ar//r=0.1  cannot  account  for  the  size  of  the 
changes  observed  here.  Hence  it  appears  that  changes  in 
bulk  parameters  near  the  surface  are  insufficient  to  ac¬ 
count  for  the  SEWS  observations. 


fuse  surface  scattering  as  the  dominant  mechanism  for 
0  >  0.42  ML. 

This  interpretation  is  corroborated  by  other  experimen¬ 
tal  results  on  this  system.  Figure  10  compares  the  SEW 
results  with  the  mainly  LEED-derived  phase  dia- 
gram23,60, 74,75  for  the  W(100)-H  surface.  (Henceforth  we 
use  H  genetically  unless  specifically  referring  to  O.)  At 
6=0.42  ML,  where  a  is  a  minimum,  the  surface  is  recon¬ 
structed  in  what  is  known  as  the  incommensu¬ 
rate— c  (2x2  )-H  phase  which  is  characterized  by  -f  order 
spots  split  into  a  quartet  along  the  ±(  10)  and  ±(01)  direc¬ 
tions.  As  more  gas  is  added  the  spots  of  the  y  order 
quartet  streak  along  directions  perpendicular  to  their 
splitting  (the  streaking  is  apparent  at  0=0.48  ML  for  our 
sample)  and  then  eventually  disappear,  indicating  a  loss  of 
order  in  the  one  dimension  followed  by  a  loss  of  order  in 
both  dimensions.  These  LEED  changes,  signaling  the  on¬ 
set  of  disorder,  are  accompanied  by  a  sharp  rise  in  the 
SEW  attenuation  until  the  maximum  at  6=1.28  ML. 
The  subsequent  fall  in  a  as  the  surface  approaches  the 
well  ordered  p  ( 1 X 1  )-H  state  at  saturation  occurs  con¬ 
currently  with  a  decrease  in  the  diffuse  LEED  back¬ 
ground.  The  high  degree  of  order  at  saturation  is  also  re¬ 
flected  in  high-energy  ion  scattering76  results  which  show 


B.  Surface-assisted  free  carrier  absorption:  6  >0.42  ML 

The  coverage  region  for  6  >0.42  ML  shows  best  the 
free-carrier  absorption  effect  for  H2  and  D2  adsorption. 
This  is  due  to  the  fact  that  for  6  <0.42  ML  the  0.3-eV  in¬ 
trinsic  surface  states  are  not  yet  quenched67,6*  and  still 
contribute  to  optical  absorption.24  For  this  reason  we  first 
concentrate  on  6  >  0.42  ML.  - 

Between  6=0.42  and  2.0  ML  the  coverage  dependence 
of  the  SEW  attenuation  coefficient  (see  Fig.  6)  is  similar 
to  results  for  resistivity  changes  upon  molecular  adsorp¬ 
tion  onto  thin  metallic  films.  71  In  a  number  of  these 
cases  Wissmann70  has  demonstrated  thht  the  Fuchs  theory 
of  diffuse  surface  scattering4*  applies.  The  general 
features  of  the  changes  in  resistivity  upon  chemisorption 
observed  in  these  studies  are  as  follows:  Initial  adsorption 
increases  the  diffuse  scattering  due  to  the  randomness  of 
adatom  placement  and  leads  to  an  initial  increase  in  the 
resistivity.  For  some  systems  the  diffuseness  is  a  max¬ 
imum  at  saturation,  but  for  others  there  is  a  peak  at  some 
intermediate  coverages).  Such  behavior  is  not  unexpected 
since  a  saturated  layer  is  often  more  uniform  than  one  at 
intermediate  coverages.  A  useful  analogy  is  the 
Nordheim  peak72  in  resistivity  of  a  binary  alloy  AXB\_X 
as  x  is  increased  from  zero  to  unity.  In  fact,  a  recent  cal¬ 
culation  including  the  effects  of  multiple  scattering73  has 
predicted  the  existence  of  such  a  maximum  for  an  atomi¬ 
cally  smooth  surface  contaminated  by  adsorbed  atoms. 
We  also  note  a  general  experimental  trend  that  the  sa¬ 
turated  coverage  resistivity  is  consistently  higher  than  the 
clean  surface  resistivity.  The  striking  similarity  between 
the  W(100)-H  (and  -D)  attenuation  coefficient  coverage 
dependence  and  the  resistivity  coverage  dependence  of 
thin  gas  covered  films  strongly  points  to  changes  in  dif- 


Coverage  9 

FIG.  10.  The  LEED  determined  phase  diagram  and  Aa 
versus  coverage  for  W(100)-D.  The  sharp  dip  occurs  for  the 
split— c(2x2)  phase  while  the  broad  maximum  occurs  for  the 
disordered  state.  The  boundaries  on  the  incommensu¬ 
rate— c(2x2)  phase  were  determined  from  our  LEED  and 
thermal  desorption  measurements. 
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that  the  top  layer  is  unreconstructed  and  in  ir  vibra'-onal 
spectra  which  show  that  the  Vj  and  2v:-/3,  absorption 
features  are  their  sharpest  at  0  =  2.0  ML.77 

The  same  basic  coverage  dependence  also  appears  in  the 
SEW  data  for  the  absorption  of  CO  and  02  as  shown  in 
Fig.  9:  an  initial  rise  in  a  followed  by  a  peak  at  some  in¬ 
termediate  coverage  and  a  decrease  towards  saturation. 
Although  the  0.3-eV  ISS  are  quenched  by  adsorption  of 
these  gases,  the  quenching  occurs  more  gradually;®7,68'78 
hence  the  loss  of  the  associated  oscillator  strength  can  be 
approximated  by  a  slowly  decreasing  background.  This 
appears  to  result  in  a  smaller  value  of  a  at  saturation 
compared  with  zero  coverage.  Nitrogen  adsorption,  on 
the  other  hand,  shows  no  decrease  in  a  after  the  initial 
rise,  suggesting  that  disorder  is  a  maximum  at  saturation 
for  this  system.  These  observations  for  the  other  gases 
further  support  the  thesis  that  changes  in  diffuse  surface 
scattering  are  being  observed  on  the  W(100)  surface. 

It  should  also  be  noticed  that  the  variations  in  a  for  H2 
and  D2  are  at  least  twice  as  large  as  for  the  other  adsor¬ 
bates.  This  seems  strange  since  it  is  expected  that  the 
larger  atoms  should  have  a  higher  cross  section  for  elec¬ 
tronic  scattering  and  suggests  that  it  is  the  direct  interac¬ 
tion  of  the  free  carriers  with  the  reconstructed  top  W 
layer  by  H  or  D  which  results  in  the  changes  in  a  for 
these  two  systems. 

In  order  to  determine  the  amount  the  specularity 
changes  we  use  the  theory  of  the  anomalous  skin  effect  as 
developed  by  Dingle.2-3  For  the  limited  frequency  range 
of  our  study  Dingle's  results  can  be  summarized  by  the 
following  equation: 

Aa/a,=  — [0.253  +  6.3xl0~5(v— 1000)]Ap  ,  (4) 

where  v  is  the  frequency  expressed  in  cm-1.  Using  this 
equation  we  have  determined  the  change  in  p  between 
0j=O.42  ML  and  ©y=l.28  ML  to  be  Ap  =  —  0.39  and 
between  0;  =  1.28  ML  and  0y  =  2.0  ML  to  be 
Ap  =  +0. 11.  The  fits  are  shown  in  Fig.  7  as  the  solid 
lines  through  the  data  labeled  A  and  D.  Although  the 
frequency  dependence  of  the  data  is  slfghtly  stronger  than 
the  theory,  both  show  an  increase  in  |  A  a  /at  |  as  to  in¬ 
creases. 

The  temperature  dependence  of  the  variations  in  a  in 
this  coverage  range  are  interesting.  The  data  labeled  D  in 
Fig.  8  shows  that  A  a  between  0,  =  1.28  ML  and  0y  =  2.0 
ML  is  temperature  independent  from  room  temperature 
down  to  165  K.  This  indicates  that  both  of  these  phases 
are  independent  of  temperature  over  this  range,  since  it 
seems  unlikely  that  the  disordered  and  pllxl)-D  states 
would  change  a  in  the  same  manner  with  temperature  if 
there  did  exist  a  change  in  either  one. 

Using  the  value  of  a  at  0y  =  1.28  ML  as  a  reference  the 
curve  labeled  A  shows  that  as  the  temperature  is  lowered 
the  relative  attenuation  at  0,  =0.42  ML  decreases  until 
approximately  200  K  where  there  begins  a  dramatic  in¬ 
crease  as  the  temperature  is  lowered  further.  The  slow  de¬ 
crease  suggests  that  the  specularity  for  this  phase  in¬ 
creases  slightly  but  that  below  200  K  it  sharply  drops  off. 
Such  a  sharp  drop  in  specularity  may  be  associated  with 
immobility  of  the  adsorbate,  as  suggested  by  Barker 
et  al .,74  which  hinders  complete  formation  of  the  surface 


phase.  This  immobility  seems  also  apparent  in  other  ir 
measurements  where  low-coverage,  low-temperature 
modes  are  observed  to  be  much  broader  than  at  room  tem¬ 
perature. 79 

C.  Surface  umklapp  scattering  and  the  intrinsic 
surface  states:  6  <0.42  ML 

For  ©<0.42  ML  it  is  apparent  that  something  else  be¬ 
sides  surface  assisted  absorption  contributes  to  the  ir 
response.  This  is  best  seen  in  the  frequency  dependence  of 
the  changes  in  a  for  0,  =O  ML  shown  in  curves  B,  C, 
and  E  of  Fig.  7.  For  0y=  1.28  and  2.0  ML  (curves  B  and 
C,  respectively)  |  \a/at\  shows  a  decrease  as  frequency 
increases  while  |Aa/cr*|  for  ©y  =  0.42  ML  (curve  E ) 
shows  a  much  stronger  dependence  on  frequency  than 
curves  A  and  D  which  show  good  agreement  with  the 
anomalous  skin  effect. 

The  resultant  strong  negative  slope  of  A a/at  versus 
frequency  from  these  three  0,  =0  ML  curves  can  be  un¬ 
derstood  in  terms  of  the  D2  induced  quenching  of  the 
0.3-eV  intrinsic  surface  states.  As  pointed  out  by  Ander¬ 
son  et  al.,2*  the  good  correlation  of  the  surface  reflec¬ 
tance  spectroscopy  (SRS)  transition  energies  with  the  ul¬ 
traviolet  photoemission  determined  energy  levels  strongly 
suggests  that  the  optical  transitions  from  surface  electron¬ 
ic  states  involve  final-state  energies  located  just  above  the 
Fermi  level.  Although  the  SRS  measurements  only 
probed  transitions  down  to  ==0.6  eV  the  data  show  a  con¬ 
tribution  from  these  intrinsic  states  up  to  =0.8  eV,  an  in¬ 
dication  that  these  states  are  rather  broad.  It  is  therefore 
reasonable  to  expect  a  contribution  from  the  same  states 
at  1000  cm-1  (0.125  eV).  The  strong  frequency  depen¬ 
dence  is  thus  due  to  the  fact  that  our  measurements  probe 
the  low-frequency  tail  of  these  states  where  the  absorption 
strength  is  rapidly  changing. 

Exactly  how  much  these  states  contribute  to  the 
response  is  an  open  question  and  cannot  be  determined 
from  this  study.  The  simplest  assumption  is  that  all  of 
the  decrease  in  attenuation,  after  the  small  initial  rise,  is 
due  to  the  quenching  of  optical  transitions  associated  with 
the  removal  of  the  surface  states.  Another,  perhaps  more 
interesting,  possibility  is  that  some  part  of  the  large  nega¬ 
tive  excursion  in  a  is  caused  by  the  surface-enhanced 
free-carrier  absorption  effect.  Although  surface  disorder 
is  an  unlikely  explanation  there  is  the  possibility  that  sur¬ 
face  umklapp  processes  play  a  role. 

Surface  disorder  is  largely  ruled  out  by  LEED  observa¬ 
tions  in  this  coverage  region.23  ®0-  4  5  80  The  clean  surface 
shows  a  mostly  p  ( 1  X  I )  order  (unreconstructed)  with  a 
hint  of  y  order  spots.  At  a  very  small  value  of  0  the  { 
order  spots  sharpen  considerably  as  the  surface  enters  the 
c(2x2)-H  phase.  With  further  coverage  the  y  order 
spots  split  into  a  quartet  as  described  above  as  the 
incommensurate— c(2x2)-H  state  is  formed,  with  the 
magnitude  of  the  splitting  linearly  proportional  to  the 
coverage  beyond  which  the  splitting  begins  to  cccur 
(0  =  0.26  ML  for  our  sample).  Al'hough  there  is  some 
debate  about  the  state  of  the  room-temperature  clean  sur¬ 
face  (ordered  versus  disordered),  both  of  the  hydrogen  sta¬ 
bilized  phases  definitely  show  long-range  order. 
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At  first  glance  one  might  assume  that  the  appearance  of 
y  order  LEED  spots  would  necessarily  be  associated  with 
the  opening  up  of  more  free-carrier  nonspecular  scattering 
channels.  A  look  at  Fig.  11,  in  which  a  slice  of  the  Fermi 
surface47  is  drawn  perpendicular  to  the  ( 100)  direction,  in¬ 
dicates  that  this  is  not  a  large  effect  for  W.  A  change  in 
A|i  corresponding  to  a  (y  y)  LEED  beam  is  drawn  and 
shows  for  the  majority  of  free  carriers  that  such  a  change 
in  parallel  momentum  is  impossible  since  this  momentum 
transfer  connects  relatively  few  states  on  the  Fermi  sur¬ 
face.  Therefore,  a  c(2x2)  structure  should  not  signifi¬ 
cantly  increase  the  diffuse  surface  scattering. 

What  may  be  occurring,  though,  is  a  change  in  the 
amount  of  surface  umklapp  scattering  caused  by 
adsorbate-induced  changes  in  the  relaxation  of  the  W(100) 
surface.  Since  changes  in  reconstruction  (periodicity 
parallel  to  this  surface)  can  change  the  outgoing  electrons’ 
Ic,,  it  seems  reasonable  that  a  change  in  the  first  inter¬ 
atomic  layer  spacing  (relaxation)  could  effect  the  outgoing 
fcj  and  hence  change  the  relative  amount  of  electrons 
which  are  specularly  reflected.  Such  a  change  in  the  re¬ 
laxation  of  the  first  layer  has  been  reported  by  King  and 
Thomas60  in  their  LEED  study  of  this  system.  As  0  in¬ 
creases  to  =0.40  ML  their  Fig.  16  shows  a  sharp  increase 
in  the  first  W  interlayer  spacing  of  —2%,  with  much 
smaller  subsequent  variations  as  the  surface  becomes  sa¬ 
turated. 

The  temperature  dependence  of  Aa  for  the  clean  sur¬ 
face  and  a  calculation  by  Fu  et  a/.*'  lend  mild  support  to 
this  idea.  Their  calculation  suggests  that  as  the  low- 
temperature  ( Vl  X  Vl  )R  45*  clean  surface  phase  [which  is 
quite  similar  to  the  c(2x2)-H  phase]  is  formed  from  the 
p(lXl)  RT  state  the  top  layer  spacing  also  increases  (by 
~4%)  to  an  unrelaxed  position.  A  look  at  curves  B  and 
C  of  Fig.  8  shows  that  this  phase  transition  is  also  accom¬ 
panied  by  a  decrease  in  the  surface  contribution  to  a,  fur- 


FIG.  It  A  slice  of  the  W  Fermi  surface  perpendicular  lo  the 
UOOt  direction,  illustrating  that  a  change  in  k  corresponding  to 
airy)  LEED  beam  is  forbidden  for  the  majority  of  free  car- 


ther  indicating  that  surface  relaxation  may  be  related  to 
specularity  through  the  surface  umklapp  mechanism. 

Although  the  results  of  Fu  et  al.  are  consistent  with  a 
surface  umklapp  scattering  explanation  of  the  linear  de¬ 
crease  in  the  clean  surface  contribution  to  a,  another 
mechanism,  related  to  the  surface  Debye  temperature  of 
W,  may  be  responsible.  The  dynamical  LEED  data  of 
Heilmann  et  al.*2  indicate  that  the  high-temperature 
p  ( 1 X  1 )  clean  surface  has  an  effective  surface  Debye  tem¬ 
perature  0d.s  =  21O±4O  K  while  the  low-temperature 
( v^2x  V2)R  45*  phase  has  an  effective  ©os=400±100 
K.  Since  in  the  temperature  range  of  interest  the 
electron-phonon  relaxation  rate  l/r^oc  1/0q,36  a  change 
in  ©Ds  should  effect  the  electron-phonon  scattering  in  a 
thin  layer  near  the  surface.  Using  Eqs.  (2)  and  (3)  to  esti¬ 
mate  the  change  in  a  caused  by  the  LEED  determined  in¬ 
crease  in  0DJ  yields  Aa  =  —  (8.0±3.0)X  10-J  cm-1  at 
1049  cm-1  for  a  thickness  equal  to  one  interlayer  spacing 
(1.6  A).  Experimentally,  the  change  in  a  for  the  clean 
surface  in  going  from  300  to  165  K  is  1.0xl0~2  cm-2. 
The  fact  that  the  experimental  and  calculated  values  of 
A  a  are  not  significantly  different  lends  support  to  the  idea 
that  the  change  in  0D  S  may  be  a  possible  explanation  for 
the  clean  surface  temperature  dependence  of  Aa. 

Lastly,  we  note  that  initial  adsorption  results  in  a  small, 
but  perceptible,  rise  in  a.  It  seems  reasonable  to  suppose 
that  this  is  due  to  an  increase  in  diffuse  surface  scattering 
caused  by  a  small  increase  in  disorder  of  the  surface  due 
to  random  placement  of  the  hydrogen  before  the  c(2x2) 
order  becomes  long  range. 

VI.  SUMMARY 

We  have  investigated  the  behavior  of  changes  in  the  at¬ 
tenuation  coefficient  of  10-/im  SEW’s  caused  by  D2  and 
Hj  adsorption  on  W(100).  The  variations  have  been  dis¬ 
cussed  in  terms  of  changing  Drude  parameters  of  the  met¬ 
al  in  a  small  region  of  the  surface,  the  quenching  of  0.3- 
eV  intrinsic  surface  states,  and  surface-assisted  free- 
carrier  absorption.  Changes  in  the  Drude  parameters  of 
W  in  a  thin  layer  near  the  surface  are  generally  too  small 
to  account  for  the  size  of  the  observed  variations. 

For  0>O.42  ML  there  are  four  main  experimental 
points  which  contribute  to  an  interpretation  that  the  vari¬ 
ations  in  a  are  mainly  the  result  of  changes  in  the  specu¬ 
larity  of  scattered  free  carriers  caused  by  adsorbate- 
induced  reconstruction  of  the  W(  100)  surface  atoms.  The 
four  points  are  as  follows:  (1)  The  maximum  variation  in 
a  is  quite  large,  i.e.,  it  is  10  times  larger  than  the  WrH 
symmetric  stretch  contribution  to  a  (Ref.  27)  and  at  least 
2  times  larger  than  variations  observed  for  CO,  02,  and 
Nj  adsorption.  12)  A  coverage  dependence  for  a  is  ob¬ 
served  which  is  similar  to  resistivity  measurements  on 
gas-covered  films. 61,-71  (3)  There  exists  a  high  degree  of 
correlation  between  the  LEED  deduced  structures60  and 
the  variations.  (4)  The  induced  changes  are  broadband: 
they  are  observed  over  the  entire  range  of  frequencies 
studied  from  886  to  1088  cm'1,  and  they  are  fit  well  by 
the  theory  of  the  anomalous  skin  effect2  1  with  realistic 
values  of  changes  in  the  Fuchs  specularity  parameter.48 

In  this  coverage  range  our  results  show  that  the 
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9=  1.28  ML  disordered  surface  phase  and  the  saturation 
p(  1  x  1  )-D  phase  are  independent  of  temperature  between 
165  and  300  K.  Also,  the  previously  reported74  immobili¬ 
ty  of  the  adsorbate  below  200  K  seems  apparent  for  D2 
adsorption. 

In  the  coverage  range  of  0<O.42  ML  it  is  shown  that 
surface  scattering  is  not  the  only  mechanism  for  these 
broadband  changes.  Although  the  exact  contribution  of 
the  0.3-eV  intrinsic  surface  states  to  the  optical  response 
at  10  p.m  is  not  known,  it  seems  clear  that  they  play  some 
role.  We  suggest  that  surface  umklapp  scattering43-46 
may  have  some  contribution  to  the  large  negative  change 


in  a  as  0  approaches  0.42  ML.  in  addition,  the  surface 
umklapp  mechanism  and/or  a  changing  surface  Debye 
temperature  are  possible  explanations  for  the  0  =  0  ML 
temperature  dependence  of  A  a. 
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Changes  in  the  surface  electromagnetic  wave  attenuation  coefficient  a  upon  chemisorption  of 
N2,  Oj,  CO.  H2,  and  D2  onto  W(100)  have  been  measured  at  room  temperature  near  10  pm.  The 
size  as  well  as  the  coverage  dependence  of  the  variations  point  to  changes  in  the  amount  of  diffuse 
free  carrier  surface  scattering  as  the  dominant  a-aliering  mechanism. 


1.  Introduction 

In  1938  Fuchs  [1]  showed  that  the  DC  resistivity  of  a  thin  metallic  Film  is 
larger  if  the  conduction  electrons  colliding  with  the  surface  are  diffusely  rather 
than  specularly  reflected.  This  result  was  later  extended  into  the  infrared  (IR) 
by  Dingle  [2]  who  calculated  that  diffuse  surface  scattering  also  leads  to  higher 
absorption  by  a  metal  than  does  specular  reflection  of  the  charge  carriers. 
With  the  subsequent  advent  of  ultrahigh  vacuum  (UHV)  techniques  came  the 
experimental  possibility  of  controlling  the  disorder  of  the  surface  by  atomic  or 
molecular  adsorption  and  thereby  changing  the  amount  of  diffuse  electronic 
scattering  at  the  surface.  Since  then  there  have  been  numerous  investigations 
which  have  measured  adsorbate  induced  resistivity  variations  of  thin  films 
[3,4],  For  many  of  these  studies  it  has  been  demonstrated  that  the  Fuchs 
theory  applies  [3J.  Conspicuously  absent,  though,  are  equivalent  IR  measure¬ 
ments.  In  this  paper  we  report  the  results  of  such  infrared  measurements  on 
the  W(100)  surface. 

The  IR  absorption  technique  utilized  in  the  present  study  consists  of 
measuring  the  change  in  the  surface  electromagnetic  wave  (SEW)  attenuation 
coefficient  a  upon  chemisorption  of  N2,  02,  CO,  H2,  and  D2.  Broad-band 
changes  in  a  are  observed  upon  adsorption  of  each  of  these  gases.  The  size  and 
coverage  dependence  of  the  measured  variations  suggest  that  adsorbate  in¬ 
duced  changes  in  electronic  surface  scattering  provide  the  dominant  contribu¬ 
tion  to  the  observed  variations  in  a. 
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After  a  brief  description  of  the  experimental  details  in  section  2  the  results 
for  room  temperature  (RT)  adsorption  of  the  five  gases  are  presented  in 
section  3.  The  discussion  of  the  results  (section  4)  is  divided  into  two  main 
parts.  In  section  4.1  we  discuss,  rather  generally,  possible  contributions  to  the 
measured  changes  in  a  in  terms  of  the  three-layer  model  of  McIntyre  and 
Aspnes  [5],  Then  in  section  4.2  the  discussion  focuses  on  the  a  variations  for 
each  of  the  W(100)-adsorbate  systems  individually.  The  last  section  sum¬ 
marizes  our  results  on  adsorbate  induced  changes  in  the  SEW  attenuation 
coefficient. 


2.  Experimental 

Experimental  details  are  more  thoroughly  described  elsewhere  [6].  Briefly, 
radiation  in  the  10  pm  region  from  a  line-tunable  gas  laser  (either  nC02  or 
UC02)  is  coupled  to  SEW’s  via  a  grating  etched  in  the  sample  surface.  The 
SEW’s  travel  a  distance  L  =  5  cm  along  the  smooth  W(  100)- vacuum  interface 
before  another  grating  couples  them  back  to  free  space  radiation  which  is 
detected  with  a  HgCdTe  photoconductive  detector.  The  absorption  measure¬ 
ment  is  performed  by  monitoring  the  intensity  at  the  detector  as  gas  is 
adsorbed  onto  the  surface.  Any  change  in  a  is  calculated  from 

da  »  —L~x  ln(///0), 

where  10  is  the  zero-coverage  detected  intensity.  Initial  adsorption  pressures 
varied  between  5  X  10  ~ 9  and  10" 7  Torr  above  a  background  base  pressure  of 
typically  8  x  10" 11  Torr. 


3.  Results 

In  fig.  1  we  have  plotted  the  change  in  a  (da)  versus  exposure  for 
adsorption  of  the  five  different  molecules,  N2,  02,  CO,  H2,  and  D2,  onto 
W(100)  at  various  frequencies  near  10  pm.  The  exposure  is  determined  from 
uncorrected  ion  gauge  pressure  readings.  Initial  adsorption  results  in  an 
increase  in  attenuation  for  all  five  gases.  For  N2  da  reaches  a  plateau  at  an 
exposure  of  0.53  ±  0.03  langmuir  (L).  CO  and  02  adsorption  results  in 
qualitatively  similar  variations  in  a:  a  peak  at  some  intermediate  exposure 
[0.36  ±  0.07  (0:6  ±  0.1)  L  for  CO  (02)]  with  a  final  value  of  a  less  than  at  zero 
coverage.  Chemisorption  of  H2  and  D2  yields  much  richer  structure  in  da. 
After  the  small  initial  rise  a  sharp  dip  appears  at  0.100  ±  0.004  (0.140  ±  0.006) 
L  for  H2  (D2)  followed  by  a  broad  maximum  at  0.45  ±  0.04  (0.59  ±  0.05)  L. 
As  the  surface  becomes  saturated  a  decreases  to  a  value  which  is  still  larger 
than  for  the  clean  surface. 
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Exposure  (10-6  To rr  sec) 


Fig.  1.  Change  in  a  (Aa)  versus  exposure  for  N2,  02,  CO,  H2,  and  D2  adsorption  onto  W(100)  at 
room  temperature  for  typical  runs  near  1000  cm"1.  Laser  frequencies  used  for  runs  shown:  N2 
andOj,  1049  cm"1;  CO.  1088  cm"1;  H2  and  D2,  1032  cm"1. 


Fig.  2  shows  the  same  Aa  data  of  fig.  1  plotted  versus  coverage  instead  of 
exposure.  For  D2  and  H2  the  coverage  is  determined  using  our  thermal 
desorption  (TD)  data  which  indicate  that  the  sharp  dip  in  a  corresponds  to 
8  =  0.21  ±  0.02  ( 0  is  the  fraction  of  saturation  coverage)  along  with  the 
previous  results  that  saturation  corresponds  to  2  X  1015  atoms/cm2  [7]  and 
that  the  sticking  coefficient  a  varies  as  1  -  8  [8,9].  Flash  desorption  results  are 
also  used  to  obtain  the  CO  coverage.  At  room  temperature  CO  adsorbs  into 
two  classes  of  states  labeled  a  and  P  [10].  Saturated  P  corresponds  to 
5.0  X  1014  molecules/cm2  [11]  while  TD  data  show  that  for  RT  adsorption 
onto  our  sample  8a  =  0.58^  which  implies  8M  =  7.5  X  10’4  molecules/cm2.  By 
assuming  the  form  of  a(8)  as  measured  by  Clavenna  and  Schmidt  [10]  the 
coverage  versus  exposure  is  determined.  Both  the  H2-D2  and  CO  TD  results 
indicate  that  during  exposure  the  pressure  at  the  sample  is  higher  than  at  the 
ionization  gauge  by  a  factor  of  -  1.4.  Using  this  fact,  previously  determined 
forms  of  a(8)  for  N2  [12]  and  02  [13],  and  saturation  coverage  values  of 
3  X  1014  and  6  X  1014  molecules/cm2  for  N2  [12,14]  and  02  [13]  respectively, 
coverage  versus  measured  exposure  is  calculated  for  these  two  gases. 

Since  the  two  classes  of  CO  states,  a  and  /?,  are  easily  identifiable  in  TD 
spectra,  an  attempt  was  made  to  try  and  separate  out  the  contribution  to  Ja 
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Fig.  2.  Change  in  a  (da)  versus  molecular  coverage  for  the  same  data  shown  ui  fig.  1. 


Exposure  (10  Torr  sec) 

Fig.  3.  Change  in  a  (da)  versus  CO  exposure  at  room  temperature  and  1088  cm"1  showing  the 
effects  of  a-CO  desorption  and  readsorption,  da  upon  CO  adsorption  onto  a  clean  surface  ts 
shown  in  curve  a.  At  6.6  L,  CO  exposure  was  slopped,  the  sample  flashed  to  -  700  K  and 
recooled  to  RT,  resulting  in  the  drop  in  da  shown  as  the  discontinuity  between  curves  a  and  b. 
CO  was  then  re-admitted  to  the  UHV  chamber  where  upon  da  decreased  as  shown  in  curve  b.  At 
12.3  L  exposure  was  again  stopped,  the  sample  thermally  cycled  to  ~  700  K  and  back  to  RT.  and 
then  re-exposed  to  CO  resulting  in  curve  c. 
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from  each  class  as  shown  in  fig.  3.  After  saturating  the  clean  surface  with  CO 
at  room  temperature  (curve  a)  the  a-CO  was  thermally  desorbed  ( TmiX  ~  700 
K)  and  the  sample  recooled  to  RT  whereupon  it  was  re-exposed  to  CO,  giving 
the  change  in  attenuation  shown  in  curve  b.  The  sample  was  again  thermally 
cycled  to  ~  700  K  and  back  to  RT  and  re-exposed  to  CO,  resulting  in  a 
change  in  attenuation  essentially  identical  to  curve  b  as  shown  in  curve  c. 


4.  Discussion 

4. 1.  Contributions  to  Aa 

For  a  metal  with  chemisorbed  molecules  the  standard  model  used  to 
describe  the  IR  response  is  the  three-layer  model  of  McIntyre  and  Aspnes  [5] 
which  describes  the  surface  region  between  vacuum  and  the  bulk  metal  by  an 
isotropic,  uniform  dielectric  layer.  The  effect  of  the  adsorbed  molecules  on  the 
IR  response,  then,  is  either  to  alter  the  dielectric  function  td  of  the  layer  by 
changing  the  vibrational  or  electronic  structure  at  the  surface,  or  to  change  the 
effective  dielectric  function  «m  of  the  metal  by  changing  the  nature  of  the  free 
carrier  surface  scattering. 

4.1.1.  Changes  in  id 

Several  possible  mechanisms  for  changing  «d  can  be  ruled  out.  Previous  IR 
and  electron  energy  loss  spectroscopy  (EELS)  measurements  eliminate  surface 
vibrations  as  a  contribution  to  the  strong  variations  in  a.  EELS  results  for  CO 
(15,16)  and  N2  (17)  adsorption  show  that  there  are  no  dipole  allowed  vibra¬ 
tional  transitions  in  the  1000  cm-1  region.  H2  adsorption  does  produce  an  IR 
allowed  vibrational  mode  in  this  region  [18,19],  but  it  has  been  previously 
shown  [19]  that  the  maximum  saturated-coverage  attenuation  due  to  this  mode 
[the  W2-H  symmetric  stretch  (iq)]  is  only  6  X  10“ 3  cm-1.  (This  accounts  for 
most  of  the  stronger  H2  adsorption  at  saturation  as  compared  to  D2.)  Oxygen 
also  produces  a  mode  at  ~  1000  cm-1,  but  EELS  studies  [20]  show  for  RT 
adsorption  that  exposures  greater  than  those  used  here  are  required  before  this 
vibration  becomes  apparent. 

Adsorption  can  also  be  expected  to  change  td  by  altering  the  bulk  IR 
parameters  near  the  surface.  In  general  -  and  especially  for  H:  and  D2  - 
possible  a  variations  resulting  from  changes  in  free  carrier  density  and  bulk 
relaxation  time  are  too  small  to  account  for  the  measured  changes  [6],  The 
same  is  also  true  for  changes  in  high  frequency  interband  structure  near  the 
surface  [6]. 

Contributions  to  which  cannot  be  ruled  out  are  variations  in  the  IR 
response  due  to  low  frequency  (mid-IR)  electronic  structure  of  the  surface 
upon  adsorption.  The  dominant  low  frequency  feature  on  W(100)  is  a  peak  in 
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the  surface  density  of  states  -  0.3  eV  below  the  Fermi  level.  Both  field 
emission  energy  distribution  (FEED)  [21,22]  and  ultraviolet  photoemission 
spectroscopy  studies  [23,24]  have  detected  the  presence  of  these  states  and 
have  shown  that  they  disappear  upon  chemisorption  of  each  of  the  molecules 
studied  here.  Although  0.3  eV  is  rather  far  above  0.125  eV  (1000  cm'1)  both 
FEED  [22]  and  surface  reflection  spectroscopy  (SRS)  results  [25]  indicate  that 
these  states  are  rather  broad.  Further,  the  SRS  data  show  that  final  state 
energies  for  transitions  from  surface  states  on  W(100)  lie  just  above  the  Fermi 
level.  It  is  therefore  reasonable  to  assume  that  they  contribute  to  optical 
absorption  in  the  1000  cm'1  frequency  region.  Adsorbate  induced  quenching 
of  these  states,  and  hence  their  oscillator  strength,  may  account  for  the 
negative  value  of  Aa  at  saturation  for  02  and  CO,  as  well  as  the  large  dip  in 
Aa  for  H2  and  D2.  In  fact,  the  frequency  dependence  of  the  D2  induced 
changes  in  a,  which  have  been  measured  over  the  frequency  range  of  the  two 
C02  lasers  (886  to  1088  cm'1),  strongly  support  this  idea  [6]. 

Aside  from  these  intrinsic  states  it  is  not  clear  what  contribution  other  low 
frequency  surface  structure  has  on  the  SEW  attenuation.  The  lack  of  calcula¬ 
tions  on  the  low  frequency  optical  response  of  these  systems  and  the  rather 
limited  possible  frequency  range  of  the  present  study  disallows  any  answer  to 
this  question  at  present. 

4.1.2.  Changes  in  i„ 

In  spite  of  the  uncertainty  of  the  contribution  to  Aa  from  Aia.  variations  in 
the  effective  bulk  dielectric  constant  cm  caused  by  changes  in  diffuse  scatter¬ 
ing  of  the  bulk  free  carriers  offer  an  explanation  for  the  size  of  Aa  as  well  as 
the  features  in  a($). 

For  W  near  10  pm  the  size  of  the  surface  scattering  contribution  can  be 
estimated  using  Holstein’s  phenomenological  expression  [26] 

-•.-'-(i-wKtK  (i) 

for  the  effective  scattering  rate  due  to  collisions  with  the  surface  along  with  the 
Drude  model  formula  [27] 

a  -  (w/wp^cr)-1  (2) 

for  the  free  carrier  contribution  to  the  SEW  attenuation  coefficient.  In  these 
equations  u  is  the  frequency,  wp  the  plasma  frequency  (both  in  rad/s),  vF  the 
Fermi  velocity,  c  the  speed  of  light,  1/r  the  free  carrier  relaxation  rate,  and  p 
the  Fuchs  specularity  parameter  [1]  which  characterizes  the  effective  fraction 
of  electrons  which  are  reflected  specularly  from  the  surface.  The  use  of  the 
appropriate  IR  parameters  for  W  [6]  and  the  replacement  of  1/r  in  eq.  (2)  by 
l/r„  [eq.  (1)]  yields  the  surface  scattering  contribution  to  the  attenuation: 
a  -0.16(1  -  p)  cm'1.  The  largest  measured  excursion  in  a.  which  occurs  for 
H2  between  2.1  and  -  6  X  1014  molecules/cm2,  corresponds  to  a  change  in 
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the  specularity  parameter  of  -0.4,  a  value  which  is  well  within  the  maximum 
variation  of  \Ap  |  *  1. 

Diffuse  surface  scattering  is  also  suggested  by  the  coverage  dependence  of  a 
for  adsorption  of  the  five  molecules.  The  signature  of  Aa  versus  CO,  02,  and 
N2  coverage  is  similar  to  adsorbate  induced  resistivity  variations  of  thin  films 
(either  a  monotonic  increase  or  a  maximum  at  some  intermediate  coverage); 
for  a  number  of  the  resistivity  studies  Wissmann  [3]  has  shown  that  the  Fuchs 
diffuse  scattering  theory  applies.  The  same  signature  is  also  apparent  for  H2 
and  D2  if  the  changes  before  the  sharp  dip  are  ignored. 

These  observations  on  both  the  size  and  coverage  dependence  of  Aa  lead  us 
to  consider  this  interpretation  of  the  changes,  in  light  of  other  experimental 
evidence,  for  each  of  the  W(100)-adsorbate  systems  measured  here. 

4.2.  Specific  W(100)-adsorbate  systems 

4.2.1.  W(100)~CO 

Among  the  systems  studied  here  W(100)-CO  has  the  possibility  of  provid¬ 
ing  the  simplest  example  of  adsorbate  induced  changes  in  diffuse  surface 
scattering  since  no  ordered  superstructures  form  from  RT  adsorption.  Chem¬ 
isorption  of  CO  results  in  approximately  sequential  filling  of  the  fi  and  then 
the  a  states.  fi-CO  is  thought  to  be  dissociated  [15]  (although  strong  remnant 
C-O  interactions  cannot  be  ruled  out  [16])  with  the  C  and  O  atoms  occupying 
separate  four-fold  hollow  sites.  o-CO  is  molecular  with  the  C  bonded  to  the 
substrate  in  the  on-top  position  [15]. 

Initial  adsorption  of  CO  (as  well  as  all  of  the  other  gases)  results  in  an 
increase  in  Aa  suggesting  that  initially  the  C  and  O  atoms  adsorb  randomly 
resulting  in  a  disordered  surface  and  hence  an  increase  in  the  diffuse  scattering 
of  the  free  carriers.  Such  disorder  has  also  been  inferred  from  low  energy 
electron  diffraction  (LEED)  studies  [28]  which,  upon  CO  adsorption,  initially 
show  an  increase  in  the  diffusely  scattered  background. 

With  adsorption  beyond  -lx  10u  cm-2  (henceforth  all  coverages  are 
molecular  unless  otherwise  stated)  Aa  decreases.  There  are  two  apparent 
possibilities  whose  relative  contributions  cannot  be  distinguished.  The  first  is 
that  adsorption  beyond  this  point  results  in  an  overlayer  which  is  more 
ordered,  resulting  in  a  peak  in  Aa  at  an  intermediate  coverage  [3,29]  (similar 
to  a  binary  alloy  AXB,_*  where  the  resistivity  shows  a  maximum  at  an 
intermediate  value  of  x  between  0  and  1  [30]).  The  other  possibility  is  that  CO 
induced  quenching  of  the  0.3  eV  intrinsic  surface  states  begins  to  dominate  the 
IR  response  at  this  point.  If,  however,  the  magnitude  of  the  quenching  is  a 
constant  per  adsorbed  molecule  (a  reasonable  assumption,  see  ref.  [21])  and 
only  these  two  processes  contribute  to  *ia.  then  it  can  be  concluded  that 
adsorption  beyond  1  x  10u  cm”2  results  in  less  disorder  per  adsorbed  mole¬ 
cule  than  initially  occurs. 
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Because  a-CO  is  top  bonded  and  therefore  further  away  from  the  W 
surface  it  should  have  less  of  an  effect  on  electronic  surface  scattering  than 
/8-CO.  The  fact  that  Aa  is  constant  for  coverages  >  5  x  1014  cm-2  lends  some 
support  to  this  idea,  although  the  data  of  fig.  3  show  that  a-CO  adsorption 
does  induce  some  change  in  the  SEW  attenuation.  Desorption  of  a-CO  and 
recooling  to  room  temperature  results  in  a  decrease  in  Aa  of  -  6.5  X  10“ 3 
cm-1  as  signified  by  the  jump  between  curves  a  and  b.  Readsorption  of  a-CO. 
however,  causes  a  further  Aa  decrease,  rather  than  an  increase  back  to  the 
previous  saturated  coverage  value.  Another  RT-700  K-RT  thermal  cycle  and 
readsorption  (curve  c)  results  in  Aa  values  essentially  the  same  as  for  curve  b. 

Two  conclusions  are  suggested  by  these  last  observations.  The  first  is  that 
thermal  annealing  of  /J-CO  to  -  700  K  results  in  an  increase  in  order  of  this 
layer.  Previous  work  function  measurements  [11]  have  shown  that  some  kind 
of  structural  rearrangement  takes  place  within  the  saturated  fi-CO  layer  upon 
heating  to  700  K.  That  this  rearrangement  results  in  more  order  is  not  too 
surprising  since  it  is  well  known  that  further  annealing  (to  -  1150  K)  results 
in  a  well  ordered  c(2  X  2)  superstructure  upon  desorption  of  half  the  /?-CO 
[28].  The  second  conclusion  is  that  a-CO  has  no  effect  on  the  amount  of 
diffuse  surface  scattering  since  initial  adsorption  into  any  state  should  increase 
the  disorder  of  the  surface  (assuming  that  no  simultaneous  reconstruction  of 
the  surface  occurs). 

Although  the  effects  of  CO  adsorption  are  not  as  simple  as  might  have  been 
expected,  they  do  point  to  surface  assisted  absorption  caused  by  diffuse 
surface  scattering  as  a  contribution.  The  main  support  comes  from  the  6 
dependence  of  J8-CO  adsorption  and  the  fact  that  initial  a-CO  adsorption 
decreases  the  attenuation  of  the  surface  wave. 

4.2.2.  W(100)-0 

The  similar  qualitative  changes  in  Aa  for  02  and  CO  are  at  first  surprising 
since  LEED  [31]  and  EELS  [20]  indicate  that  02  produces  ordered  superstruc¬ 
tures  on  W(100).  Initially,  diffuse  LEED  streaks  arc  formed  along  the  {10} 
directions  as  02  is  adsorbed.  As  the  coverage  approaches  2.5  x  10u  cm"2  the 
diffuse  streaks  break  apart  into  a  p(4x  1)  pattern.  With  further  coverage  a 
complicated  series  of  LEED  patterns  is  formed  which,  along  with  EELS, 
suggests  that  the  oxygen  atoms  are  serially  forming  double  (at  2.5  X  10M 
cm-2),  triple  (at  3.75  X  1014  cm"2),  and  quadruple  (at  5.0  X  1014  cm"2)  rows. 

As  expected,  Aa  increases  as  oxygen  initially  bonds  to  the  surface;  this 
continues  until  a  coverage  of  2.5  X  1014  cm"2.  Even  though  a  well  ordered 
structure  forms  at  this  coverage  it  is  not  inconsistent  that  there  exists  a  peak  in 
Aa  due  to  surface  scattering.  This  is  because  free  carrier  reflections  equivalent 
to  ([2«  +  l]/2,  0)  and  ([2«  +  l]/4,  0)  LEED  spots  are  effectively  diffuse  since 
such  reflections  result  in  nonzero  changes  in  the  conduction  electron  velocity 
parallel  to  the  metal-vacuum  interface  [6.26].  The  decrease  in  Aa  after  this 
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coverage  indicates  that  the  effective  disorder  per  adsorbed  atom  is  decreasing 
up  to  -  5  X  1014  cm-2,  probably  due  to  less  scattering  into  these  effectively 
diffuse  LEED  channels. 

A  comparison  of  the  CO  and  02  results  is  interesting.  Initially,  both 
molecules  are  thought  to  break  apart  with  individual  atoms  occupying  the 
four-fold  hollow  sites  [IS, 20].  If  this  is  the  case  then  our  results  indicate  that 
the  carbon  atoms  interact  more  strongly  with  the  free  carriers  than  do  the 
oxygen  atoms.  In  addition,  the  fact  that  nearly  saturated  02  coverages,  which 
exhibit  pt.4  x  l)-like  structures,  have  a  higher  value  of  Aa  than  comparable  CO 
coverages  further  implies  that  the  ([2n  +  l]/2,  0)  and  ([2n  +  l]/4,  0)  reflec¬ 
tions  are  effectively  diffuse. 

So,  02  adsorption,  in  spite  of  the  complicated  superstructures  which  form, 
affects  da  in  a  manner  similar  to  what  would  be  expected  for  a  simpler 
order-disorder-order  type  of  coverage  dependence. 

4.2.3.  W(100)-H  and  W(100)-D 

Although  H2  and  D2  adsorption  result  in  the  most  complicated  structure  in 
a(9),  these  two  W(100)-adsorbate  systems  provide  the  best  examples  of 
adsorbate  induced  changes  in  electronic  surface  scattering  among  those  studied 
here.  Partly  this  is  due  to  the  fact  that  a  fairly  small  fractional  coverage 
(0  =  0.2)  of  either  adsorbate  (as  compared  to  -  saturated  coverages  of  02  and 
N2,  and  -  half-saturated  coverage  of  CO)  is  sufficient  to  remove  the  0.3  eV 
intrinsic  surface  states  [21-25].  Because  of  this  we  concentrate  on  changes  in  a 
after  the  sharp  minimum  at  2.1  X  1014  cm~2  (0  *  0.21). 

The  coverage  dependence  of  a  between  0  =  0.21  and  2.0  is  again  suggestive 
of  changes  in  diffuse  surface  scattering  as  the  dominant  a-changing  mecha¬ 
nism.  This  is  corroborated  by  other  experimental  results  on  this  system.  At 
0  =  0.21,  where  a  is  a  minimum,  the  surface  is  reconstructed  in  what  is  termed 
the  incommensurate-c(2  x  2)-H  (or  -D)  phase  which  is  characterized  by  1/2 
order  spots  split  into  a  quartet  along  the  ±(10)  and  ±(01)  directions  [9,32]. 
As  more  gas  is  added  the  spots  of  the  1/2  order  quartet  streak  along  directions 
perpendicular  to  their  splitting  t  id  then  eventually  disappear,  indicating  a  loss 
of  order  in  one  dimension  followed  by  a  loss  of  order  in  both  dimensions. 
These  LEED  changes,  signalling  the  onset  of  disorder,  are  accompanied  by  a 
sharp  rise  in  the  SEW  attenuation  until  the  maximum  at  0  =  0.6.  The  subse¬ 
quent  fall  in  a  as  the  surface  approaches  the  well-ordered  p(l  X  1)-H  state  at 
saturation  occurs  concurrently  with  a  decrease  in  the  diffuse  LEED  back¬ 
ground.  The  high  degree  of  order  at  saturation  is  also  reflected  in  high  energy 
ion  scattering  results  [7]  which  show  that  the  top  layer  is  unreconstructed  and 
in  IR  vibrational  spectra  [33]  which  show  that  the  v,  and  2v2  fi}  W2~H 
vibrational  absorption  features  are  their  sharpest  at  6  =  1. 

A  comparison  with  02  and  CO  shows  that  Aa  for  H2  and  D2  is  at  least 
twice  as  large.  This  seems  strange  since  it  is  expected  that  the  larger  atoms 
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should  have  a  higher  cross  section  for  electronic  scattering  and  suggests  that  it 
is  the  direct  interaction  of  the  free  carriers  with  the  reconstructed  top  W  layer 
by  H  or  D  which  results  in  the  changes  in  a  for  these  two  adsorbates. 

The  frequency  dependence  of  Aa  after  6  =  0.21,  which  has  been  measured 
between  886  and  1088  cm'1,  also  supports  the  surface  scattering  thesis.  We 
find  that  Aa  between  6  =  0.21  and  0.6  and  between  0  =  0.6  and  2.0  are  fit  well 
using  the  expressions  calculated  by  Dingle  (2). 

For  6  <  0.21  it  is  not  clear  what  the  sharp  dip  in  Aa  is  due  to.  Contrary  to 
the  region  for  6  >  0.21  the  frequency  dependence  of  the  a  changes  for  lower 
coverages  does  not  comform  to  Dingle’s  expressions,  but  instead  point  to  the 
quenching  of  the  intrinsic  surface  states  as  providing  some  contribution. 
Surface  disorder  is  largely  ruled  out  by  LEED  observations  [32]  which  show 
long  range  c(2  X  2)  order  for  these  coverages.  However,  unlike  the  changes  in 
crystal  momentum  associated  with  the  p(4  x  l)-0  structures,  the  c(2  X  2)  1/2 
order  LEED  reflections  connect  relatively  few  states  on  the  W  Fermi  surface 
and  therefore  cannot  directly  contribute  in  a  large  way  to  changing  the 
character  of  the  free  carrier  scattering  from  the  surface.  An  interesting 
possibility,  discussed  in  detail  elsewhere  [6],  is  that  a  decrease  in  relaxation  of 
the  top  layer  W  atoms  as  0  approaches  0.2  [9]  induces  changes  in  free  carrier 
interband  scattering  at  the  surface,  thus  altering  the  1R  absorption. 

So,  for  W(100)-H  (and  -D)  Aa  after  0  =  0.21  is  well  explained  by  the 
surface  assisted  absorption  process.  For  0  <  0.21  the  removal  of  the  intrinsic 
surface  states  seems  to  play  a  role  although  it  is  not  clear  that  this  is  the  only 
a-altering  mechanism  operating  in  this  region. 

4.2.4.  W(I00)-N 

The  simple  behavior  of  a  with  N  coverage  is  a  bit  puzzling  in  light  of  the 
data  for  the  other  four  adsorbates.  At  all  stages  of  RT  adsorption  c(2  x  2) 
structure  is  evident  in  LEED  measurements  [14,34],  Analysis  of  the  LEED 
data  has  resulted  in  the  interpretation  that  the  N  atoms  are  bound  in  the 
four-fold  hollow  sites  and  that  the  atoms  are  grouped  into  islands  [34], 

This  island  formation  may  be  the  key  to  the  difference  between  N: 
adsorption  and  that  of  the  other  gases.  As  shown  in  figs.  1  and  2  Aa  is  strictly 
increasing  for  N2  while  all  the  other  gases  have  regions  of  negative  slope  as 
well  as  negative  Aa.  This  difference  is  somewhat  surprising  since  N2  adsorp¬ 
tion  is  also  effective  in  removing  the  W  intrinsic  surface  states  [21,24],  Since 
the  c(2  x  2)  1/2  order  spots  should  contribute  negligibly  to  effective  diffuse 
scattering,  as  mentioned  above,  this  suggests  that  diffuse  scattering  from  the 
island  edges  is  quite  strong  -  strong  enough,  in  fact,  to  cancel  any  negative  Aa 
tendencies  caused  by  intrinsic  surface  state  quenching. 

A  recently  proposed  model  [35]  for  N2  adsorption  onto  W(100)  offers  some 
explanation  for  why  scattering  at  the  island  edges  should  be  so  strong.  The 
model  involves  displacements  of  the  top  layer  W  atoms  with  the  largest 
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displacements  occurring  at  the  edges  of  the  islands.  As  in  the  case  of  H2  and 
D2,  the  interaction  of  the  free  carriers  with  the  reconstructed  W  atoms  is 
expected  to  be  quite  large,  and  thus  could  negate  the  loss  of  oscillator  strength 
from  the  surface  states. 

A  possible  problem  with  this  explanation  is  that  beyond  some  coverage  Aa 
should  still  decrease  somewhat  since  the  total  length  of  the  island  boundaries 
should  decrease  as  the  islands  begin  to  coalesce.  Although  not  apparent  in  the 
data  of  figs.  1  and  2  some  N2  runs  do  seem  to  show  a  slight  decrease  in  a 
beyond  -  2x  10u  cm';  before  the  plateau,  but  the  measured  decrease  is 
always  less  than  the  long  term  fluctuations  (=  2  X  10“ 3  cm-')  in  the  detected 
signal.  LEED  measurements  do  indicate,  though,  that  such  coalescence  is  not 
very  pronounced  unless  the  sample  is  annealed  above  room  temperature. 
(Although  it  has  not  yet  been  studied,  the  effect  on  a  of  such  annealing  should 
provide  a  good  test  of  this  idea  of  strong  scattering  from  the  island  edges.)  In 
addition,  at  the  highest  coverages  there  exist  1/5  order  features  in  the  ±(10) 
and  ±(01)  directions  around  the  integral  order  beams  which  should  contribute 
to  effectively  diffuse  scattering. 

For  N2,  then,  we  propose  that  strong  scattering  from  the  displaced  W 
atoms  on  the  boundaries  of  the  N-induced  islands  is  responsible  for  the 
increase  in  diffuse  scattering  upon  adsorption. 


5.  Summary 

Changes  in  the  SEW  attenuation  coefficient  a  upon  chemisorption  of  N2, 
02,  CO,  H2,  and  D2  have  been  presented  and  discussed.  Although  the 
adsorption  behavior  of  each  species  on  the  W  surface  is  quite  different  (N2, 
island  formation;  02,  row  formation;  CO,  no  superstructures;  H2  and  D2,  W 
reconstruction),  the  changes  in  a  can  be  fairly  well  understood,  in  large  part, 
in  terms  of  changes  in  the  diffuse  scattering  of  the  free  carriers  of  the  metal. 

A  useful  extension  of  the  present  work  would  involve  measuring  the 
adsorbat*e  induced  1R  changes  over  a  broader  frequency  range,  both  higher 
and  lower  in  energy.  Such  a  study,  which  would  necessitate  the  use  of  a  more 
broad- banc,  albeit  less  sensitive,  technique  such  as  IR  absorption  spec¬ 
troscopy,  should  help  in  sorting  out  the  importance  of  the  different  contribu¬ 
tions,  such  as  surface  states  and  surface  scattering,  to  the  broad-band  IR 
response  of  these  W(100)-adsorbate  systems. 
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Cyclotron  resonance  of  a  two-dimensional  electron  gas  at  a  GaAs/GaAIAs  interface  is 
measured  in  the  far  infrared  at  intensities  of  up  to  10  kW/cm’.  Both  the  cyclotron  mass  and 
the  carrier  density  are  independent  of  intensity  but  the  relaxation  time  changes  with  a 
dependence  which  is  similar  to  that  observed  earlier  in  high  dc  field  mobility  studies. 


The  high  electron  mobility  in  a  two-dimensional  elec¬ 
tron  gas  ( 2DEG )  system  decreases  at  large  dc  electric  fields 
because  of  field  induced  electron  heating  above  the  lattice 
temperature.'"  When  a  magnetic  field  is  applied  perpendic¬ 
ular  to  the  conducting  plane,  the  high  mobility  translates 
into  a  narrow  resonance  line  shape  centered  at  the  cyclotron 
frequency '  and  the  nonlinear  properties  of  the  2DEG  can  be 
probed  in  more  detail  with  intense  far-infrared  (FIR)  radi¬ 
ation.  For  example,  with  circularly  polarized  radiation 
tuned  to  the  resonance  condition,  the  analogue  of  the  dc 
electric  field  experiment  can  be  carried  out.  In  addition,  by 
mapping  out  the  complete  line  shape  at  different  intensities, 
the  electric  field  dependence  of  the  effective  carrier  density, 
the  cyclotron  mass,  and  the  relaxation  time  can  be  measured 
to  provide  a  more  complete  picture  of  the  underlying  dy¬ 
namics. 

In  this  letter  we  report  on  the  intensity  dependence  of 
the  cyclotron  resonance  transmissivity  profile  of  a  high  mo¬ 
bility  (5x10s  cm:/V  s)  GaAs/Ga,,  7  Al<j  j  As  2DEG  at 
three  FIR  frequencies:  8.2, 27.8,  and  53.4  cm" '.  For  all  sam¬ 
ples,  the  resonant  field  to  zero  field  transmissivity  ratio 
( termed  normalized  transmissivity )  is  observed  to  increase 
for  laser  intensities  above  10  W/cm!  and  approach  1  at  the 
highest  intensities  of  10  kW/cnr.  OfT  resonance,  in  the 
wings  of  the  cyclotron  mode,  the  transmissivity  decreases 
with  increasing  intensity.  Fitting  these  cyclotron  resonance 
line  shapes  using  the  Drude  conductivity  relation  shows  that 
the  dramatic  change  in  the  transmissivity  occurs  because  of 
FIR  field  induced  changes  in  the  cyclotron  resonance  relaxa¬ 
tion  time. 

The  experimental  arrangement  is  shown  schematically 
in  Fig.  1.  A  high  power  pulsed  transversely  excited  atmo¬ 
spheric  (TEA)  CO,  laser  C  is  used  at  a  repetition  rate  of  0.2 
Hz  to  pump  the  superfluorcscent  FIR  transitions  in  CHjF 
gas  at  27.8  and  53.4  cm  * 1  and  in  '  'CH,F  at  8.2  cm  '.  The 
gas  is  enclosed  in  an  8-m-long  Pyrcx  tube  and  the  pressure 
optimized  to  obtain  maximum  FIR  power  at  each  frequency. 
FIR  pulses  of  40  ns  with  up  to  12  mJ  of  energy  are  obtained  4 
The  output  FIR  laser  pulse  polarization  is  fixed  by  a  metal 
wire  polarizer  P  About  109?-  of  the  FIR  laser  output  reflect¬ 
ed  by  a  Mylar  beamsplitter  B  is  focused  by  a  parabolic  mirror 
onto  a  py  roelectric  detector  R  as  a  reference  signal,  and  the 
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remainder  is  sent  through  an  oversized  stainless-steel  wave¬ 
guide  W  to  the  polymethylpentene  (TPX)3  cryostat  input 
window  Q.  Inside  the  liquid  helium  cryostat  H  the  radiation 
is  focused  through  a  linear  polarizer  onto  the  sample  by  a 
TPX  lens,  recollimated  after  the  sample  by  a  second  TPX 
lens  and  reflected  back  to  the  cryostat  output  window  by  a 
copper  mirror.  The  energy  detector  D  is  either  a  Golay  cell 
or  a  Ge:Ga  bolometer  operated  at  2  K.  The  detector  signal  is 
digitized  and  read  by  a  computer.  To  obtain  good  statistics 
20  shots  in  which  the  reference  signal  does  not  vary  by  more 
than  20%  are  taken  at  each  magnetic  field  value  and  zero 
field  normalizing  measurements  are  made  periodically  to 
cancel  out  slow  changes  in  detector  sensitivity.  The  pulse 
energy  at  the  sample  position,  measured  with  a  pyroelectric 
detector,  divided  by  the  area  of  the  aperture  and  by  the  pulse 


FIG  I  Experimental  arrangement.  Legend:  M.  mirror;  C.  TEA  CO.  laser; 
/.  inpul  NaCI  window.  F,  FIR  laser;  O.  crystal  quartz  and  TPX  oulpul  win¬ 
dow;  4,  digitizer  and  microcomputer;  R.  reference  pyrometer  deteclof;  P. 
linear  polarizing  grid.  B.  Mylar  beamsplitter;  T.  metal  mesh  attenuators;  D. 
Ge  Ga  bolometer  or  Golay  cell  detector,  W,  oversize  waveguide;  Q.  TPX 
cryostat  window;  II  liquid  helium  cryostat,  L,  TPX  lens;  S.  superconduct¬ 
ing  magnet  and  2  DEG  sample. 
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duration  is  the  average  incident  intensity  quoted  here.  Al¬ 
though  it  is  difficult  to  estimate  reliably  the  degree  of  lattice 
heating  in  the  2DEG  sample  since  both  the  nature  of  the 
electron-phonon  coupling  and  the  dissipation  of  heat  in  het¬ 
erostructure  in  a  magnetic  field  are  not  well  understood,  for 
a  given  power  the  shorter  pulse  length  the  smaller  will  be 
such  heating.  Our  FIR  pulses  are  an  order  of  magnitude 
shorter  than  those  used  in  the  high  dc  field  mobility  stud¬ 
ies.1' 

The  2DEG  sample  with  the  GaAs  substrate  wedged  at 
2*  to  minimize  interference  effects  is  mounted  in  a  4  T  super¬ 
conducting  magnet  in  the  Faraday  configuration.  Black 
polyethylene  on  each  side  of  the  sample  compartment  elimi¬ 
nates  room  light.  The  electron  density  is  monitored  by  Shub- 
nikov-de  Haas  resistivity  measurements  and  can  be  in¬ 
creased  by  activating  a  light-emitting  diode  near  the  sample. 
Both  the  magnet  and  the  sample  are  immersed  in  liquid  heli¬ 
um  at  4.2  K.  To  measure  the  intensity -dependent  transmissi¬ 
vity,  calibrated  metal  mesh  attenuators  7  are  placed  before 
and  after  the  waveguide  cryostat  system  so  that  the  detector 
signal  remains  about  the  same  throughout  the  experiment. 
In  this  manner  the  intensity  at  the  sample  can  be  varied  over 
five  orders  of  magnitude  without  saturating  the  detector. 
Two  types  of  experiments  have  been  performed:  ( 1 )  at  fixed 
FIR  laser  frequency  the  normalized  transmissivity  is  mea¬ 
sured  as  a  function  of  incident  intensity  and  (2)  at  fixed 
intensity  the  magnetic  field  is  scanned  stepwise  to  map  out 
the  cyclotron  resonance  line  shape. 

The  intensity  dependence  of  the  normalized  transmissi¬ 
vity  at  the  cyclotron  resonance  condition  has  been  studied  at 
three  FIR  frequencies:  8.2, 27.8,  and  53  4  cm  ’.  The  experi¬ 
mental  results  are  shown  in  Fig.  2.  At  8.2  and  53.4  cm  1  the 
density  is  2.5  X  10"  electrons/cm'.  At  27.8  cm  1  three  dif¬ 
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FIG  2  Normalized  transmissivity  «  intensity  at  three  different  laser  fre¬ 
quencies  (at  53  4cm  '.  2.5  *  10"  elcctrons/cm  tb)  27  8cm  1  squares. 
15*  10"  electrons/cm  .  circles.  2  5  ■  10"  elecirons/cm  .  triangles. 
4  5x  10"  electrons/cm  (cl  8  2  cm  2  5  *  10"  electrons/cm  Typical 
esperimental  error  is  represented  by  (he  error  bar  in  ic) 


ferent  densities  have  been  studied:  1.5,  2.5.  and  4.5  -  10" 
electrons/cnr.  Within  the  uncertainty  of  the  measurement 
the  transmissivity  begins  to  rise  for  all  curves  in  Fig.  2  at  the 
same  intensity,  independent  of  laser  frequency  and  electron 
density.  When  the  magnetic  field  is  scanned  at  a  fixed  inten¬ 
sity  to  map  out  the  cyclotron  resonance  line  shape,  results 
like  those  shown  in  Fig.  3  are  obtained  at  the  three  frequen¬ 
cies.  The  cyclotron  mass  is  essentially  independent  of  laser 
intensity. 

In  order  to  quantify  these  experimental  results  the  data 
as  typified  by  Fig.  3  are  fitted  to  the  frequency-dependent 
transmissivity  equation  for  linearly  polarized  light1’ 


/ha,.)  =2/7: 


_ I _ 

1  +  4-/7  (t .  ( t-j  i/c 


i  I  4-  4,7/7  (T  ( (■>  i  /<-  j 

where/?  -  ( 1  -t-  fl(ijA, )  being  the  refractive  index 

of  the  GaAs  substrate  and  c  being  the  speed  of  light.  The 
2 DEG  conductivity  a  (pzt  )  is  assumed  to  be  of  the  Drude 
form,  namely. 


tt.rr  I 

<z,  (*>,)= - - - .  (2) 

m*  1  -  t((a  -t-u,  )r 

where  n,  is  the  2DEG  electron  density,  e  is  the  electron 
charge,  m*  is  the  electron  effective  mass,  r  is  the  electron 
relaxation  time,  and  ra,  =  eB  /m*c  is  the  cyclotron  frequen¬ 
cy  for  the  magnetic  field  B.  The  values  n0aA,  =  3.5  and 
m*  =  0.067 m,  are  used  in  the  calculation. 


HG  .V  Variation  with  FIR  intensity  of  the  2DEG  cyclotron  resonance 
transmissivity  profile  (a)  0  1  W/cm',  (b)  0  6  W7cm\  (c )  4.3  W/cm\  { d  > 
20  W/em  .  le)  74  W'/cm  .  I  f)  360  W/cnv  Laser  frequency  is  27  8  cm  *. 
temperature  is  4  2  K.  density  is  2  5  x  10"  elect rons/cm 
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FIG.  4.  Two-dimensional  electron  gas  density  n,  and  relaxation  time  rvs 
FIR  intensity  at  the  sample.  The  values  are  obtained  by  a  least-square  lit  of 
the  intensity-dependent  line  shape  data,  e  g.,  Fig.  3,  to  the  transmissivity 
equation  [  Eq.  ( I )  J .  For  each  set  of  data  n,  and  r  are  normalized  by  their 
values  at  low  intensity.  Open  squares:  =  2.3  x  10",  r0=  11  ps,  8.2 

cm  filled  circles:  n„  =  2.9x  10"  electrons/cm:,  r„  =  11  ps.  27.8  cm-1; 
open  circles:  =  3.0X  10".  r„  =13  ps.  27.8  cm'1;  triangles: 

=2.2x10"  electrons/cm:,  r„  =  9  ps.  27.8  cm'1;  filled  squares: 
Hu  =  2.6 X  10"  electrons/cm',  ru=  8  ps.  53.4  cm'1.  Solid  and  dashed 
lines:  Refs.  I  and  2.  respectively. 


At  each  laser  intensity  the  two  fit  parameters  n,  and  r 
are  determined.  The  solid  lines  in  Fig.  3  show  the  best  fits  of 
this  particular  set  of  data  to  Eq.  ( 1 ).  The  low  power  values 
for  the  electron  density  n and  relaxation  time  t0  obtained 
from  the  fit  compare  well  with  the  Shubnikov-de  Haas  de¬ 
rived  electron  density  and  with  the  electron  relaxation  time 
derived  from  the  measured  dc  mobility.  Figure  4  shows  the 
behavior  of  njn^  and  r/r0  (o r/z/p0)  as  a  function  of  the 
FIR  electric  field  [  ( 8n7c )  ( /  / 2 )  ] 1 ' 2,  where  1  is  the  measured 
linearly  polarized  incident  intensity  so  that  /  /2  is  the  intensi¬ 
ty  of  the  circularly  polarized  component  that  causes  Landau 
level  transitions.  As  Fig.  4  shows  above  20  V/cm,  n ,  is 
roughly  constant  while  r  decreases  with  increasing  electric 
field. 

Since  the  effective  density  is  independent  of  intensity, 
our  data  do  not  appear  to  be  consistent  with  an  incoherent 
saturation  mechanism.  Our  results  do  not  agree  with  the 
observation  by  Helm  et  al?  of  partial  saturation  of  a  2DEG 
with  laser  intensities  of  100  mW/cm’  in  0.3  ms  pulses  at  a 
frequency  of  84. 2  cm'1.  Although  the  pulse  times  in  the  two 
experiments  differ  by  104,  a  comparison  can  be  made  because 


both  times  are  much  longer  than  the  relaxation  time  for  the 
2DEG. 

The  field  dependence  of  r  obtained  from  the  FIR  data 
appears  to  show  nearly  the  same  behavior  as  has  been  found 
in  the  earlier  dc  mobility  studies.  The  lines  superimposed  on 
top  of  our  results  in  Fig.  4  show  the  electric  field  dependence 
of  the  normalized  dc  mobility  measured  by  Shah  et  al.  on  two 
different  GaAs/GaAlAs  multiple  quantum  well  hetero¬ 
structures.  1,2  Because  of  the  different  way  in  which  the  two 
kinds  of  measurements  have  been  made  this  agreement  is 
quite  surprising.  Whereas,  at  resonance  in  the  FIR,  electron 
heating  could  proceed  in  the  same  way  as  has  been  proposed 
for  the  dc  case;  off  resonance,  this  same  heating  mechanism 
must  become  ineffective.  Inspection  of  curves  (a)  and  (f)  in 
Fig.  3  shows  that  off  resonance  the  absorption  is  greater  at 
high  FIR  intensities  than  it  is  at  low  ones,  an  observation 
which  is  not  consistent  with  a  resonant  process. 

We  conclude  that  electron  heating  does  occur  when  cy¬ 
clotron  resonance  is  probed  at  large  FIR  laser  intensities  and 
that  the  FIR  electric  field  dependence  of  the  relaxation  time 
is  the  same  as  that  previously  found  for  dc  mobility  measure¬ 
ments.  Since  the  heating  does  not  follow  the  resonance  shape 
of  the  conductivity,  the  good  agreement  between  these  two 
very  different  kinds  of  measurements  indicates  that  a  nonre¬ 
sonant  process  must  be  important. 
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Proton  Tunneling  with  Millielectrovolt  Energies  at  the  Be-H  Acceptor  Complex  in  Silicon 
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High-resolution  temperature-dependent  measurements  of  the  infrared  spectra  associated  with  the 
acceptor  complexes  Be-H  and  Be-D  in  Si  together  with  direct  far-infrared  studies  of  the  low-lying 
transitions  provide  an  unambiguous  identification  of  motional  tunneling  at  electronic  defects  in  ele¬ 
mental  semiconductors.  This  novel  system  consists  of  a  proton  or  deuteron  tunneling  around  the 
substitutional  Be  double  acceptor. 

PACS  numbers:  71 ,55.Fr 


A  dynamical  tunneling  model  has  been  proposed  to 
account  for  the  unusual  electronic  states  found  for 
shallow  acceptor  and  donor  complexes  containing  hy¬ 
drogen  in  ultrapure  germanium.1"4  Tunneling  of  the 
H  ion  around  the  substitutional  impurity  should,  in 
general,  alter  the  multiplicity  and  symmetry  of  the 
electronic  ground  state  so  that  its  behavior  is  com¬ 
pletely  different  from  that  found  for  a  monatomic  sub¬ 
stitutional  impurity.  Because  these  special  defects 
have  only  been  observed  with  the  sensitive  photother- 
mal  ionization  spectroscopic  technique,  a  complete  ex¬ 
perimental  determination  of  the  tunneling  properties 
has  not  been  possible;  hence  the  model  fits  to  the 
available  experimental  data  do  not  provide  a  unique 
explanation.4  If  dynamic  tunneling  does  play  an  im¬ 
portant  role  in  the  electronic  structure  of  these  de¬ 
fects,  then  this  internal  degree  of  freedom  may  also 
appear  in  the  electronic  spectra  of  other  complex  im¬ 
purity  centers  which  involve  the  lightest  element.  We 
have  searched  for  such  effects  in  two-component  com¬ 
plexes  in  Si. 

In  this  paper  the  first  observation  of  proton  tunnel¬ 
ing  at  deep-level  acceptor  complexes  in  silicon  is 
described.  High-resolution  temperature-dependent 
measurements  of  the  infrared  (ir)  spectra  associated 
with  the  acceptor  complexes  Be-H  and  Be-D  together 
with  direct  far-infrared  (fir)  studies  of  the  low-lying 
transitions  are  combined  with  analysis  based  on  the 
tunneling  model2  to  produce  an  unambiguous  picture 
of  a  proton  or  deuteron  tunneling  around  the  substitu¬ 
tional  Be  double  acceptor.  Large  isotope  shifts  in  the 
electronic  and  protonic  spectra  demonstrate  that  a 
strong  interaction  occurs  between  the  hole  in  its  elec¬ 
tronic  ground  state  and  the  proton  which  tunnels 
between  four  ( 11 1 )  symmetry  directions  near  the  ac¬ 
ceptor  site. 

Pure  fioat-zone-refined  silicon  wafers  have  been 
doped  with  use  of  diffusion  techniques.  Procedures 
for  double  doping  which  are  somewhat  different  from 
those  described  in  earlier  works5  6  are  outlined  here. 
Samples  for  the  ir  and  fir  have  been  treated  differently 
because  of  the  discrepancy  in  size  (2.5  x6x  17  mm3  vs 
5x10x25  mm3).  Typically,  for  the  ir,  a  Si  sample 
with  about  1  Mm  of  Be  evaporated  onto  two  sides  is 


sealed  off  in  a  carbonized  quartz  ampule  together  with 
a  SOO-Torr-He  +  200-Torr-H  gas  mixture  and  then 
heated  to  1300°C  for  20  min  and  air  quenched.  For 
the  fir,  the  standard  sandwich-sample  technique5  is 
employed,  but  the  quartz  ampule  is  charged  with  I  atm 
of  H2  and  maintained  at  1300*0  for  1  h  to  permit  dif¬ 
fusion  deep  into  the  sample.  The  optical  absorption 
spectra  are  measured  with  fast  Fourier-transform 
Michelson  and  Lamellar  grating  interferometers. 

Low-temperature  ir  spectra  for  three  acceptor  com¬ 
plexes,  /4(Li,Be),  y4(D,Be),  and  /l(H,Be),  are  shown 
in  Fig.  1.  The  acceptor  /l (Li, Be),  an  example  of  a 
nontunneling  center,5  produces  eight  spectral 
features7 ;  the  strongest  of  these  are  numbered  in  the 
figure.  With  increasing  temperature  a  second  se- 
que>  :e  of  lines  appears  shifted  to  lower  frequency  by 
1’..3  cm"1  with  respect  to  the  trace  shown  in  Fig.  1. 
These  data  demonstrate  that  the  T*  ground  state  is 
split  into  two  Kramers  doublets  (A E  —  11.3  cm-1)  as 
a  result  of  symmetry  lowering  ( T4  —  C3»)  produced 
by  the  nearby  interstitial  Li+  ion.5 

The  absorption  spectra  of  the  /4(D,Be)  and  ^(H,Be) 
centers  which  are  also  presented  in  Fig.  1  display  new 
lines  in  addition  to  the  main  sequence.  The  ratio  of 
the  line  strengths  of  the  secondary  to  the  main  transi¬ 
tions  is  observed  to  be  constant  for  many  samples 
doped  with  a  variety  of  diffusion  procedures.  A  simple 
sum  of  the  two  absorption  spectra  has  been  obtained 
for  a  sample  heated  in  a  mixture  of  H}  and  D2  gas, 
demonstrating  that  the  complex  contains  the  hydrogen 
isotope  in  the  monatomic  form.  Further  analysis  of 
the  spectra  shows  two  additional  distinguising  features: 

(1)  When  D— *  H  the  eight  transitions  of  the  main 
sequence  (some  of  the  strong  transitions,  labeled  I,  are 
recorded  in  Table  I)  all  display  a  large  isotope  shift  of 
7.8  cm-1  to  lower  frequencies.  This  isotope  effect  is 
an  order  of  magnitude  larger  than  that  previously  re¬ 
ported  for  shallow  dynamical  tunneling  centers  in  ger¬ 
manium.1 

(2)  The  low-temperature  fine  structure  is  a  weak 
high-frequency  replica  (labeled  T  in  Table  I)  of  the 
main  sequence  with  a  shift  of  16.2  cm"1  for  ^(D,Be) 
and  38.7  cm"1  for  /4(H,Be).  Note  that  the  complete 
replica  spectrum  can  be  identified  for  the  >4(D,Be) 
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Frequency  (cm*1) 

FIG.  1.  Infrared  transmittance  of  deep-level  acceptor 
complexes  in  Si.  The  spectra  measured  with  the  sample  at 
1.7  K  are  normalized  to  70-K  data  so  that  the  bulk-Si  pho¬ 
non  absorption  divides  out.  In  each  case  the  ordinate  ex¬ 
tends  from  0  to  1.  Instrumental  resolution  is  0.5  cm-1. 


center  in  Fig.  1,  but  only  one  replica  line  can  be  clearly 
observed  for  the  /4(H,Be)  center  because  of  accidental 
coincidences  with  the  strong  lines  of  the  main  se¬ 
quence.  The  apparent  simplicity  of  the  low-temper¬ 
ature  /<(H,Be)  spectrum  makes  it  clear  why  earlier 
workers,6  who  did  not  investigate  the  heavy  isotope, 
missed  this  effect. 

The  multiplicity  of  the  low-lying  energy-level 
scheme  can  be  seen  most  easily  with  temperature-de¬ 
pendent  measurements.  The  spectra  for  the  A(D, 
Be)  center  at  a  number  of  temperatures  are  shown  in 
Fig.  2.  As  the  temperature  is  increased  from  1.7  K 
12(a))  to  15  K  (2(d)!  two  additional  replicas  (II*  and 


Frequency  (cm'1) 

FIG.  2.  Temperature  dependence  of  the  ir  transmittance 
for  the  /<(D,Be)  center.  The  temperatures  are  (a)  1.7  K,  (b) 
5  K.  (c)  8  K.  and  (d)  15  K.  Again  the  bulk  phonon  contri¬ 
bution  is  removed  by  normalizing  to  the  70-K  spectrum. 
The  main  acceptor  absorption  series  at  low  temperature  is 
identified  in  (a)  by  the  markers  labeled  I.  while  the  replica 
of  this  series,  which  also  occurs,  is  labeled  1*.  With  increas¬ 
ing  temperature  other  absorption  lines  appear  which  follow 
the  series  pattern.  These  have  been  labeled  II*,  III*,  and  IV* 
in  succession.  The  absorption  series  II*  is  accompanied  by  a 
lower-frequency  replica  II.  The  instrumental  resolution  is 
0.5  cm"1. 

Ill*)  appear  on  the  high-frequency  side  of  the  main 
lines.  This  blue  shift  is  quite  unusual  and  indicates 
that  splittings  in  the  electronic  excited  states  must  oc- 


TABLE  I.  ir  acceptor  plus  tunneling  transition  frequencies  (cm"1)  for  strong  lines  1,  2, 
and  4  of  /HD, Be)  and  ,4(H,Be).  Since  the  transition  labeled  4  is  a  doublet  (A£  -  1.9  ±0.1 
cm-1),  the  recorded  values  are  for  the  lower  of  the  two  lines  unless  they  are  not  resolved 
in  which  case  the  center  of  gravity  is  given. 


Electronic 

level 

i 

I* 

II 

II* 

III* 

1  (r,) 

633.3 

—  649 

D 

—  642 

2  (T.) 

663.9 

680.0 

658.0 

673.9 

670.3 

4  (r«  +  r;) 

704.3 

720.5 

698.3 

714.4 

710.6 

i  (r,) 

625  4 

Hidden 

H 

Hidden 

2  (r,) 

655.7 

Hidden 

642.1 

680.6 

f>b8  3 

4  (r*+r7) 

696  5 

736.3* 

Hidden 

721.5* 

- ;  0  5* 

•Not  resolved. 
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FIG.  3.  Far-infrared  absorption  spectrum  of  ^(H,Be)  and 
/4(D.Be).  The  absorption  coefficient  difference  f/t(H,Be) 
-  4(D,Be)l  is  shown  vs  frequency  for  a  sample  temperature 
of  1.2  K.  The  instrumental  resolution  is  0.8  cm-1  and  the 
/4(D,Be)  lines  are  unresolved. 

cur.  Absorption  series  II*  and  II  [shown  in  Fig.  2 (b)  1 
have  the  same  thermal  activation  energy  and  hence 
start  from  the  same  excited  state.  Inspection  of  Table 

I  indicates  that  the  energy  difference  between  II*  and 

II  is  the  same  as  that  between  I*  and  I.  At  still  higher 
temperatures  the  series  III*  appears  without  its  low- 
frequency  partner  and  series  IV*,  which  overlaps  II, 
also  may  exist. 

Ir  spectra  for  /l(H,Be)  show  quite  similar  tempera¬ 
ture  dependences,  but  the  spectral  changes  occur  at 
somewhat  higher  temperatures.  The  corresponding 
transitions  are  recorded  in  Table  I. 

The  fir  absorption  spectrum  produced  by  direct  exci¬ 
tation  of  the  tunneling  states  is  shown  in  Fig.  3  where 
the  difference  in  absorption  coefficient  [/4(H,Be) 
— /l(D,Be)l  is  plotted  versus  frequency.  All  of  the 
transition  frequencies  which  are  given  in  Table  II  agree 
with  those  obtained  from  the  ir  measurements.  Note 
that  both  ratios  of  H/D  for  the  transition  frequencies 
A  or  B  give  —  2.5  which  is  too  large  to  stem  from  vi¬ 
brational  or  rotational  motion  of  the  defect  but  is  con¬ 
sistent  with  dynamical  tunneling.  A  new  observation 
is  the  rapid  dependence  of  the  tunneling  linewidth 
with  frequency.  Since  both  spectral  widths  for  the  hy¬ 
drogen  isotope  are  resolved  in  this  figure,  one  can 
determine  that  the  inverse  lifetime  varies  as  w2,  and 
hence  is  proportional  to  the  Si  phonon  density  of 
states.  Higher-resolution  measurements  on  the  deu¬ 
terium  isotope  indicate  that  the  9.8-cm"1  mode 
linewidth  is  (FWHM)  —0.3  cm-1  while  the  5.9-cm"1 
mode  has  a  width  <  0.2  cm"1.  In  a  magnetic  field  of 
6.7  T  all  of  these  lines  show  a  weak  magnetic  field 
dependence,  i.e.,  the  field-induced  frequency  shift  is 
only  resolvable  for  the  lowest-frequency  deuterium 
line. 

What  kind  of  center  can  produce  both  the  giant  iso- 


TABLE  II.  Far-infrared  tunneling  transition  frequencies 
A,  B,  and  C  in  inverse  centimeters  for  /t(D,Be)  and 
/4(H,Be). 


fir 

D 

H 

D 

ir 

H 

c 

12.9  (weak) 

B 

9.75 

25.5 

9.8 

26 

A 

5.9 

14.5 

6.0 

-14 

tope  effects  and  the  complex  electronic  spectra?  Be¬ 
cause  Be  in  Si  acts  as  an  acceptor  and  belongs  to  group 
II,  it  is  considered  to  be  a  substitutional  double  accep¬ 
tor.5  Since  Li  is  an  interstitial  donor  then  as  a  result  of 
compensation  the  complex  ^(Li,Be)  can  be  pictured5 
as  the  single  acceptor  (Be2-  +Li+  +h+).  Although 
hydrogen  usually  acts  as  an  acceptor  in  elemental  sem¬ 
iconductors,8  it  is  natural  when  combining  H  with  sub¬ 
stitutional  Be  to  view  it  as  a  donor,  just  like  the  Li 
case.  The  resultant  defect  configuration  is  (Be2- 
+  p+  +  h+ ),  which  should  be  contrasted  with  the  pre¬ 
viously  reported2  tunneling  acceptors  (Si  +  H"  +h+) 
and  (C  +  H"+h+)  in  Ge.  Of  course,  the  general 
theoretical  formalism  for  a  dynamic  tunneling  system 
developed  by  Haller,  Joos,  and  Falicov2  is  still  valid 
and  we  now  outline  its  application  to  the  system  of  in¬ 
terest. 

Our  data  together  with  the  theoretical  model  of  an 
acceptor  with  a  tunneling  proton  produces  energy  level 
schemes  such  as  the  one  shown  in  Fig.  4.  This  ls-like 
electronic  ground  state  of  a  substitutional  acceptor  is  a 
fourfold-degenerate  Tg  state  while  a  p-like  excited 
state  has  r«,  T7,  or  r8  symmetry.7  A  noninteracting 
proton  tunneling  between  the  four  interstitial  sites 
around  the  Be  gives  rise  to  a  singlet  T ]  and  a  triplet  r5 
level.  Shown  on  the  left-hand  side  of  Fig.  4  is  the  en¬ 
ergy  scheme  with  the  proton-hole  interaction  turned 
off  so  that  each  electronic  state  (0,1, 2,4)  has  available 
two  proton-tunneling  configurations.  This  weak- 
coupling  limit  is  expected  to  describe  the  p-Iike  excited 
states  since  the  hole  wave  function,  which  vanishes  at 
the  defect  origin,  has  a  small  amplitude  in  the  proton 
neighborhood;  however,  for  this  deep  center  a  strong 
interaction  between  the  1  s-like  electronic  ground  state 
(Tg)  and  the  tunneling  proton  should  occur.  The 
resultant  level  structure  is  pictured  on  the  right-hand 
side  of  Fig.  4.  According  to  Haller,  Joos,  and  Falicov2 
the  three  T t  levels  have  mixed  proton-tunneling  confi¬ 
gurations  while  the  r6  and  r7  states  have  only  Ts  tun¬ 
neling  character;  hence  two  sets  of  optical  transitions 
are  possible  from  a  Tg  state  but  only  one  from  r4  or 

r7. 

The  assignments  for  the  observed  absorption  series, 
designated  by  I,  I*.  II,  II*,  and  III*,  are  shown  in  Fig. 
4.  The  asterisk  identifies  the  ir  transition  to  the  excit¬ 
ed  state  with  the  higher-energy  proton-tunneling  con- 
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FIG.  4.  Energy-level  diagram  of  the  ir  and  fir  transitions 
observed  for  YH.Be)  and  YD, Be).  The  numbers  on  the 
left-hand  side  identify  a  few  of  the  strong  acceptor  transi¬ 
tions.  Next  is  shown  the  multiple!  arrangement  for  the  ac¬ 
ceptor  GY  r,f  r,,  and  rt  +  r7)  and  tunneling  proton 
GY  Tj)  with  no  interaction.  A  strong  coupling  between  the 
hole  and  motional  degrees  of  freedom  for  the  ground-state 
manifold  is  shown  at  the  lower  right.  The  vertical  transi¬ 
tions  are  described  in  the  text. 


figuration.  The  deuteron  or  proton  tunnel  splitting, 
called  4r  in  Ref.  2,  is  found  directly  from  the  energy 
difference  between  the  series  IM  or  II*-II  to  be  16.2 
cm-1  for  /KD.Be)  and  38.8  cm-1  for  ,4(H,Be).  The 
ground-state  splittings  calculated  from  the  difference 
between  1  and  II,  I*  and  II*,  etc.,  are  both  consistent 
with  the  temperature-dependent  profile  and  the  fir 
measurements  as  shown  in  Table  II. 

In  conclusion,  three  new  distinguishing  features 
have  been  found  for  this  novel  dynamical  tunneling 
system:  (i)  The  observation  of  the  low-temperature  ir 
replica  spectrum  makes  possible  the  first  direct  deter¬ 


mination  of  the  tunneling  parameter,  (ii)  The  appear¬ 
ance  of  other  ir  replicas  at  elevated  temperatures  and 
the  direct  observation  of  fir  transitions  provide  confir¬ 
mation  of  the  tunnel  splittings  in  the  ground-state 
manifold,  (iii)  The  large  ir  and  fir  isotope  shifts 
demonstrate  that  a  strong  interaction  exists  between 
the  electronic  ground  state  of  the  hole  and  the  tunnel¬ 
ing  proton  for  this  deep-level  defect.  Previous  pho- 
tothermal  ionization  spectroscopic  studies1"4  of  shal¬ 
low  dynamical  tunneling  centers  in  Ge  have  been 
made  with  very  low  defect  densities  (  —  1011  cm-2)  so 
that  the  small  isotope  effect  (  <0.5  cm-1)  can  be 
resolved  but  no  direct  determination  of  the  tunneling 
energy  for  these  delicate  centers  has  been  possible. 
Our  discovery  that  large  isotope  effects  can  be  ob¬ 
served  for  a  deep-level  tunneling  center  in  Si  (density 
~  1016  cm-3)  opens  up  the  experimental  possibilities. 
Now  a  variety  of  spectroscopic  tools  can  be  employed 
to  probe  more  completely  the  dynamical  behavior  of 
these  unusual  yet  robust  electronic  defects. 
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Our  experimental  attenuated-total-internal-reflection  curves  for  Ag  films  are  reanalyzed  to  in¬ 
clude  the  thin  contamination  layer  which  has  been  proposed  by  W.  H.  Weber  to  account  for  the 
observed  differences  in  the  optical  constants  at  the  metal-air  and  metal-liquid  interfaces.  In  con¬ 
tradiction  to  Weber's  claim,  the  Orude  parameters  of  the  Ag  film  still  show  the  same  behavior  as 
when  this  adlayer  is  ignored. 


It  has  been  proposed1  that  the  reported  differences  in 
optical  constants  of  Ag  films  measured  at  metal-air  and 
metal-liquid  interfaces2  can  be  attributed  to  the  thin  con¬ 
tamination  layer  that  exists  on  Ag  films  which  have  been 
exposed  to  air.  In  this  Comment,  we  describe  our  test  of 
Weber’s  hypothesis  and  show  that  the  inclusion  of  such  an 
adlayer  into  the  data  analysis  does  not  alter  our  original 
interpretation  of  the  experimental  results. 

The  values  of  the  complex  dielectric  function  of  Ag  are 
calculated  by  a  nonlinear  least-squares  fit  of  the  exact 
Fresnel  reflection  formulas  to  the  experimental  attenuat¬ 
ed-total-reflection  curves.  The  three-layer  model  used  in 
the  original  analysis2  is  now  replaced  by  a  fojir-layer 
model.  The  added  layer  has  a  thickness  of  12  A  and  a 
dielectric  constant  «0— 4,  values  chosen  by  Weber  to  mini¬ 
mize  the  apparent  difference  in  the  real  part  of  the  Ag 
dielectric  function  for  the  metal-air  and  metal-liquid  inter¬ 
face.  The  calculated  values  of  «2,  and  r,-1  for  this 
four-layer  model  are  given  in  Table  I.  The  values  obtained 
with  the  three-layer  model  are  given  in  Ref.  2. 

The  change  in  the  Drude  parameters  are  most  clearly 
seen  with  plots  of  «,  versus  wavelength  squared,  and  rf-1 
versus  frequency  squared  for  each  model.  Since 


and 

*e~'’"To'+Po0)1  ,  (2) 

the  slope  of  the  straight-line  fits  to  the  data  in  Fig.  1  is 
directly  related  to  m2pt  and  the  slope  in  Fig.  2  to  p0. 

Figure  l  shows  our  original  «,  values  for  a  three-layer 
least-squares  fit  to  A2,  and  the  new  calculated  values  of  a 
four-layer  system  with  the  adlayer  parameters  assumed  by 
Weber.  Inspection  of  this  figure  demonstrates  that  the  ad¬ 
dition  of  the  layer  to  the  Ag  surface  does,  indeed,  alter  the 
relative  intercept  and  change  the  relative  values  of  the  two 
slopes  for  air  and  for  liquid.  The  result  overall  is  approxi¬ 
mately  the  same  as  Weber’s  correction  equation  for  «, 
would  imply.1  Note  that  for  the  four-layer  system,  al¬ 
though  the  two  curves  cross  at  some  wavelength,  the  slopes 
are  different  with  copt  still  larger  in  the  metal-liquid  inter¬ 
face  than  for  the  air-metal  interface  (9.56  ±  0.03  eV  com¬ 
pared  to  9.43  ±  0.03  eV).  At  this  point  it  may  appear  that 
the  four-layer  model,  with  perhaps  slightly  different  pa¬ 
rameters,  could  make  the  proposed  indcx-of-refraction- 
induced  effect  vanish.  We  now  examine  the  imaginary 
part  of  the  metal  dielectric  function  and  demonstrate  that 
this  is  not  the  case. 

By  means  of  Eq.  (6)  in  Ref.  2,  the  inverse  relaxation 
time  is  obtained  from  the  values  of  c,  and  t2  for  the  three- 
and  four-layer  model.  In  Fig.  2  the  triangles  identify  the 


TABLE  I.  Experimentally  determined  dielectric  function  of  Ag  and  related  optical  parameters  for 
the  four-layer  models.  The  adlayer  has  a  thickness  of  12  A  and  a  dielectric  constant  of  4. 


CCI, 
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A  (A) 

«2 

r“'  (I014  sec'1) 

«i 

*2 

r"'  (I014  sec-1 

4416 

-6.881 

0.413 

1.538 

-6.971 

0.565 

2.171 

-9.639 

0.533 

1.446 

-9.626 

0.630 

1.770 

5145 

-11.148 

0.653 

1.520 

-11.238 

0.721 

1.719 

5309 

-12.153 

0.692 

1.467 

-12.191 

0.759 

1.650 

5682 

-14.613 

0.843 

1.456 

-14.587 

0.830 

1.471 

6328 

-18.925 
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1.375 

—  19.032 

1  148 

1.477 

6571 

-20.385 

1.136 

1.324 

-20.130 

1.127 

1.353 

6764 
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-22.107 
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FIG.  1.  Measured  real  part  of  the  dielectric  function  of  silver 
vs  wavelength  squared.  For  air  (A,  dashed  line)  and  CC14  (a, 
solid  line)  half  spaces  using  the  three-layer  model  described  in 
Ref.  2.  For  air,  (O,  dashed)  and  CC14  ( •,  solid)  with  the  four- 
layer  model  described  in  the  text. 


three-layer  case  and  the  circles  the  four  layer.  For  each 
model,  the  change  is  identified  by  the  appropriately  la¬ 
beled  arrow.  In  contrast  with  the  result  found  in  Fig.  1, 
now  the  presence  of  the  same  adlayer  increases  the  differ¬ 
ence  between  the  two  slopes  over  that  given  by  the  three- 
layer  case.  For  both  models,  f)0  is  smaller  for  the  metal- 
liquid  interface. 

These  comparisons  in  the  same  figures  of  the  four-  and 
three-layer  results  demonstrate  that  both  models  produce 
the  same  qualitative  behavior  for  the  Ag-film  Drude 
parameters;  namely,  ^'(liquid)  >  r0-1  (air),  /J(liquid) 
</J(air),  and  cop (liquid)  >  affair).  Also  note  that  when 
the  adlayer  parameters  are  chosen  in  such  a  way  as  to  min¬ 
imize  the  difference  between  ^(liquid)  and  cop( air),  the 
difference  between  ^(liquid)  and  /J( air)  is  actually  in¬ 
creased.  The  adlayer  does  not  produce  the  required  effect 
of  decreasing  the  difference  for  both  Drude  parameters  at 
the  same  time. 

Our  complete  data  analysis  using  both  the  three-  and 
four-layer  models  shows  that  the  differences  in  the  Ag- 
film  optical  constants  which  occur  when  the  air  half  space 
is  replaced  by  the  liquid  one,  are  not  eliminated  when 
Weber's  hypothetical  thick  adlayer  is  included.  It  should 
be  mentioned  at  this  point  that  electron-spectroscopy  for 
chemical  analysis  measurements  on  our  Ag  films  showed 


FIG.  2.  Measured  relaxation  frequency  vs  frequency  squared. 
For  air  (A,  dashed  line)  and  CC14  (▲,  solid  line)  half  spaces  us¬ 
ing  the  three-layer  model.  For  air  (O,  dashed)  and  CCI4  (•, 
solid)  with  the  four-layer  model. 


about  3  A  of  physisorbed  oxygen,  less  than  a  monolayer  of 
chemisorbed  oxygen,  and  no  sulfur.  Whether  such  a  layer 
should  be  added  to  the  three-layer  model  without  at  the 
same  time  treating  surface  roughness  is  a  moot  point  since 
both  components  are  clearly  present  on  the  Ag-film  sur¬ 
face  and  they  produce  corrections  of  opposite  sign  for 
Re(«Ag)  w^en  >ndex  of  the  neighboring  half  space  is 
changed. 

Weber  has  argued  qualitatively  that  an  absorbing  ad¬ 
layer  can  account  for  the  systematic  changes  which  we 
have  found  for  Im(«Ag)  between  the  air  and  liquid  inter¬ 
face.  The  argument  seems  reasonable  until  one  acknowl¬ 
edges  the  restrictions  imposed  on  the  optical  properties  of 
this  adlayer.  First,  at  frequencies  below  the  interband 
transitions  the  optical  constants  of  noble  metals  obey  the 
frequency  dependencies  described  by  Eqs.  (1)  and  (2) 
both  in  vacuum  and  in  air.2  Second,  irrespective  of  the 
mechanism,  the  parameters  must  obey  Kramers-Kronig 
(KK)  relations.  In  this  Comment  we  have  shown  that  the 
transparent  adlayer,  which  automatically  satisfies  KK  re¬ 
lations,  preserves  the  known  frequency-dependent  proper¬ 
ties  of  Ag,  i.e.,  a  plasma  frequency  which  is  constant  and  a 
relaxation  rate  which  varies  quadratically  with  frequency 
(see  Figs.  I  and  2).  Forcing  the  adlayer  to  absorb  strongly 
in  the  visible  region  causes  its  index  of  refraction  to  be  fre¬ 
quency  dependent  as  well.  To  account  for  our  experimen¬ 
tal  results  on  Ag,  these  frequency  dependences  must  con¬ 
spire  in  such  a  way  to  give  the  same  frequency-dependent 
properties,  as  found  for  the  nonabsorbing  adlayer  de¬ 
scribed  above.  This  line  of  reasoning  leads  us  away  from 
the  absorbing  adlayer  concept. 

The  recent  attenu^ted-total-reflection  measurements  on 
silver  films  at  5145  A  by  Kurosawa,  Pierce,  Ushioda,  and 
Hemminger’  provide  additional  experimental  information 
on  the  index-of-refraction-induced  change.  From  their 
measurements  they  conclude  that  the  polar  liquid-induced 
changes  in  the  fig  optical  constants  can  be  explained  by 
assuming  a  10-A-ihick  nonabsorbing  adlayer  on  the  silver 
surface  with  dielectric  constant  3;  however,  the  changes  in 
optical  constants  were  not  measured  as  a  function  of  fre- 
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quency.  In  this  Comment  we  have  shown  that  frequency- 
dependent  measurements  are  required  to  demonstrate  that 
an  adlayer  alone  cannot  account  for  the  index-of- 
refraction-induced  changes  in  both  optical  constants. 

In  Weber’s  conclusions,  it  is  argued  that  Table  I  in  the 
earlier  work  by  Eagen  and  Weber4  shows  that  an  index- 
of-refraction-induced  change  in  the  optical  constants  does 
not  occur.  This  table  presents  values  of  «,  and  t  with  es¬ 
timated  uncertainties  at  a  number  of  wavelengths  for  one 
Ag  film  measured  at  a  metal-vacuum  interface  and  for 
another  film  measured  at  a  metal-substrate  interface.  Our 
least-squares  fits  to  these  data  give  Drude  plasma  frequen¬ 
cies  which  differ  by  more  than  the  uncertainties  of  the 
measurement.  These  Drude  parameters  are  listed  in  Table 
II.  Although  these  differences  are  comparable  in  magni¬ 
tude,  but  opposite  in  sign,  to  the  differences  reported  in 
our  measurements  on  a  single  dirty  film,'  the  variations  in 
the  optical  properties  from  sample  to  sample  are  large 
enough  so  that  the  differences  evident  in  Table  II  are  not 
deemed  significant1,4  by  the  authors.  We  conclude  that 
because  of  these  sample-to-sample  limitations,  a  compar¬ 
ison  of  the  optical  properties  of  the  two  films  listed  in 
Table  I  of  Ref.  4  cannot  be  used  to  argue  for  or  against  the 
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TABLE  II.  Drude  parameters  for  the  two  Ag  films  described 
in  Ref.  4. 

Metal-vacuum  Metal-substrate 

o)  (eV)  9.77  ±0.03  9.39  ±0.03 

P0  (1012  sec~leV J)  8.2 ±  1.4  3.40±0.5 


index-of-refraction-induced  changes  which  we  have 
described. 

In  conclusion,  our  analysis  demonstrates  that  the  pres¬ 
ence  of  a  thin  contamination  layer  on  Ag  films  does  not 
qualitatively  alter  our  original  interpretation  of  the  experi¬ 
mental  results.  An  index-of-refraction-induced  change 
still  is  observed  in  the  Drude  parameters  of  Ag  films. 
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An  IR  grating  on  a  clean  W(100)  surface  is  shown  to  generate  both  homogeneous  and  inhomogeneous  surface 
electromagnetic  waves.  An  observed  interference  between  these  two  components,  which  can  be  described  in  terms 
of  a  two-beam  interferometer  with  variable  arm  amplitude  and  fixed  optical  path,  is  used  to  measure  the  plasma 
frequency  accurately  in  the  IR. 


The  attenuation  of  surface  electromagnetic  waves 
(SEW’s)  propagating  on  a  clean  W(100)  surface  has 
been  used  to  monitor  adsorbate  vibrational  modes,1 
surface  reconstruction,2  and  chemisorption-induced3 
changes  in  the  free-carrier  behavior  of  the  metal  in  the 
room-temperature  regime.  We  report  on  the  mea¬ 
sured  properties  of  SEW’s  at  elevated  temperatures. 
At  temperatures  greater  than  1000  K  the  SEW  signal 
is  attenuated  to  such  a  large  extent  that  low-intensity 
surface  skimming  plane  electromagnetic  waves 
(PEW’s),  which  are  also  generated  at  the  grating  cou¬ 
pler,  can  be  detected.  Since  both  kinds  of  wave  are 
generated  coherently  at  the  input  but  travel  across  the 
surface  with  different  velocities,  interferometry  is  pos¬ 
sible.  This  interferometer  has  been  used  to  measure 
the  plasma  frequency  of  W  in  the  10-/im  wavelength 
region. 

A  description  of  the  experimental  apparatus  and 
procedures  for  making  SEW  attenuation  measure¬ 
ments  in  ultrahigh-vacuum  UHV  conditions  has  been 
given  in  Ref.  2.  Because  the  SEW  wave  vector  is 
greater  than  that  of  light,  gratings  etched  into  the 
surface  are  used  to  couple  CO2  laser  radiation  into  and 
out  of  the  SEW  spectrum. 

To  identify  the  interference  signature  between 
SEW’s  and  PEW’s,  the  temperature  dependence  of 
the  signal  from  the  output  coupler  is  measured  at 
many  laser  frequencies  across  the  ranges  of  the  12C02 
and  13C02  laser  gases.  Each  of  the  eight  data  traces 
shown  in  Fig.  1  represents  an  intensity-versus-tem- 
perature  run  at  a  fixed  frequency.  Note  the  strong 
minimum  in  log  (intensity)  that  occurs  near  T  =  400°C 
and  v  =  1000  cm"1.  For  all  frequencies  the  tempera¬ 
ture  dependence  at  low  temperatures  agrees  with  that 
calculated  for  SEW’s  from  the  temperature  depen¬ 
dence  of  the  dc  resistivity  p(T),  namely, 

ls(T)  =  IoS  exp(-o)2Lp(T)/4Tc],  (1) 

where  w  is  the  angular  frequency,  L  the  propagation 
distance,  and  c  the  velocity  of  light.  For  high  tem¬ 
peratures  an  additional  component,  which  tracks  the 
SEW  intensity  when  the  laser’s  beam  angle,  input 
power,  or  polarization  is  changed,  comes  from  PEW’s 
generated  at  the  input  grating. 

Next  it  is  shown  that  the  phase  difference  between 
SEW’s  and  PEW’s  in  the  IR  is  large  enough  over  the 


sample  length  to  account  for  the  interference  effect. 
The  difference  in  phase  between  the  two  beams  can  be 
written  as 

6  =  qsL  -  qpL  +  <t> ,  (2) 

where  qs  is  the  SEW  wave  vector,  qp  the  PEW  wave 
vector,  and  <p  the  phase  shift  between  the  two  waves, 
which  is  assumed  to  be  a  constant  over  the  frequency 
region  of  interest.  For  a  Drude  metal  in  the  10-pm- 
wavelength  region 

qs  ~  («/c)[l  +  w2/2wp2  +  (c0  +  3/4)(l  -  l/u>2r2)w4/2u>p4], 

(3) 

where  up  is  the  plasma  frequency  in  the  IR  and  «0  the 
low-frequency  contribution  to  the  dielectric  constant 
from  the  interband  transitions.  Since  the  most  effec¬ 
tive  PEW’s  are  those  within  a  wavelength  of  the  sur¬ 
face  at  the  output  coupler,  qp  =  w/c;  hence  the  phase 
difference  between  the  two  arms  for  a  Drude  metal  is 

6  *  (a>L/c)(cir/2«p2)  +  4>-  (4) 

The  resultant  intensity  at  the  detector  now  takes  on 
a  familiar  form,  namely, 

/  =  Is  +  IP  +  2 (IsIP)in  cos  d.  (5) 

Destructive  interference  occurs  when  6  -  (2 n  +  l)x, 
with  n  a  nonnegative  integer,  and  the  largest  effect 
appears  when  the  two  arms  of  the  optical  bridge  are 
balanced.  A  temperature  sweep  at  a  fixed  laser  fre¬ 
quency  (the  data  shown  in  Fig.  1)  corresponds  to  vary¬ 
ing  the  relative  amplitudes  for  the  two  components  of 
Eq.  (5)  while  keeping  6  fixed. 

To  obtain  accurate  values  for  8(w)  each  data  set 
must  be  fitted  to  Eq.  (5),  so  knowledge  of  the  expected 
temperature  and  frequency  dependences  of  Is  and  Ip 
is  required.  The  SEW  intensity  dependence  is  given 
in  Eq.  (1).  Although  an  accurate  expression  for  the 
PEW  intensity  dependence  requires  the  solution  of  a 
deep  grating  diffraction  problem,  we  show  here  that 
the  appropriate  temperature-dependent  form  can  be 
constructed  from  a  simple  phenomenological  model. 

We  assume  that  the  near  fields  at  the  source  that 
evolve  into  interfering  PEW’9  and  SEW’s  extend  a 
distance  &n  above  the  grating  and  that  the  mean  SEW 
amplitude  height  is  a  measure  of  this  range.  Now 
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Temperature  (°C) 

Fig.  1.  Transmitted  IR  intensity  as  a  function  of  tempera¬ 
ture  and  frequency.  Eight  SEW  signal-versus-temperature 
scans  are  shown  at  approximately  equal  frequency  intervals 
across  the  CO2  laser  spectrum. 

consider  a  PEW  traveling  along  the  sample  surface. 
In  analogy  with  a  parallel-plate  transmission  line  ge¬ 
ometry  for  fixed  plate  separation  5q,  the  transmitted 
intensity  would  be 

IP  =*  I  op  exp(-2ry/50),  (6) 

where  r  is  the  normalized  real  part  of  the  surface  im¬ 
pedance  of  the  metal2  and  y  the  distance  along  the 
surface. 

Because  of  diffraction  the  height  of  the  PEW  actual¬ 
ly  changes  with  distance  from  the  source;  hence  from 
the  geometry 

Hy)  =  «o  +  y  tan(a),  (7) 

where  i0  is  the  height  at  y  *  0  and  a  =  sin_1(X/5o),  the 
angle  to  the  first  minimum  in  a  single-slit  diffraction 
pattern.  The  appropriate  generalization  of  Eq.  (6)  is 

IP  -  IQP  expj— I  dy[2r/5(y)]j  •  (8) 

Another  loss  mechanism  for  this  component  occurs  at 
the  output  coupler,  y  =  L,  since  only  the  fraction  [5o/ 
5(L)|  of  the  PEW  will  be  coupled  out  to  the  detector. 
The  resultant  expression  is 

IP  =  ;0P(1  +  /L/50)"u+2r/n,  (9) 

where 

f  =  A/(502  —  X2)l/2.  (10) 

Three  adjustable  parameters  8,  /os,  and  Iqp  in  Eqs. 
(1),  (5),  and  (9)  are  used  to  fit  the  temperature-depen¬ 
dent  data  with  the  data  near  Is  *  Ip  weighted  most 
heavily.  Figure  2  shows  three  data  sets  with  the  corre¬ 


sponding  best  fit  given  by  this  model.  The  destruc¬ 
tive  interference  between  the  SEW’s  and  the  PEW’s 
increases  in  strength  as  the  laser  frequency  is  de¬ 
creased  from  1072  cm-1  [Fig.  2(a)]  to  995  cm-1  [Fig. 
2(c)  ] .  Sample  movement  produces  the  high-tempera- 
ture  difference  between  the  data  and  the  model  shown 
in  Figs.  3(b)  and  3(c).  The  fits  show  that  the  PEW 
intensity  at  room  temperature  is  ~1%  of  the  SEW 
value.  Even  here  the  PEW  component  cannot  be  ne¬ 
glected  since  at  frequencies  where  a  r  phase  shift  oc¬ 
curs  the  resultant  intensity  is  ~20%  smaller  than  for 
the  SEW  component  alone. 

In  order  to  compare  the  results  with  relation  (3)  the 
measured  phase  angle  8  is  plotted  against  the  frequen¬ 
cy  cubed  as  shown  in  Fig.  3.  A  linear  dependence  is 
observed  over  the  measured  range  covering  nearly  x 
rad.  A  least-squares  fit  to  these  data  gives  L/(2cwp2) 
=  1.04(±0.04)  X  1025  sec"3  and  <t>  =  -0.51(±0.05)t. 

Our  measurements  show  that  the  output  coupling  of 
both  the  SEW’s  and  the  PEW’s  takes  place  over  the 
full  extent  of  the  3.8-mm-wide  output  grating,  so  half 
of  this  distance  is  added  to  the  separation  between 
gratings  to  obtain  the  total  propagation  distance. 


Fig.  2.  Comparison  of  three  sets  of  data  with  the  two-beam 
interferometer  model.  The  open  circles  represent  the  data 
and  the  solid  line  the  model  fit.  The  three  frequencies  are 
(a)  1072  cm"1,  (b)  1027  cm'1,  and  (c)  995  cm"1. 
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Pig.  3.  Frequency  dependence  of  the  phase  difference  be¬ 
tween  the  two  arms  of  the  interferometer.  The  filled  circles 
represent  the  data;  the  open  circles  are  alternative  values  for 
the  data  points  directly  above  them  [since  cos(t  +  x )  =*  cos(r 
—  x)|.  The  solid  line  shows  a  linear  least-squares  fit. 


Substitution  of  L  *  4.94  cm  gives  hup  =*  7.0(±0.3)  eV 
for  the  plasma  frequency  in  the  IR. 

By  measuring  the  transmitted  intensity  as  a  func¬ 
tion  of  placement  of  the  input  beam  on  the  input 
grating,  it  has  been  possible  to  show  that  the  PEW's 
occur  not  because  of  the  impedance  change  at  the 
grating-smooth-metal  boundary4  but  instead  because 
of  the  finite  number  of  grating  lines  intercepted  by  the 
input  beam  (about  20  in  the  actual  experiment).  For  a 
plane  wave  incident  upon  an  infinite  grating  at  the 
angle  for  maximum  SEW  generation  the  first-order 
PEW  beam  is  forbidden,  but  when  the  grating  consists 
of  a  small  number  of  lines  this  selection  rule  against 
the  first-order  beam  is  weakened,  so  some  PEW’s  do 
appear.  The  end  result  is  that  a  narrow  grating  cou¬ 
pler  or  a  small  spot  size  automatically  sends  PEW’s 
along  with  the  desired  SEW’s  across  the  surface. 

Another  feature  of  our  data  that  needs  examination 
is  the  finite  intensity  at  the  destructive  interference 
condition.  At  a  frequency  near  1000  cm-1,  where  8  ■ 
x,  the  resultant  intensity  should  dip  down  to  the  back¬ 
ground-noise  level  when  Is  m  Ip;  however,  the  smallest 
measured  value  is  ~4%  of  the  SEW  (or  PEW)  value  at 
that  temperature.  [This  is  for  the  data  at  995  cm"1 
shown  in  2(c).]  There  are  two  contributions  to  this  4% 
minimum.  The  first  is  that  the  SEW  beam  is  attenu¬ 
ated  more  strongly  than  the  PEW  beam  along  the  deep 
grating  couplers.  This  leads  to  less  than  complete 
destructive  interference  since  the  two  components’  in¬ 
tensities  cannot  be  matched  across  the  entire  length  of 
the  output  coupler  (and  hence  across  all.  the  detector). 
The  second  contribution  comes  from  the  difference  in 
SEW  and  PEW  wave  vectors.  This  difference  in  wave 
vector  (which  causes  the  whole  interference  effect  in 
the  first  place)  also  leads  to  less  than  complete  de¬ 
structive  interference  since  the  phase  angle  8  is  not  a 
constant  across  the  whole  output  coupler.  From  mea¬ 
sured  values  of  the  difference  in  attenuation  and  wave 
vector  of  the  two  components  we  calculate  minimum  1/ 


Is  (or  I/Ip)  »  2.5%,  in  reasonable  agreement  with  the 
measured  value  of  4%. 

The  experimentally  determined  value  of  the  con¬ 
stant  phase  factor  0  ■  -*72  cannot  be  explained  with¬ 
in  the  framework  of  this  simple  phenomenological  pic¬ 
ture.  A  rigorous  diffraction  calculation  may  be  re¬ 
quired  for  the  relative  phases  of  the  two  components 
at  the  input  and  output  coupler  to  be  explored. 

Our  direct  determination  of  the  real  part  of  the 
inverse  dielectric  function  of  a  clean  W(100)  surface 
that  we  have  characterized  with  an  effective  plasma 
frequency  (hup  *  7.0  eV)  in  the  Drude  model  approxi¬ 
mation  can  be  compared  with  the  value  deduced  for 
the  same  frequency  region  from  reflectivity5  and  ellip- 
sometric6  measurements  (hwp  =  6.4  eV).7  The  errors 
in  these  earlier  measurements  are  large  enough  that 
the  two  values  are  within  the  uncertainties. 

This  new  measurement  technique  is  not  limited  to 
W.  From  relation  (4)  the  criterion  for  destructive 
interference  (assuming  that  0  =  — *■/ 2)  can  be  rewrit¬ 
ten  in  terms  of  the  SEW  attenuation  coefficient  in  Fq. 
(1)  to  give 

oi2Lp(T)/4icc  —  3t/o)t.  (11) 

Inspection  of  the  right-hand  9ide  of  this  expression 
shows  that  the  IR  frequency  should  be  chosen  to  en¬ 
sure  that  ur  >  1  for  a  particular  metal  so  that  the 
attenuation  of  the  SEW  will  not  obscure  the  interfer¬ 
ence  effect 

The  analysis  given  here,  which  uses  the  interference 
between  SEW’s  and  PEW’s  to  determine  the  plasma 
frequency  of  the  metal,  should  have  general  applica¬ 
bility  in  the  IR.  Moreover,  to  produce  such  an  effect 
gratings  are  not  required,  since  the  interference  has 
already  been  detected  with  an  aperture-excitation 
technique8  on  a  smooth  metal  surface. 

Discussions  with  J.  J.  Quinn,  V.  A.  Yakovlev,  and  G. 
Zhizhin  by  A.  J.  Sievers  have  been  particularly  stimu¬ 
lating.  This  research  is  supported  by  the  National 
Science  Foundation  under  grant  DMR-84-09823  and 
by  the  U.S.  Air  Force  Office  of  Scientific  Research 
under  grant  AFOSR-85-0175.  Additional  support  has 
been  received  from  the  Materials  Science  Center 
(MSC)  at  Cornell  University,  MSC  Rep.  No.  5848. 
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The  technique  of  incoherent  laser  saturation  has  been  used  to  probe  the 
ls(Ai )  -*  ls(Ti )  transition  of  the  Se  donor  in  AlSb  at  liquid  helium 
temperatures.  The  intensity  dependence  of  the  absorptivity  is  described  by 
a  cascade  excitation  process  in  which  at  high  flux  two  photons  promote 
the  electron  into  the  conduction  band.  The  ls(T|)  lifetime  is  determined 
to  be  between  1  and  5  psec  which  is  about  three  orders  of  magnitude 


smaller  than  found  for  shallow  levels 

ALTHOUGH  HIGH  POWER  LASER  saturation  studies 
of  shallow  levels  in  semiconductors  show  that  impurity 
lifetimes  in  the  range  1  nsec  to  1  /jsec  can  be  obtained  at 
low  temperatures  [  1  -4] ,  no  similar  measurements  on 
deep  levels  have  been  reported.  Yet  incoherent  saturation 
measurements  on  such  transitions  should  be  quite 
informative  since  the  data  would  provide  some  evidence 
as  to  how  the  increase  in  electron-phonon  coupling 
constant  for  these  localized  states  is  balanced  against  the 
highly  forbidden  multiphonon  decay  required  to  relax 
them. 

In  this  paper  the  technique  of  incoherent  saturation 
[5J  has  been  used  to  investigate  the  dynamics  of  the 
ls(Ai)-»  1j(T()  transition  of  the  Se  donor  in  AlSb. 
This  transition  was  chosen  because  linear  spectroscopy 
[6]  shows  that  it  is  coincident  with  both  the  10P(18) 
and  10P(20)  COj  gas  laser  lines.  At  a  sample  tempera¬ 
ture  of  1.7  K  the  saturation  of  the  absorptivity  com¬ 
mences  at  a  beam  intensity  of  a  few  hundred  kWcm'2. 
Although  the  rate  of  change  of  absorptivity  with  incident 
intensity  is  much  faster  than  can  be  explained  by  a 
homogeneously  or  inhomogeneously  broadened  two 
level  system,  the  data  can  be  understood  with  a  three 
level  cascade  excitation  model.  The  measured  excited 
state  lifetime  is  ~  1-5  psec  which  is  about  103  times 
smaller  than  the  values  observed  for  shallow  defects 

[1-41 

The  sample  of  Se  doped  AlSb  is  from  the  batch 
described  in  [6] .  Our  spectroscopically  determined 
values  of  the  linewidth  and  center  frequency  of  the 
Is(Ai)-*  ls(T|)  transition  as  a  function  of  tempera¬ 
ture  are  shown  in  Fig.  1 .  At  low  temperature  the  tran¬ 
sition  which  occurs  at  944.71cm'1  is  straddled  by  two 
CO,  User  lines,  10P(18)  at  945.98cm*1  and  !0P(20) 
at  944.19  cm'1.  The  temperature  independence  of  the 
linewidth  at  low  temperatures  can  be  interpreted  in  one 


in  three-five  compounds. 

of  two  ways;  either  the  line  is  homogeneously  broadened 
and  the  low  temperature  linewidth  is  a  measure  of  the 
intrinsic  multiphonon  decay  lifetime  or  the  line  is 
inhomogeneously  broadened  due  to  crystal  strains  so  the 
low  temperature  line  width  obtained  from  linear  spectro¬ 
scopy  provides  no  intrinsic  information  about  the 
lifetime.  For  the  former  case  one  can  readily  estimate 
the  User  intensity  where  saturation  should  commence 
[5]  since 

ty  . 

I,  ~  ~  (^W),  0) 

where  w  is  the  center  frequency.  Aw  is  the  linewidth 
(FWHM)  and  a  is  the  absorption  cross  section.  For  this 
transition  the  linewidth  is  5.6  cm-1  and  a  —  2.1  x 
I0*,s  cm'2  so  /f  =  3.7  MW  cm'2.  Since  the  saturation 
parameter  will  be  smaller  for  the  case  of  inhomogeneous 
broadening  [5] ,  the  I,  calculated  above  represents  the 
maximum  possible  value  for  this  2  level  system. 

The  incoherent  saturation  measurements  have  been 
made  with  the  experimental  setup  shown  in  Fig.  2.  A 
line  tunable  C02  laser  Q-switched  at  a  frequency  of 
1 70  Hz  produces  200  nsec  pulses  with  peak  power  of 
1 .5  kW.  The  beam  is  focused  to  a  spot  size  of  9  x 
l  O'4  cm2  on  the  sample  which  is  immersed  in  pumped 
He  at  1 .7  K  in  an  optical  access  cryostat.  Sample  trans¬ 
mission  as  a  function  of  incident  laser  intensity  is 
obtained  by  moving  the  CaF2  calibrated  attenuators 
(labeled  A  in  Fig.  2)  one  by  one  from  one  side  of  the 
cryostat  to  the  other  so  that  errors  assocUted  with 
nonlinear  detector  response  are  minimized.  The  User 
output  is  monitored  with  a  pyroelectric  detector  which 
is  also  used  to  normalize  the  transmission  by  means  of 
a  two  channel  box  car  integrator. 

The  experimental  incoherent  saturation  results  are 
shown  in  Fig.  3.  The  data  in  the  top  picture  are  for 
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Fig.  1.  Temperature  dependence  of  the  line  width  and 
center  frequency  of  the  ls(Ai )  -*  lsfT, )  transition  in 
AlSb:Se.  The  value  of  the  low  intensity  absorption 
coefficient.  Oo,  at  line  center  is  4.2  cm'1 . 


Fig.  2.  Experimental  arrangement  for  incoherent  satura¬ 
tion  measurement.  Legend:  B,  beamsplitter;  A,  CaFj 
calibrated  attenuators;  D,  HgCdTe  detector;  R,  pyro¬ 
electric  reference  detector;  and  S,  sample  cryostat. 

10P(18)  and  the  bottom  picture  are  for  10P(20). 
Inspection  of  these  data  show  that  they  are  qualitat¬ 
ively  similar  but  quantitatively  different.  The  saturation 
intensity,  I„  taken  here  to  be  equal  to  the  incident 
intensity  value  at  which  the  absorption  coefficient  is 
equal  to  one  half  of  its  low  intensity  value,  is  1 20  kW 
cm'1  for  10P(  18)  and  400kWcm_l  for  10P(20).  The 
fact  that  the  two  curves  do  not  look  identical  may  be 
an  indication  of  heating  since  with  increasing  tempera¬ 
ture  the  absorption  line  center  frequency  in  Fig.  1  moves 
to  lower  frequencies  hence  away  from  the  1 0P(  1 8)  line 
and  towards  the  10P(20)  line.  Such  an  effect  does 
occur  in  the  data  shown  in  Fig.  3 ;  however,  the  exper¬ 
iment  described  next  indicates  that  it  is  electron  heating 
and  not  sample  heating  that  is  important.  By  running  the 
laser  in  a  CW-chopped  mode,  the  energy  per  pulse  at  the 
sample  can  be  increased  by  an  order  of  magnitude  over 


Fig.  3.  Absorption  coefficient  vs  laser  intensity.  Top: 
10P(18)  laser  line  on  the  high  frequency  side  of  the 
1j(Aj) -►  ls(Ti)  donor  line.  Bottom:  10P(20)  on  the 
low  frequency  side  of  the  donor  line.  The  sample  tem¬ 
perature  is  1 .7  K. 


hw 
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Fig.  4.  Schematic  diagram  of  the  cascade  excitation 
model. 

that  in  the  Q-switched  mode.  No  change  in  sample 
transmission  occurs  in  this  configuration  indicating  that 
sample  heating  is  not  important.  Another  unusual 
features  of  the  data  in  Fig.  3  is  that  the  rate  of  decrease 
in  absorptivity  with  increasing  incident  laser  intensity  is 
much  faster  than  can  be  explained  by  either  the  homo¬ 
geneous  or  inhomogeneous  two  level  system  (5].  Because 
of  the  electron  heating  contribution  to  the  intensity 
dependent  transmission  only  the  !0P(20)  data  which  is 
closest  to  line  center  will  be  analyzed  below  . 
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Both  the  rapid  intensity  dependence  of  the  absorpt¬ 
ivity  shown  in  Fig.  3  and  the  electron  heating  result  can 
be  understood  in  terms  of  a  three  level  cascade  excita¬ 
tion  model  which  promotes  carriers  into  the  conduction 
band.  This  idea  is  illustrated  in  Fig.  4.  At  high  laser 
intensity  the  first  Ixo  takes  the  carrier  from  the  ls(At) 
to  ls(T,)  and  the  second  fxo  takes  it  from  li(Tt )  to 
the  conduction  band.  If  the  relaxation  time  from  the 
conduction  band  back  to  the  bound  states  is  long 
enough  then  the  ls(A()-*  ls(Tt)  transition  will  appear 
to  saturate  at  high  intensities  but  the  saturation  intensity 
will  not  be  related  to  the  homogeneous  line  width 
obtained  from  linear  spectroscopy. 

To  estimate  the  functional  dependence  of  this  effect 
with  intensity  the  problem  is  simplified  to  three  levels.  If 
the  saturation  intensity  between  the  first  two  states  is 
called  It  and  the  parameter  for  the  next  two  states, 
Ih .  then  the  steady  state  solution  of  the  rate  equations 
for  the  three  level  model  gives  for  the  absorption 
coefficient 


a  =  a0/[l  +  (///,,  )+(///)*]. 


T  =  [4/., /«,]*'*. 


I_ 
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electron  to  go  from  the  conduction  band  to  the  ls(Ti) 
excited  state.  The  ls(Ti)  level  ionization  energy  [6 J  of 
29.4 meV  (237cm"1)  is  close  to  the  effective  mass 
theory  (EMT)  ground  state  binding  energy  for  AlSb  of 
40meV  (323  cm'1 )  so  the  EMT  value  [ 7 J  of  the  photo¬ 
ionization  cross  section,  o2  *5  x  10" 17  cm':.  is  used  in 
the  Ih  expression,  the  excited  state  analogue  of  equation 
(1).  Although  a  rather  long  relxxation  time  of  4  to  20 
nsec  is  obtained,  inspection  of  the  phonon  density-of- 
states  of  AlSb  indicates  that  the  phonon  gap  between 
the  optic  and  acoustic  phonon  branches  occurs  for  these 
transition  energies  [8] . 

In  conclusion,  incoherent  saturation  measurements 
on  the  ls(Ai  )-►  ls(Tj )  transition  show  that  the  satura¬ 
tion  intensity  is  orders  of  magnitude  larger  than  pre¬ 
viously  measured  for  shallow  states.  The  step  function 
turn  off  of  the  absorptivity  with  increasing  laser  intensity 
is  ascribed  to  a  two  step  excitation  process  which 
promotes  the  electron  into  the  conduction  band  coupled 
with  a  slow  relaxation  time  back  to  the  ls(T i )  level. 
This  process  masks  a  precise  determination  of  the 
ls(T!)  excited  state  lifetime  from  the  saturation  profile; 
however,  the  value  is  between  1  and  5  psec,  about  three 
orders  of  magnitude  smaller  than  found  for  shallow 
states  in  three-five  compounds  [4] . 


The  effective  saturation  intensity  in  the  high  intensity 
limit  is  the  geometric  mean  of  the  saturation  intensity 
for  each  of  the  transitions.  In  this  derivation  it  has  been 
assumed  that  rapid  relaxation  occurs  within  the  con¬ 
duction  band  and  also  that  the  time  for  relaxation  from 
the  conduction  band  directly  to  the  ground  state  is 
much  larger  than  any  of  the  other  times  in  the  problem. 
The  solid  curve  for  10P(20)  in  Fig.  3  is  obtained  with 
/,  =3.7MWcm'J  and  /  =  420kWcm"3.  The  reason- 
able  fit  to  the  data  demonstrates  that  the  cascade 
excitation  process  completely  masks  the  saturation 
behavior  of  the  IslA^-*  1j(T()  transition  so  it  is  not 
possible  from  this  incoherent  saturation  measurement  to 
determine  conclusively  whether  or  not  the  transition  is 
homogeneously  broadened;  however,  the  absorption 
linewidth  and  the  large  intensity  range  for  which  satura¬ 
tion  is  not  observed  to  occur  provide  definite  bounds  on 
the  ls(T, )  excited  state  lifetime.  The  value  is  between 
I  and  5  psec. 

The  measured  saturation  intensity  /  together  with 
/.  determined  from  the  linewidth  can  be  used  in 

•l 

equation  (3)  to  estimate  the  relaxation  time  for  the 
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